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Symmetry protected topological (SPT) phases are gapped short-range-entangled quantum phases 
with a symmetry G, which can all be smoothly connected to the same trivial product state if 
we break the symmetry. The Haldane phase of spin-1 chain is the first example of SPT phase 
which is protected by 50(3) spin rotation symmetry. The topological insulator is another exam- 
ple of SPT phase which is protected by U{1) and time reversal symmetries. It has been shown 
that free fermion SPT phases can be systematically described by the K-theory. In this paper, we 
show that interacting bosonic SPT phases can be systematically described by group cohomology 
theory: distinct d-dimensional bosonic SPT phases with on-site symmetry G (which may contain 
anti-unitary time reversal symmetry) can be labeled by the elements in 'H^^'^[G, Ut{^)\ - the Borel 
(1 -|- d)-group-cohomology classes of G over the G-moduIe ?7t(1). Our theory, which leads to ex- 
plicit ground state wave functions and commuting projector Hamiltonians, is based on a new type 
of topological term that generalizes the topological 6'-term in continuous non-linear cr-models to 
discrete non-linear a-models. The boundary excitations of the non-trivial SPT phases are described 
by continuous/discrete non-linear cr-models with a non-local Lagrangian term that generalizes the 
Wess-Zumino-Witten term for continuous non-linear cr-models. As a result, the symmetry G must 
be realized as a non-on-site symmetry for the low energy boundary excitations, and those boundary 
excitations should be gapless if the symmetry is not broken on the boundary. As an application of our 
result, we can use 'H^^'*[(7(l) xi Z2, (7t(1)] to obtain interacting bosonic topological insulators (pro- 
tected by time reversal Z2 and boson number conservation), which contain one non-trivial phases 
in ID or 2D, and three in 3D. We also obtain interacting bosonic topological superconductors (pro- 
tected by time reversal symmetry only), in term of 'H^^''[Z2 , I7t(1)], which contain one non-trivial 
phase in odd spatial dimensions and none for even. Our result is much more general than the above 
two examples, since it is for any symmetry group. For example, we can use 'H^'^''[U{1) x Zj , Urii)] 
to construct the SPT phases of integer spin systems with time reversal and 17(1) symmetry, which 
contain three non-trivial SPT phases in ID, none in 2D, and seven in 3D. Even more generally, 
we find that the different bosonic symmetry breaking short-range-entangled phases are labeled by 

the following three mathematical objects: (Gh,G*,^i+'*[G*,C/t(1)]) , where Gh is the symmetry 
group of the Hamiltonian and G* the symmetry group of the ground states. 

PACS numbers: 71.27.+a, 02.40.Re 
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I. INTRODUCTION 
A. Background 

Understanding phases of matter is one of the central 
issues in condensed matter physics. For a long time, we 
believed that all the phases and phases transitions were 
described by Landau symmetry breaking theory. ^^'^ In 
1989, it was realized that many quantum phases can con- 
tain a new kind of orders which are beyond the Landau 
symmetry breaking theory. A quantitative theory of the 
new orders was developed based on robust ground state 
degeneracy and the robust non-Abelian Berry's phases 
of the degenerate ground states, which can be viewed 
as new "topological non-local order parameters" .'^^^ The 
new orders were named topological order. Topologically 
ordered states contain gapless edge excitations and/or de- 
generate sectors that encode all the information of bulk 
topological orders. ^'^ The nontrivial edge states provide 
us a practical way to experimentally probe topological 
order and illustrate the holographic principle which was 
introduced later. The excitations in those topologi- 
cally ordered states in general carry fractional charges^ ^ 
and obey fractional statistics. ^^"^^ 

Since its discovery, we have been trying to obtain 
a systematic and deeper understanding of topological 
orders. Recently, it was realized that topological or- 
ders are related to long range entanglements. ^^'^^ In 
fact, we can regard topological orders as patterns of 
long range entanglements^* defined through local unitary 
(LU) transformations. ^^"^^ 

The notion of topological orders and long range entan- 
glements leads to the following more general and more 
systematic picture of phases and phase transitions (see 
Fig. 1).^* For gapped quantum systems without any 
symmetry, their quantum phases can be divided into two 
classes: short range entangled (SRE) states and long 
range entangled (LRE) states. 
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SRE states are states that can be transformed into 
direct product states via LU transformations. All SRE 
states can be transformed into each other via LU trans- 
formations. So all SRE states belong to the same phase 
(see Fig. la). 

LRE states are states that cannot be transformed into 
direct product states via LU transformations. It turns 
out that, many LRE states also cannot be transformed 
into each other. The LRE states that are not con- 
nected via LU transformations belong to different classes 
and represent different quantum phases. Those differ- 
ent quantum phases are nothing but the topologically or- 
dered phases. Fractional quantum Hall states^^'^^, chiral 
spin liquids, ^"^'^^ spin liquids, ^^"^^ non-Abelian frac- 
tional quantum Hall states, ^^"'^^ etc are examples of topo- 
logically ordered phases. The mathematical foundation 
of topological orders is closely related to tensor category 
theory^^'^'^'"^'^''^* and simple current algebra. ^^''^^ 

For gapped quantum systems with symmetry, the 
structure of phase diagram is even richer (see Fig. lb). 
Even SRE states now can belong to different phases. The 
Landau symmetry breaking states belong to this class 
of phases. However, there are more interesting exam- 
ples in this class. Even SRE states that do not break 
any symmetry and have the same symmetry can belong 
to different phases. The ID Haldane phases for spin- 
1 chain'^^''^^ and topological insulators'^*"'*'^ are examples 
of non-trivial SRE phases that do not break any symme- 
try. Those phases are beyond Landau symmetry breaking 
theory since they do not break any symmetry. We will 
call those phases Symmetry Protected Topological (SPT) 
phases. 

For gapped quantum systems with symmetry, the cor- 
responding LRE phases will be much richer than those 
without symmetry. We may call those phases Symme- 
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try Enriched Topological (SET) phases. Projective sym- 
metry group (PSG) was introduced to study the SET 
phases. ^^'"'^ Many examples of this kind of states can be 
found in Ref. 44,46-48,50, but a systematic understand- 
ing is still lacking. 



B. Motivation 

The notion of topological order and long range entan- 
glements deepens our understanding of quantum phases 
and guides our research strategy. This allows us to make 
significant progress. 

For example, there is no long range entanglement in 
gapped ID states. ^^'^^ So, without symmetry, all gapped 
ID quantum states belong to the same phase. For sys- 
tems with a certain symmetry, all gapped ID phases are 
either SPT phases protected by symmetry or symmetry 
breaking states. Since both SPT phases and symmetry 
breaking states are short range entangled, it is easy to 
understand them. As a result, a complete classification 
of all ID gapped bosonic/fermionic quantum phases for 
any symmetry can be obtained. ^*"''''* (A special case of 
the above result, a classification of ID fermionic systems 
with T"^ — \ time reversal symmetry can also be found 
in Ref. 53,55,56). Using the idea of LU transformations, 
we also developed a systematic and quantitative theory 
for non-chiral topological orders in 2D interacting boson 
and fermion systems. **'^'''^'' We would like to mention 
that symmetry protected Berry phases have been used 
to study various topological phases. ^^'^* 

Motivated by the ID classification result, in this pa- 
per and in Ref. 59 we would like to study SPT phases in 
higher dimensions. Since SPT phases are short range en- 
tangled, it is relatively easy to obtain a systematic under- 
standing. (Another way to make the classification prob- 
lem easier is to consider only free fermion systems which 
are classified by K-theory.^°'^^) In Ref. 59, we study some 
simple but highly non-trivial examples. Those non-trivial 
examples lead to the generic and systematic results dis- 
cussed in this paper. Some other examples of 2D gapped 
SPT phases are given in Ref. 48,49,51,52,62. 



C. Summary of results 



FIG. 1: (Color online) (a) The possible gapped phases for 
a class of Hamiltonians H(cj\,g2) without any symmetry re- 
striction, (b) The possible gapped phases for the class of 
Hamiltonians -ffsymm(3i, 32) with symmetry. Each phase is 
labeled by its entanglement properties and symmetry break- 
ing properties. SRE stands for short range entanglement, 
LRE for long range entanglement, SB for symmetry break- 
ing, SY for no symmetry breaking. SB-SRE phases are the 
Landau symmetry breaking phases, which are understood by 
introducing group theory. The SY-SRE phases are the SPT 
phases, and we will show that they can be understood by 
introducing group cohomology theory. 



Using group theory, we can obtain a systematic un- 
derstanding of symmetry breaking phases (or more pre- 
cisely, short-range entangled symmetry breaking phases) . 
In this paper, we will show that, using group co- 
homology theory, we can obtain a systematic under- 
standing of short-range entangled symmetric phases of 
bosons/qubits, even with strong interactions. Those 
phases are called bosonic SPT phases. In particular, we 
have obtained the following results for bosonic systems: 

1. From each element in (1 -f (i)-Borel-cohomology 
group ■H^^''[G', C/t(1)], we can construct a dis- 



4 



tinct SPT phase that respects the on-site sym- 
metry G in d-spatial dimensions. Here G may 
contain anti-unitary time reversal transformation 
and H^+'^[G, C/t(1)] is introduced in appendix D. 
Note that 'H}+'^[G, Ut{1)] is an Abehan group that 
can be calculated from the symmetry group G. 
The identity element in 'H^^'^[G^ Ut{^)] correspond 
to trivial SPT phases while other elements corre- 
sponds to non-trivial SPT phases. For example, 
'H'^[SO{-i),U(l)] = Z2. So a ID integer spin chain 
with SO(2>) spin rotation symmetry (but no trans- 
lation symmetry) has two kinds of SPT phases: one 
is the trivial 5* = phase and the other is the Hal- 
dane phase.^^'^'^ 

2. The low energy effective theory of a SPT phase with 
symmetry G is given by a topological non-linear tr- 
model that contains only a 27r-quantized topologi- 
cal 9-ieim. The 27r-quantized topological 6'-term in 
(d-h 1)D is classified by 1-L^+'^[G, Ut{1)] which gen- 
eralizes the topological term for non-linear cr-model 
with continuous symmetry. 

3. We argue that a d-dimensional non-trivial SPT 
phase has gapless boundary excitations if the on- 
site symmetry is not broken at the boundary. This 
is similar to the holographic principle for intrin- 
sic topological orders.^'* The boundary excitations 
of the SPT phase are described by a non-linear 
cr-model with a non-local Lagrangian (NLL) term 
that generalizes the Wess-Zumino-Witten (WZW) 
term^'^'®'^ for continuous non-linear cr- models. In 
(1+1)D and for continuous symmetry group, it is 
shown that non-linear cr-model with WZW term is 
gapless which is described by Kac-Moody current 
algebra. We generalize the gaplessness result to 
any symmetry groups and to any dimensions, pro- 
vided that the symmetry is not broken. 

4. The SPT phases that respect on-site symmetry 
G and translation symmetry can be obtained in 
the following way: In 1-dimension, those phases 
are labeled by H^[G, C/t(1)] x ^^[G, Ut{1)].^^-'^'^ In 
2-dimension, they are labeled by 'H^[G,Ut{^)] x 

|-H2[G,C/t(1)]} X ■H3[G,C/j,(1)]. a partial re- 
sult ■HHG,L/t(1)] X |-H2[G,[/t(1)]}^ was ob- 
tained in Ref. 51. In 3-dimension, they are 

labeled by n^lG^Uril)] x {-H^iG, [/t(1)]} x 
{H3[G,C/T(l)]}'xH4[G,t/T(l)]. 

This paper is organized as the following: In section 
II, we list the bosonic SPT phases for many symmetry 
groups in dimension 0,1,2,3, and discuss some examples 
of those SPT phases. In section III, we give a brief re- 
view of local unitary transformations. In section IV, we 
discuss canonical form of the ground state wave function 



for SPT phases. In section V, we study on-site sym- 
metry transformations that leave the canonical ground 
state wave function unchanged. In section VI, we con- 
struct the on-site symmetry transformations through the 
cocycles of the symmetry group. In section VII, we intro- 
duced topological non-linear cr-model and discuss their 
SPT phases. We also argue that the boundary states 
of the topological non-linear cr-model are gapless if the 
symmetry is not broken. In section VIII, we construct 
and classify topological non-linear cr-model through the 
cocycles of the symmetry group. In section IX and X, we 
show that the ground states of the topological non-linear 
CT-model all have trivial intrinsic topological orders, and 
the same SPT order if constructed from equivalent cocy- 
cles. In section XI, we discuss the relation between the 
cocycles in the topological non-linear cr-model and the 
Berry's phase. In section XII, we study SPT phases with 
both on-site and translation symmetries. 

II. EXAMPLES OF BOSONIC SPT PHASES 

In table I, we list the SPT phases for some simple sym- 
metry groups. In the following, we discuss some of those 
phases in detail. We also give some simple examples for 
some of the listed SPT states. 



A. SO(3) SPT states 

For integer spin systems with the full S'0(3) spin ro- 
tation symmetries, the symmetry group is 5*0(3). From 
H^+'^[SO{'i), C/(l)], we find one non-trivial SPT phase in 
ID and infinite many in 2D. Those 2D SPT phases la- 
beled by fc G Z have a special property that they can 
be described by continuous non-linear cr-model with 27r- 
quantized topological 0-term: 

S = j drd^a: (^Tr(a^.9ta,,g) (1) 
+ ^ 2^ ^Tr[(.g-ia,g) {g^'d^g) {g~'d>.g)]) , 

where g{x,t) is a 3 x 3 matrix in 50(3) and 9 — 27rfc, 
fc e Z. This is because the topological term 

I drd'x -^'-^^Tr[{g-%g){g-^d,g){g-'d,g)] 

(2) 

corresponds to H^[SO{3), U{1)] whose elements are also 
labeled by fc £ Z. 

At the boundary, the topological term reduces to the 
well known WZW term which gives rise to gapless edge 
excitations.^'* We would like to point out that g~^{d^g) 
and {dfig)g~^ create excitations on the edge that move in 
the opposite directions. Since the 50(3) symmetry acts 
as g{x,t) — )■ hg{x,t), h e S'0(3), only {d^g)g~^ carries 
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TABLE I: (Color online) SPT phases of interacting bosonic systems in d-spatial dimensions protected by on-site symmetry 
G. In absence of translation symmetry, the above table lists 'H^^'^fG, (/t(1)] whose elements label the SPT phases. Here Zi 
means that our construction only gives rise to the trivial phase. Z„ means that the constructed non-trivial SPT phases plus the 
trivial phase are labeled by the elements in Z„. Z^ represents time reversal symmetry, "trn" represents translation symmetry, 
f/(l) represents f/(l) symmetry, Z„ represents cyclic symmetry, etc . Also (m, n) is the greatest common divisor of m and n. 
The red rows are for bosonic topological insulators and the blue rows bosonic topological superconductors. The red/blue rows 
without translation symmetry correspond to strong bosonic topological insulators/superconductors and the red/blue rows with 
translation symmetry also contain weak bosonic topological insulators/superconductors. 



non trivial 5*0(3) quantum numbers while g^^{d^g) is a 
50(3) singlet. So on the edge, the excitations with non- 
trivial 50(3) quantum numbers all move in one direction, 
which breaks the time reversal and parity symmetry. In 
fact, under time reversal or parity transformations, 9 — > 
—9 and k — >■ —k. 

In the above example, we see that a 27r-quantized topo- 
logical 0-terms in a non-linear cr-model gives rise to a non- 
trivial SPT phase. However, the topological 6'-terms in 
continuous non-linear cr-model do not always correspond 
to non-trivial SPT phases. For example, 7r2(50(3)) — 
and the continuous 50(3) non-linear cr-model has no 
topological 0-terms in (l-l-l)D. But we do have a ID non- 
trivial SPT phase protected by 50(3) symmetry. Also, 
7r4(50(3)) = Z2 and the continuous 50(3) non-linear cr- 
model has non-trivial topological 0-terni in (3-|-l)D. But 
such topological 0-term cannot produce non-trivial SPT 
phase protected by 50(3) symmetry, since the topologi- 
cal term becomes trivial once we include the cut-off. In 
this paper, we will show that non-trivial 27r-quantized 
topological 0-terms can even be defined for discrete non- 
linear CT-models on discretized space-time. It is the 2tt- 
quantized topological 0-terms in discrete non-linear a- 
models that give rise to non-trivial SPT phases. 



It is also possible that the above 50(3) SPT phases 
labeled by k and fc -I- 1 are connected by a continuous 
phase transition that do not break any symmetry. The 
gapless critical point is likely to be described by eqn. (1) 
with6l = 27r(fc-h5). When 6* < 27r(fc-H i), it may flow to 
27rfc at low energies and when 9 > 27r(fc + |), it may flow 
to 2TT{k + 1). Since all the SPT phases are described by 
27r-quantized topological 0-terms in discrete non-linear a- 
modcls, the above picture about the transitions between 
SPT phases may be valid for generic SPT phases. 



For integer spin systems with time reversal and the full 
50(3) spin rotation symmetries, the symmetry group is 
50(3) X Zj. From n^+'^[SOi3) x Z^,Ut{1)], we find 
one non-trivial SPT phase in 2D and seven in 3D. Note 
that on systems with boundary, the topological 0-term in 
eqn. (1) breaks the time reversal symmetry. So we cannot 
use those Z classified topological 0-terms to produce 2D 
SPT phases with time reversal and 50(3) spin rotation 
symmetries. As a result, the 2D SPT phases with time 
reversal and 50(3) spin rotation symmetries are only 
described by Z2. 
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B. SU{2) SPT states 

For bosonic systems with SU{2) symmetry, the SPT 
phases are labeled by n'^+'^[SU{2), U{1)]. We find infinite 
many non-trivial SU{2) SPT phases in (2 + 4n) spatial 
dimension. Those SU{2) SPT phases labeled by A; € 
Z. There is no non-trivial SU{2) SPT phase in other 
dimensions. Similarly, those SU{2) SPT phases in 2- 
dimensions can be described by continuous non-linear a- 
model with 27r-quantized topological 0-term: 

j ATA^x{^Tv{d^g^d^9) (3) 

+ i ^ ^ \^A{9~'d,g){g-'d,g){g-'d>.g)]) , 
where g{x,t) is a 2 x 2 matrix in SU{2) and 9 — 27rfc, 

fc e z. 



did not have the U{1) symmetry, since we could include 
one of the corresponding quasiparticle operators in the 
Hamiltonian which would gap the edge excitations.''^ In 
the presence of U{1) symmetry, the quasiparticles carry 
U{1) charges q^K~^ni — 1 — k and q^K~^n2 — 1. We 
see that when k ^ 1, both quasiparticles that correspond 
to the null vectors carry non-zero U{1) charges. Thus, the 
quasiparticle operators cannot be included in the Hamil- 
tonian, and they do not gap the gapless edge excitations. 
The correspond state will have U{1) protected gapless ex- 
citation and correspond to a non-trivial U{1) SPT state. 
We see that the X-matrix and the charge vector q de- 
scribe a non-trivial U{1) SPT state when fc ^ 1 and a 
trivial state when fc = 1. The 2D U{1) SPT states are 
labeled by an integer. 



D. Bosonic topological insulators/superconductors 



C. [/(I) SPT states 

From H^+'^[U{1),U{1)] = Z for even d and 
n^+'^[U{l),Uil)] = Zi for odd d, we find that spin/boson 
systems with U{1) on-site symmetry have infinite non- 
trivial SPT phases labeled by non-zero integer in d = 
even dimensions. This generalizes a result obtained by 
Levin for d = 2.^9 We note that n^[SU{2),U{l)] = 
n^[U{l),U{l)] = Z. The SPT states with SU{2.) sym- 
metry can also be viewed as SPT states with U{1) sym- 
metry. We know that an SU{2) SPT state labeled by 
fc € Z is described by eqn. (3) with 9 = 2i:k. One may 
wonder if such an SU{2) SPT state is also an non-trivial 
U{1) SPT state labeled by /c e Z. 

We like to point out that it is believed that all 2D 
gapped phases with Abelian statistics are classified by K- 
matrix and the related U{1) Chern-Simons theory.^'^"^^ 
All the quasiparticles in the 2D SPT phases are bosons. 
So the SPT phases are also described by i4'-matrices. We 
just need to find a way to include symmetry in the K- 
matrix approach, which is done in Ref. 48. In particular, 
Michael Levin^^ pointed out that a 2D C/(l) SPT phase 
can be described by a U{1) x U{1) Chern-Simons theory 
(or a double-layer quantum Hall state) (see also Ref. 69, 
70) 

C = ^Kjjai,,d,aj^e'^-'^ + ^qjA^d^aj^e^''^ + ... (4) 
An 2n 



with the if-matrix and the charge vector q 



,65-67 



K = 



1 

1 2k 



Q = 



(5) 



We note that such a i^T-matrix has two null vectors 
rii = ,n2 = that satisfy nfK^^rii = 0. The 

null vectors correspond to quasiparticles with Bose statis- 
tics. Such null vectors would destabilize the state if we 



The C/(l) X fine in Table I describes the SPT 
phases for interacting bosons with time reversal symme- 
try Z2 and boson number conservation (symmetry group 
— t/(l) X Z2 , where time reversal T and U{1) transforma- 
tions Ug satisfy TUg = U-gT). Those phases are bosonic 
analogues of free fermion topological insulators protected 
by the same symmetry. From H^^'^[U{1) x Zj, U{1)], we 
find one kind of non-trivial bosonic topological insulators 
in ID or 2D, and three kinds in 3D. The only non-trivial 
topological insulator in ID is the same as the Haldane 
phase. 

The Z2 line in Table I describes interact- 
ing bosonic analogues of free fermion topological 
superconductors''^"'^^ with only time reversal symme- 
try, . Since n^+''[Z^ ,U{1)] = Z2 for odd d and 
n^+'^[Z^,U{l)] = Zi for even d, we find one kind of 
"bosonic topological superconductors" or non-trivial 
SPT phases in every odd dimensions (for the spin/boson 
systems with only time reversal symmetry). 



E. Other SPT states 

The U{l)x Z2 line describes the SPT phases for integer 
spin systems with time reversal and U{1) spin rotation 
symmetries (symmetry group — U{1) x Zj, where time 
reversal T and U{1) transformations Ug satisfy TUg — 
UgT). From %^+'^[U{l) x Zj,[/(1)], we find three non- 
trivial SPT phases in ID, none in 2D, and seven in 3D. 

We also find that '}V-+'^[Zn, U{1)] = Z„ for even d and 
'H^+'^[Zn,U{l)] = Zi for odd d. So spin/boson systems 
with Z„ on-site symmetry have n—\ kinds of non-trivial 
SPT phases in d = even dimensions. 

For integer spin systems with D2h symmetry but no 
translation symmetry, we discover 15 new SPT phases 
in 1D,53,77 g3 gprp phases in 2D, and 511 new SPT 
phases in 3D. 
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FIG. 2; (a) A triangular lattice. The Hamiltonian term (6) 
acts on the seven sites in the shaded area, (b) A geometric 
representation of the the phase factors in eqn. (6). 



F. Ideal ground state wave functions and exactly 
soluble Hamiltonians for SPT phases 

We can construct the idea ground state wave func- 
tions and exactly soluble Hamiltonians for all the SPT 
phases described by 'H^^'^[G,Ut{1)]- The elements in 
H^^'^[G,Ut{1)] are complex functions of d + 2 vari- 
ables iyd+iigo,--,9d+i), 9i e G. Vd+i{9o,--,9d+i) 
is a pure phase \v d+i{gQ gd+i)\ = 1 that sat- 
isfy certain cocycle conditions (IC) and (17). From 
each element Vd+i{9oi ■■■T9d+i) we can construct the d- 
dimensional ground state wave function for the corre- 
sponding SPT phase. In 2D, we can start with a tri- 
angle lattice model where the physical states on site- 
i are given by \gi), gi € G (see Fig. 2a). The ideal 
ground state wave function is then given by ^{{gi\) = 
]\A^3i)-:9i.g3,9k)]\yV^^{y,gt,g3,gk), where Ha and 
multiply over all up- and down-triangles, and the or- 
der of ijk is clockwise for up-triangles and anti-clockwise 
for down-triangles (see Fig. 2a). 

To construct exactly soluble Hamiltonian H that re- 
alizes the above wave function as the ground state, 
we start with an exactly soluble Hamiltonian Hq = 
-Y.^\^^){'l^^l = EgiGcl^i), whosc grouud statc 

is i'o({5i}) = 1. Then, using the local unitary trans- 
formation U = Y{^Vi{l,gi,gj,gk)Y{^v^^{i,gi,gj,gk), 
we find that the above ideal ground state wave func- 
tion is given by $ — [/$o and the corresponding ex- 
actly soluble Hamiltonian is given hy H = Hi , where 
Hi = U\(f)i){(f)i\W . Hi acts on a seven-spin cluster la- 
beled by i, 1 - 6 in shaded area in Fig. 2a 

Hi\gi,gig2g3g4g5ge) = ^ I5i,5i525354ff556)x 

Mg4, gs, g», g'd'^ajgb, g», g6)t^3(g», g[, ge, gi) 
i^3(gi,gj,g2,gi)j^3(g3,g»,g^,g2)i'3(g4,g3,gi,gD 

The above phase factor has a graphic representation as in 
Fig. 2b. (For a detailed explanation of the graphic repre- 
sentation see Fig. 10.) H has a short ranged interaction 
and has the symmetry G: \{gi}) — > \{ggi}}, g G G, if 
i^sigo, gs) satisfies the 3-cocycle conditions eqn. (IG) 
and eqn. (20). 

For symmetry G = Z2 and using the 3-cocycle calcu- 



lated in section J 2, we find that the Hamiltonian that 
realize the non-trivial Z2 SPT state in two dimensions is 
given by 

(7) 

where 













or 


^ (0 





which act on site-i. Also rf^- are operators acting on site-i 
and site-j: 

^+|a),®|/3), = |a),®|/3)„ (a,/3)^(0,l) 

^^^■|i)^^|o), = -|i).®|o)„ 

V~JW)^®\P), = W)^<^\P),. (a,/3) 7^ (1,0). (9) 

G. A classification of short-range-entangled states 
with or without symmetry breaking 

The above results are for bosonic states that do not 
break any symmetry of the Hamiltonian. Combin- 
ing group theory (that describe the symmetry break 
states) and group cohomology theory (that describes the 
SPT states), we can obtain a theory for more general 
short-range-entangled states that may break the sym- 
metry Gh of the Hamiltonian down to the symmetry 
G^ of the ground states. We find that the different 
symmetry breaking short-range-entangled phases are de- 
scribed/labeled by the following three mathematical ob- 
jects: (Gh,G*,H1+''[G*,[/t(1)]). 

Landau symmetry breaking theory tries to use 
(Gif,G>i.) to describe/label all the symmetry break- 
ing short-range-entangled phases. We see that 
Landau symmetry breaking theory misses the 
third label H^+'^[G^ ,Ut{1)]. The SPT phases 
do not break any symmetry and are described by 

(Gh. Gh,H^+''[Gh, Ut{1)]) U^+^IGh. Ut{1)]. 



III. LOCAL UNITARY TRANSFORMATIONS 

In the rest of the paper, we will explain the ideas, the 
way of thinking, and the detailed calculations that allow 
us to obtain the results described in the above section. 
We will start with a short review of local unitary (LU) 
transformation . 1 1 

LU transformation is an important concept which is 
directly related to the definition of quantum phases.^® In 
this section, we will explain what it is. Let us first intro- 
duce local unitary evolution. A LU evolution is defined 



1 2 



(b) 



(a) 



FIG. 3: (Color online) (a) A graphic representation of a 
quantum circuit, which is form by (b) unitary operations on 
blocks of finite size I. The green shading represents a causal 
structure. 



as the following unitary operator that act on the degrees 
of freedom in a quantum system: 



(10) 



where T is the path-ordering operator and H{g) — 
^iOi{g) is a sum of local Hermitian operators. Two 
gapped quantum states belong to the same phase if and 
only if they are related by a LU evolution. ^^'''^'^^ 

The LU evolutions is closely related to quantum cir- 
cuits with finite depth. To define quantum circuits, let us 
introduce piecewise local unitary operators. A piecewise 
local unitary operator has a form 



pwl 



where {U^} is a set of unitary operators that act on non 
overlapping regions. The size of each region is less than 
some finite number I. The unitary operator Upwi defined 
in this way is called a piecewise local unitary operator 
with range I. A quantum circuit with depth M is given 
by the product of M piecewise local unitary operators: 



r/^^ = f/^^' f/(^^ ■■■u 



(M) 
''pwl 



We will call J7*^c ^ transformation. In quantum in- 
formation theory, it is known that finite time unitary 
evolution with local Hamiltonian (LU evolution defined 
above) can be simulated with constant depth quantum 
circuit (?e a LU transformation) and vice-verse: 



M 

circ 



(11) 



So two gapped quantum states belong to the same phase 
if and only if they are related by a LU transformation. 

In this paper, we will use the LU transformations to 
simplify gapped quantum states within the same phase. 
This allows us to gain a deeper understanding and even 
to classify gapped quantum phases. 



IV. CANONICAL FORM OF MANY-BODY 
STATES WITH SHORT RANGE 
ENTANGLEMENTS 

A generic many-body wave function $(mi, to^v) is 
very complicated. It is hard to see and identify the quan- 



(a) 



(b) 



FIG. 4: (Color online) Transforming a SRE state to a tensor- 
network state which take simple canonical form, (a) A SRE 
state, (b) Using the unitary transformations that act within 
each block, we can transform the SRE state to a tensor- 
network state. Entanglements exist only between the degrees 
of freedom on the connected tensors. 



m 



(a) 



a 



m 



(b) 



m 
X 



.a 



(c) 



(d) 



FIG. 5: (Color online) Graphic representations of tensors: 
(a) yl", (b) A"^, and (c) A^^^.,. (d) A corner represents a 
special rank-2 tensor Aap = 5^/3 ■ 



tum phase represented by a many-body wave function. In 
this section, we will use LU transformations to simplify 
many-body wave functions in order to understand the 
structure of quantum phases. 

Such an approach is very effective in ID^^^'"''^ which 
leads to a complete classification of gapped ID phases. 
In two dimensions, the approach allows us to classify non- 
chiral topological orders. In this paper, we will 
study another problem where such an approach is ef- 
fective. We will use LU transformations to study SRE 
quantum phases with symmetries and study SPT phases 
that do not break any symmetry. 



A. Cases without any symmetry 

Without any symmetry, we can always use LU trans- 
formations to transform a SRE wave function into a prod- 
uct state. In the following, we will describe how to choose 
such LU transformation and what is the form of the re- 
sulting product state. 

We first divide our system into patches of size I as in 
Fig. 4a. If / is large enough, entanglement only exists be- 
tween regions that share an edge or a corner. In this case, 
we can use LU transformation to transform the state in 
Fig. 4a into a state with many unentangled regions (see 
Fig. 4b). For example, some degrees of freedom in the 
middle square in Fig. 4a may be entangled with the de- 
grees of freedom in the three squares below, to the right, 
and to the lower-right of the middle square. We can use 
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the LU transformation inside the middle square to move 
all those degrees of freedom to the lower-right corner of 
the middle square. Similarly, we can use the LU trans- 
formation to move all the degrees of freedom that are 
entangled with the three squares below, to the left, and 
to the lower-left of the middle square to the lower-left 
corner of the middle square, etc. Repeat such operation 
to every square and we obtain a state described by Fig. 
4b. For stabilizer states, such reduction procedure has 
been established explicitly.^" 

Fig. 4b is a graphic representation of a tensor-network 
description of the state. In the graphic representa- 
tion, a dot with n legs represents a rank n tensor (see 
Fig. 5). If two legs are connected, the indexes on those 
legs will take the same value and are summed over. In 
the tensor-network representation of states, we can see 
the entanglement structure. The disconnected parts of 
tensor-network are not entangled. In particular, the 
tensor-network state Fig. 4b is a direct product state. 

If there is no symmetry, we can transform any direct 
product state to any other direct product state via LU 
transformations. So all SRE states belong to one phase. 



B. Cases with an on-site symmetry 

However, when we study phases of systems with cer- 
tain symmetry, we can only use the LU transformations 
that respect the symmetry to connect states within the 
same phase. In this case, even SRE states with the same 
symmetry can belong to different phases. 

Let us consider d-dimensional systems of N sites that 
have only an on-site symmetry group G. We also assume 
that the states \m) on each site form a linear representa- 
tion Umm'ig), 5 G G of the group G. 

To understand the structure of quantum phases of the 
symmetric states that do not break the symmetry G, we 
can only use symmetric LU transformation that respects 
the on-site symmetry G to define phases. Two gapped 
symmetric states are in the same phase if and only if they 
can be connected by a symmetric LU transformation.-'^^ 

We have argued that generic LU transformations can 
change a SRE state in Fig. 4a to a tensor-network state 
in Fig. 4b. The LU transformations rearrange the spa- 
tial distributions of the entanglements which should not 
be affected by the on-site symmetry G. So, in the fol- 
lowing, we would like to argue that symmetric LU trans- 
formations can still change a SPT state in Fig. 4a to a 
symmetric tensor-network state in Fig. 4b (although a 
generic proof is missing) . 

We first assume that symmetric SRE states have tensor 
network representation as shown in Fig. 6. The linked 
dots represent the entangled degrees of freedom. The 
dots in each shaded circle represent a site, which forms 
a linear representation of the on-site symmetry group G. 
We then divide the systems into large squares (see Fig. 
6). The size of the square is large enough such that en- 
tanglement only appears between squares that share an 




FIG. 6: (Color online) A tensor network representation of a 
SRE state with on-site symmetry G. The all the dots in each 
shaded circle form a site. The degrees of freedom on each site 
(ie in each shaded circle) form a hnear representation of G. 
However, the degrees of freedom on each dot may not form a 
hnear representation of G. 



FIG. 7: (Color online) The canonical tensor network repre- 
sentation for ID SRE state I'&psre). The two dots in each 
rectangle represent a physical site. 



edge or a vertex. Now we view the degrees of freedom 
in each square as a large effective site. The degrees of 
freedom on each effective site form a linear representa- 
tion of G. Now, we can use an unitary transformation in 
each square to rearrange the degrees of freedom in that 
square (which corresponds to change basis in the large 
effective site). This way, we can transform the SPT state 
in Fig. 6 into the canonical form in Fig. 4b, where the 
degrees of freedom on each shaded square form a linear 
representation of G. So Fig. 4b is a symmetric tensor- 
network state. We would like to point out that although 
in Fig. 4b, we only present a 2D tensor-network state in 
canonical form, the similar reduction can be done in any 
dimensions. 



V. CLASSIFY SYMMETRY 
TRANSFORMATIONS OF SPT STATES 

After the symmetric state being reduced to the canon- 
ical form in Fig. 4b, the on-site symmetry transforma- 
tion is generated by the following matrix on the effective 
site-i: [/* , , , , which forms a linear represen- 

tation of the on-site symmetry group G. The symmetry 
transformation , , / / keeps the SRE state 

aia20t3a4,a^a20t^a^ ^ 
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a- 



A, 4 "^3 M 



"4 Pa 



FIG. 8: (Color online) The canonical tensor network repre- 
sentation for 2D SRE state I^I/pSre)- The four dots in each 
square represent a physical site. 



I^pSRe) in Fig. 7 or Fig. 8 invariant: 

«)iC/i*pSRE) = I^'pSRe) 



(12) 



for any lattice size. 

Eqn. (12) is one of the key equations. It describes 
the condition that the on-site symmetry transformations 
t/* must satisfy so that those on-site symmetry trans- 
formations can represent the symmetry of a SRE state. 
So to classify all possible symmetry transformations of 
SPT states, we need to find all the pairs (f/*, I^'psre)) 
that satisfy eqn. (12). Those different solutions can cor- 
respond to different SRE symmetric phases. 

However, two different solutions may not corre- 
spond to different phases. They may be "equivalent" 
and can correspond to the same phase. So to understand 
the structure of SRE symmetric phases, we also need to 
find out those "equivalent" relations. Clearly one "equiv- 
alent" relation is generated by unitary transformations 



ai 020304, /3i/32/33 ^4 



on each effective physical site: 



V , , , , 



where the repeated indices are summed over. Here 
^ ^ „ a R a a are not the most general on-site uni- 

010:20:304, P1P2P3P4 o 

tary transformations. They are the on-site unitary trans- 
formations that map IvE'pSre) to another state I^psre) 
having the same form as described by Fig. 7 or Fig. 8. 
The second "equivalent" relation is given by 




u. 



axa20304,o\o'^o'^o\ ^ ai/3iQ2&a3&a404,ai/3jQ2/'2"3^3"4' 

- ^ 0^020304,0'^0'^0'^0'^ ^^PuPi h,l3'2 Pz,^'^ Pa^P'^ 

where W^''^^, , a — 1,2,3,4 are linear representations of 
the on-site symmetry group G which satisfy a condition 
that the direct product representation VF^'* (g) H^^^i+a; ^ 
^/3,i+x+y ^ y^i,i+y contains a trivial ID representation: 



W 



w 



2,i+x 



(15) 



FIG. 9: (Color online) Adding four local degrees of freedom 
that form a ID representation does not change the phase of 
state. 



Such an "equivalent" relation arises from the fact that 
adding local degrees of freedom that form a ID represen- 
tation does not change the phase of state (see Fig. 9) It 
is clear that if the transformations t/* satisfy eqn. (12), 
from the second "equivalent" relation also satisfy 
eqn. (12). 

The solutions of eqn. (12) can be grouped into classes 
using the equivalence relations eqn. (13) and eqn. (14). 
Those classes should correspond different SPT states. 

We note that the condition eqn. (12) involves the whole 
many-body wave function. In appendix A, we will show 
that the condition eqn. (12) can be rewritten as a local 
condition where only a local region of the many-body 
wave function is used. Although we only discuss the 2D 
case in the above, similar result can be obtained in any 
dimensions. 

The discussions in the last a few sections outline some 
ideas that may lead to a classification of SPT phases. In 
this paper, we will not attempt to directly find all the 
solutions of eqn. (12) and to directly classify all the SPT 
phases. Instead, we will try to explicitly construct, as 
general as possible, the solutions of eqn. (12). Our goal 
is to find a general construction that produces all the 
possible solutions. 



VI. CONSTRUCTING SPT PHASES THROUGH 
GROUP COCYCLES 

In this section, we will construct solutions of eqn. (12) 
through the cocycles of the symmetry group G. The dif- 
ferent solutions will correspond to different SPT phases. 



A. Group cocycles 

The cocycles, cohomology group, and their graphic 
representations on simplex with branching structure are 
discussed in appendix D and E. Here we just briefly in- 
troduce those concepts. A d-cochain of group G is a 
complex function Vdigoj ffij gd) of 1 -f d variables in G 
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that satisfy 

Wd{go,9i,-;9d.)\ = l, (16) 

^d^\9o,gi, -^gd) = vd{ggo,ggi, -^ggd), g <^G 

where s{g) = 1 if g contains no anti-unitary time rever- 
sal transformation T and s{g) = —1 if g contains one 
anti-unitary time reversal transformation T. [When G is 
continuous, we do not require the cochain VdiSQigii ■■■i9d) 
to be a continuous function of gi. Rather, we only require 
Vd{goi giT--i gd) to be a so called measurable function of 
f/i.^'' A measurable function is not continuous only on a 
measure zero space.] 

The d-cocycles are special d-cochains that satisfy 



Il'^d"^^ (5o,--,ff»-i,5i-i-i,---,.9<i+i) = 1- 

i=0 

For d = 1, the 1-cocycles satisfy 

vi{gi,g2)vi{go,gi)/vi{ga,g2) = 1 
The 2-cocycles satisfy 

t^2(gi,g2,g3)'^2(go,gi,g3) ^ ^ 
'^2(go,g2,g3)f^2(go,gi,g2) 

and the 3-cocycles satisfy 

i^sigi, g2, 93, 94)1^3(90, gi> gsi gd'^sjgo, gi, g2, gs) 
(go, g2, g3, g4)i^3(go, gi, g2, gd 



(17) 



(18) 



(19) 



= 1 

(20) 



The d-coboundaries A,^ are special d-cocycles that can 
be constructed from the (d — l)-cochains ^id-i- 



^d{go,--,gd) = n^d-i^ (go,--,gj-i,gi+i,---,gd)- 



i=0 



(21) 



For d = 1, the 1-coboundaries are given by 

Ai(go,gi) = Mo(gi)/Mo(go) (22) 

The 2-coboundaries are given by 

h{go,gi,g2) = Mi(gi,g2)Mi(go,gi)//^i(go,g2), (23) 

and the 3-coboundaries by 

, / X M2(gi,g2,g3)M2(go,gi,g3) 

A3(go,gi,g2,g3) = — 7 ^ — -, r- (24) 

M2(go,g2,g3)Ai2(go,gi,g2) 

Two d-cocycles, Vd and v'^, are said to be equivalent iff 
they differ by a coboundary A^: Vd = v'd^d- The equiv- 
alence classes of cocycles give rise to the d-cohomology 
group H''[G,C/t(1)]. 

A d-cochain can be represented by a d-dimensional 
simplex with a branching structure (see Fig. 10). A 
branching structure (see appendix E) is represented by 





(b) 



FIG. 10: (Color online) (a) The line from go to g\ is a graphic 
representation of vi{go,gi). The triangle (go, 31,52) with a 
branching structure (see appendix E) is a graphic represen- 
tation of i^2(go, gi, 52). Note that, for the first variable, the 
go-vertex is connected to two outgoing edges, and for the last 
variable, the (;2-vertex is connected to two incoming edges, (a) 
can also be viewed as the graphic representation of eqn. (18) 
and eqn. (D25). The triangle corresponds to (di;^i)((?o, gi, 52) 
in eqn. (D25) and the three edges correspond to i^i (51,52), 
!/i(5o,gi) and v^^igo^gi)- (b) The tetrahedron (50,51,52,53) 
with a branching structure is a graphic representation of 
^^3 (50, 51, 52, 53). (b) can also be viewed as the graphic repre- 
sentation of eqn. (19) and eqn. (D26). The tetrahedron cor- 
responds to (d2t'2)(5o, 51, 52, 53) in eqn. (D26), and the four 
faces correspond to 1/2 (51 , 52 , 53 ) , V2 {go , 5i . 53 ) , i^2~ ^ (5o , 52 , 53 ) , 
and J^2~^(5o, 51.52)- 



arrows on the edges of the simplex that never form an 
oriented loop on any triangles. We note that the first 
variable go in Vdigo, gi, gd) corresponds to the vertex 
with no incoming edge, the second variable gi to the ver- 
tex with one incoming edge, and the third variable g2 
to the vertex with two incoming edges, etc . The condi- 
tions eqn. (18) and eqn. (19) can also be represented as 
in Fig. 10. For example, Fig. 10a has three edges which 
correspond to h'i{gi,g2), i^i(go,gi) and J^r^(go,g2)- The 
evaluation of a 1-cochain vi on the complex Fig. 10a is 
given by the product of the factors z/i((7i, (72), '^i(goigi) 
and i^i^igo, g2)- Such an evaluation will be 1 if vi is a 
cocycle. In general, the evaluations of cocycles on any 
complex without boundary are 1. 

Such a geometric picture will help us to obtain most 
of the results in this paper. 



B. (1 + 1)D case 

Let us discuss the ID case first. We will choose the 
ID SPT wave function to have a fixed form of a "dimer 
crystal" (see Fig. 7): 

|*psre) = ... <E) (^ 1^2 = g, A = g))® 

(El/52 = g,7i=g))®... (25) 

geG 

where we have assumed that physical states on each dot 
in Fig. 7 are labeled by the elements of the symmetry 
group G: ai,Pi e G. The dimmer in Fig. 7 corresponds 
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a 



g-lg* 



FIG. 11: (Color online) The graphic representation of 
eqn. (27). /2(ai, 02, 5, ff*) is represented by the polygon 
with a branching structure as represented by the arrows on 
the edge which never form a oriented loop on any triangle. 
1/2(01, g~^g*,g*) and 1/2(02, g~^g*,g*) are represented by the 
two triangles as in Fig. 10a. The value of the cocycle 1/2 on 
a triangle (say 1/2(01, g~^ g* , g*)) can be viewed as flux going 
through the corresponding triangle. 



to a maximally entangled state J2geG l'-'^2 = 9, Pi — g) ■ 

Next, we need to choose an on-site symmetry trans- 
formation (12) such that the state I^'psre) is invariant 
(where the two dots in each shaded box represent a site) . 
We note that U'^{g) acts on the states on the i site which 
are linear combinations of \ai,a2) in Fig. 7. Note that 
ai,a2 G G. So we can choose the action of ^(g) to be 
(see Fig. 11) 



W(g)\ai,a2) = 12(0^1,^2, g, g*)\gai, ga2) 



(26) 



where f 2(0-1, a2,g,g*) is a phase factor 
\f2(cti,a2,g,g*)\ — 1. We will use a 2-cocycle 
V2 & 'H^[G,Ut(1)] for the symmetry group G to con- 
struct the phase factor f2- (A discussion of the group 
cocycles is given in the appendix D.) 

Using a 2-cocycle 1^2, we construct the phase factor /2 
as the follows (see Fig. 11): 



f2(ai,a2,g,g* 



^2(«i,ff ^g*,g*) 
1^2(02, g^^g*,g*)' 



(27) 



Here g* is a fixed element in G. For example we may 
choose g* = 1. In appendix F, we will show that U'^(g) 
defined above is indeed a linear representation of G that 
satisfies eqn. (12). In this way, we obtain a SPT phase 
described by |*pSre) that transforms as U'^(g)- 

Note that here we only discussed a fixed SRE wave 
function. If we choose different cocycles in eqn. (27), the 
same wave function (25) can indeed represent different 
phases. One may wonder how a fixed SRE wave function 
can represent different quantum phases. 

To see this, let us examine how the state varies un- 
der the symmetry group. Notice that the phase factor 
f 2(011, a2,g,g*) is factorized, the basis \ai) varies as 

M(g)\ai) = i/2(ai,g^'^g* ,g*)\gai). 

The states \ai) form a representation of G itself, and 
the operator g transforms a state into another. The 
representation matrix element is given as M(g)a^^ga-^ — 
i^2(9~^g* , g* ,cti), and eqn. (27) can be rewritten as 



eqn. (26) we have U\g) = M(g)(»[M(g)]^ . Actually, this 
matrix M(g) is a projective representation of the group 
G, corresponding to the 2-cocycle 1^2- 

DiiTerent classes of cocycles V2 correspond to differ- 
ent projective representations. In the trivial case, where 
^2(01, g^^ g* ,g*) — 1, M(g) can be reduced into linear 
representations, and the corresponding SPT phase is a 
trivial phase. 

We will also show, in appendix F, that on a finite seg- 
ment of chain, the state l^'pSRs) has low energy excita- 
tions on the chain end. The excitations on one end of 
the chain form a projective representation described by 
the same cocycle 1/2 that is used to construct the solution 
U'^(g)- The end states and their projective representation 
describe the universal properties of bulk SPT phase. 

The different solutions of eqn. (12) constructed from 
different 2-cocycles do not always represent different SPT 
phases. If 1^2(501 5ii 52) satisfies eqn. (16) and eqn. (19), 
then 



^'2(90, 91,52) = V2(9o,gi,g2 



^ ^J■ligl,92)^J■l(go,gl) 
Ml (90,52) 



(28) 



f2(ai,a2,g,g* 



M(g) [M(g) c 



From 



also satisfies eqn. (16) and eqn. (19), for any /ii(5o,5i) 

satisfying Mi(9ffo,g5i) = Mi^^^(5o, ffi), g & G. So 
^^2 (.90, 52, 53) also gives rise to a solution of eqn. (12). 
But the two solutions constructed from 1^2(50,51,52) 
and 1^2(50,51,52) are related by a symmetric LU trans- 
formations (for details, see discussion near the end of 
appendix I). They are also smoothly connected since 
we can smoothly deform Mi (50,51) to /ii (50,51) = 
1. So we say that the two solutions obtained from 
^^2(50, 52, 53) and 1^2(50,52,53) are equivalent. We 
note that 2^2(50,52,53) and 1^2(50,52,53) differ by a 2- 
coboundary t^^^9^'9^^t^'-^3o,9i) ^ ^^^^ equivalence 

classes of 1^2(50,52,53) is nothing but the cohoniology 
group W^IG, Ut(1)]. Therefore, the different SPT phases 
are classified by 'H'^[G,Ut(1)]. 

We see that, in our approach here, the different SPT 
phases are not encoded in the different wave functions, 
but encoded in the different methods of fractionalizing 
the symmetry transformations U'^(g). 



C. (2-|-l)D case 

The above discussion and result can be generalized to 
higher dimensions. Here we will discuss 2D SPT state 
as an example. We choose the 2D SPT state to be a 
"plaquette state" (see Fig. 8) 

I^pSRe) = «'squaros(X! I"l ^9,^2^ 5, 73 = 5, ^4 = 5)) 

see 

(29) 

where we have assumed that physical states on each 
dot in Fig. 8 are labeled by the elements of the sym- 
metry group G: ai,Pi,... € G. The four dots in a 



13 




g-lg* 

FIG. 12: (Color online) The graphic represen- 

tation of the phase factor /3(ai, Q2, Q3, Q4, 5, g*) in 
eqn. (30). The arrows on the edges that never form 
a oriented loop on any triangle represent the branch- 
ing structure on the complex. The four tetrahe- 
drons give rise to V3{ai, 02, g* , g"), ^3(02,03, g~^g* , g*), 
iy^^{a4,a3,g~^g',g*), and v^^{ai,a4,g~^g*,g*) 

linked square in Fig. 8 form a maximally entangled state 
EgsG l"i = 9,^2 = 9,13 = 5, = g). We require that 
the state I^&pSre) is invariant under an on-site symmetry 
transformation (12) (where the four dots in each shaded 
square represent a site). 

To construct an on-site symmetry transformation (12), 
in 2 dimensions, the action of C/* is chosen to be 

C/*(.g)|ai,a2,a3,a4) (30) 
= /a (ai , "2 , as , "4 , 5i ff * ) , ga2 , gas , ga^) . 

Here /3(q!i, 012, as, 04, 5, g*) is a phase factor that cor- 
responds to the value of a 3-cocycle 1^3 e ■H'^[G, C/t(1)] 
evaluated on the complex with a branching structure in 
Fig. 12: 

/3(ai,a2,a3,a4,5:.9*) 
^ t^3(ai, «2, g~^g*,g*)i^3{a2,a3, 9~^9*,9*) ^g^x 
1^3 ("4, "3, 9~^9*,9*)'^3iai,a4, g-^g*,g*) ' 

In appendix G, we will show that W^ig) defined above 
is indeed a linear representation of G that satisfies 
eqn. (12). We will also show that (see appendix 1 and 
Ref. 59) in a basis where the many-body ground state 
is a simple product state, although ®iU'^{g) is an on-site 
symmetry transformation when acting on the bulk state, 
it cannot be an on-site symmetry transformation when 
viewed as a symmetry transformation acting on the effec- 
tive low energy degrees of freedom on the boundary when 
the 3-cocycle 1^3 is non-trivial. 

VII. SPT PHASES AND TOPOLOGICAL 
NON-LINEAR cr-MODELS 

A. The fixed-point action that does not depend on 
the space-time metrics 

In the above, we have constructed SPT states and their 
symmetry transformations using the cocycles of the sym- 



metry group. We can easily find the Hamiltonians such 
that the constructed SPT states are the exact ground 
states. In the following, we are going to discuss a La- 
grangian formulation of the construction. We will sys- 
tematically construct models in d -I- 1 space-time dimen- 
sions that contain SPT orders characterized by elements 
in T-L^'^'^[G, ^7^(1)]. It turns out that the Lagrangian for- 
mulation is simpler than the Hamiltonian formulation. 

A SPT phase can be described by a non-linear u-model 
of a field n{x, r), whose imaginary-time path integral is 
given by 

Z = j Dn e-/d''xdr£[ri(x,r)l_ ^^2) 

We will call the term q- I <i''^'ir c[n(x,r)\ ^-^^ action- 
amplitude. The imaginary-time evolution operator from 
Ti to T2, U[n2{x), ni[x),T2,Ti], can also be expressed as 
a path integral 

U[n2{x),n,{x),T2,n]= J Dn / ^''x/- d. £[r.(x..)] 

(33) 

with the boundary condition n(x,Ti) = ni{x) and 
n{x, T2) = n2{x). 

If the model has a symmetry, the field n transforms 
as n — >■ 5 • n under the symmetry transformation g g G. 
The action-amplitude has the G symmetry 

g-/ d'^xdr £[n.(cc,r)] ^ / d''cc dr £[g-n(x,r)] _ ^3^^ 

To understand the low energy physics, we concentrate on 
the "orbit" generated by G from a fixed ng: {5 • no;g € 
G}. Such an "orbit" is a symmetric space G/H where 
H is the subgroup of G that keeps n-o invariant: H = 
{h;h ■ no = no,h G G}. We can always add degrees 
of freedom to expand the symmetric space G/H to the 
maximal symmetric space, which is the whole space of 
the group G. So to study SPT phase, we can always 
start with a non-linear cr-model whose field takes value 
in the symmetry group G, the maximal symmetric space. 
Such a non-linear a-model is described by a path integral 

Z ^ J Dg ^ '^"'"^^ ^^si^'^'^l geG. (35) 

We would like to consider non-linear cr-models that 
describe a SRE phase with finite energy gap and fi- 
nite correlations. So a low energy fixed point action- 
amplitude e"-^'^ aijdr £[g(cc,r)] jj^^g^ j^q^ depend on the 
space-time metrics. In other words, the fixed-point 
non-linear cj-model must be a topological quantum field 
theory.®^ We will call such non-linear (T-model a topolog- 
ical non-linear tr-model. A trivial topological non-linear 
CT-model is given by the following fixed-point Lagrangian 
Cfix[g{x,T)] = which describes the trivial SPT phase. 

A non-trivial topological non-linear cr-model has 
a non-zero Lagrangian Cfix[g{x,t)] ^ 0. How- 
ever, the corresponding fixed-point action-amplitude 



14 



^- J d''xdT c[g{x,T)] (Joes not depend on the space-time 
metrics. One possible form of the fixed-point Lagrangian 
^{ix[gix,T)] is given by a pure topological 6'-term. As 
stated in section XI, the origin of the topological 0-term 
may be the Berry phase in coherent state path integral. 
For a continuous non-linear cr-model whose field takes 
values in a continuous group G, the topological 0-term 
is described by the action-amplitude c-l^ ^'^'^ Ctopa[g{x,T)] 
that only depends on the mapping class from the space- 
time manifold M to the group manifold G. Such 
kind of topological term is given by a closed (1 + d)- 
form wi+rf on the group manifold G which is classified 
by H^~^'^{G,R). The corresponding action is given by 
/ d'^xdr £topo[9{x,T)] = / uji+d- 

The other possible form of the fixed-point Lagrangian 
Cfi4g{x,T)] is given by a WZW term.63.64 rpj^^ 
term is described by the action 5wzw that cannot be 
expressed as a local integral on the space-time manifold 
M. That is to say, we cannot express ^wzw as 5wzw = 
J d'^xdr £yi/zw[g{x,T)]. We have to view the space-time 
manifold M as a boundary of another manifold M^xt in 
one higher dimensions, M = 9Afextj and extend the field 
on M to a field on Mext- Then the WZW term S'wzw can 
be expressed as a local integral on the extended manifold 
Mext: 

S'wzw = / d^+'^xdr Cwzw[gix,T)] (36) 

such that S'wzw mod 27ri does not depend on how we 
extend M to Afoxt- A WZW term is given by a quantized 
closed {d + 2)-form LOd+2 on the group manifold G: 

S'wzw = / (37) 

which clearly does not depend on the space-time metrics. 
Later, we will show that WZW terms in (c?-|-l)-dimension 
space-time and for group G are classified by the elements 
in H'^+2[G,[/t(1)]. 

We see that both the topological 0-term and the WZW 
term do not depend on the space-time metrics. So the 
fixed-point Lagrangian may be given by a pure topolog- 
ical 0-term and/or a pure WZW term. 

B. Discrete and continuous non-linear it- model 

We would like to stress that the topological 0-term and 
the WZW term discussed above require both continuous 
group manifold and continuous space-time manifold. 

On the other hand, in this paper, we are considering 
quantum disordered states that do not break any symme- 
try. So the field g{x,T) fiuctuates strongly at all length 
scale. The low energy effective theory has no smooth 
limit. Therefore, the low energy effective theory must be 
one defined on discrete space-time. 

For discrete space-time, we no longer have non-trivial 
mapping class from space-time to the group G, and we no 




FIG. 13: (Color online) The graphic representation of the 
action-amplitude e"^*-^^'*' on a complex with a branching 
structure represented by the arrows on the edge. The vertices 
of the complex are labeled by i. Note that the arrows never 
for a loop on any triangle. 

longer have topological 0-term and WZW term. In this 
section we will show that although a generic topological 
O-term cannot be defined for discrete space-time, we can 
construct a new topological term on discrete space-time 
that corresponds to a quantized topological O-term in the 
limit of continuous space-time. Here a quantized topolog- 
ical 0-term is defined as a topological 0-term that always 
gives rise to an act ion- amplitude e-^^ xdr CR^[g{x,T)] _ ^ 
on closed space-time. We will also call the new topolog- 
ical term on discrete space-time a quantized topological 
6'-term. 

To understand the new topological term on discrete 
space-time, let us start with a continuous non-linear cr- 
model whose field takes values in a group G: g{x, t). The 
imaginary-time path integral of the model is given by 

Z = J Dg e^/d'^^drrlgCx.r)] (3g) 

with a symmetry described by G: 

^- f d^xdr C[g{x,T)] ^ ^-f d''xdr C[gg{x,r)] ^^q 

(39) 

If we discretize the space-time into a complex with a 
branching structure (such as the complex obtained by 
a triangularization of the space-time manifold) , the path 
integral can be rewritten as (see Fig. 13) 

Z-|Gr~"5]e-^«f'» 
^-S({g.})^ -Q ..^^--(^^g^,...,^,) (40) 

{i3...k} 

where e"'^*^^^'^^ is the action-amplitude on 
the discretized space-time that corresponds to 
^~ J d xdT c[g(x,T)] q£ -(-j^p continuous non-linear cr- 
model, and ''(ffi, <?j, ...,5?:) corresponds to the 
action-amplitude q~ hi,i. --.k) ^ ^'^^ ^^s(x,t)] ^ single 
simplex (i, j, k). Also, Sij,,,k = ±1 depending on the 
orientation of the simplex (which will be explained in 
detail later). 
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Here on each vertex of the space-time complex, we have 
a gi £ G. Qi corresponds to the field g{x,t) and J2{gi} 
corresponds to the path integral / Dg in the continuous 
non- linear o"- model. \G\ is the number of elements in G, 
Ny is the number of vertices in the complex. 

We note that, on discrete space-time, the model can 
be defined for both continuous group and discrete group. 
When G is a continuous group, the above model is a 
continuous non-linear ct- model (in this case, 
should be interpreted as an integral over the group man- 
ifold / dgi). When G is a discrete group, the above model 
is a discrete non-linear (T-modcl. 

We see that a non-linear cr-model on {d + 1)D dis- 
crete space-time is described by a complex function 
'^i+d{gi,g],—,gk)- Different choices of vi+d{gi,gj, ...,gk) 
give different theories/models. 

So we would like to ask: what vi+digiigj, ■■■i9k) 
will give rise to a quantized topological 0-term on dis- 
crete space-time? Very simply, we need to choose 
'^i+d{gi,gj,-,9k) so that 

n ^l+d''i9i,9j,-,gk) = l (41) 

{ij...k} 

on every closed space-time complex without boundary. 

There are uncountably many choices of 
i^i+digiigji ■■■,gk) that satisfy the above condition 
and give rise to quantized topological ^-terms. However, 
we can group them into equivalent classes, and each 
class corresponds to a type of quantized topological 
6'-terms. We will show later that the types of quantized 
topological ^-terms are classified by ^^+'^[G, Ut{1)]- So 
we can have non-trivial quantized topological O-tCTias 
discrete space-time only when 'H^+'^fG, J7t(1)] is non 
trivial. The number of equivalence classes of non-trivial 
quantized topological 0-terms is given by the number of 
the non-trivial elements in H^^'^[G, C/t(1)]. 

Prom the above discussion, it is also clear that we can- 
not generalize un-quantized topological 0-terms to dis- 
crete space-time. So on discretized space-time complex, 
the only possible topological 6-terms are the quantized 
ones. 

After generalizing quantized topological ^-terms to dis- 
crete space-time, we can now generalize WZW term to 
discrete space-time. We will call the generalized WZW 
term a non-local Lagrangian (NLL) term. To construct 
a NLL term on a closed (d + l)D space-time complex 
Md+i, we first view Md+i as a boundary of a (d + 2)D 
space-time complex Md+2- We then choose a function 
^i+digi^gj^ ■■■igk) that defines a quantized topological 9- 
term on the {d + 2)D space-time complex Md+2- Then 
the action-amplitude of i'2+d{9i,gj^ ■■■,gk) on Md+2 only 
depends on the gi on Md+i = dMd+2, the boundary of 
Md+2- Thus such an action-amplitude actually defines a 
theory on the {d + 1)D space-time complex Md+i- Such 
an action-amplitude is the NLL term on the space-time 
complex Md+i- We will see that the types of NLL terms 
on {d+ 1)D space-time complex and for group G are clas- 
sified hy H'^+'^[G,U t {I)] - 



We would like to stress that the proper topological non- 
linear cj-models are for disordered phases, and they must 
be defined on discrete space-time. Only quantized topo- 
logical 0-terms can be defined on discrete space-time. On 
the other hand, the WZW term can always be general- 
ized to discrete space-time, which is called NLL term. 
Both quantized topological 0-terms and NLL terms on 
discrete space-time can be defined for discrete groups. 



C. Quantized topological S-terms lead to gapped 
SPT phases 

We know that the action-amplitude defines a physi- 
cal model, in particular, defines imaginary-time evolu- 
tion operator U{ti,T2). For a SPT phase, its fixed point 
action-amplitude must have the following properties (on 
a closed spatial complex): 

(a) The singular values of the imaginary-time evolution 
operator ?7[gi(Ti),5i(r2),Ti,T2] are I's or O's. 

(b) The singular values of the imaginary-time evolution 
operator C/[(7i(Ti), gj(r2), ti, T2] contain only one 1. 

Usually, the imaginary-time evolution operator is given 
by U{ti,T2) = e~('^^~'^i'-^. One expects that the log of 
the eigenvalues of U{ti,T2) correspond to the negative 
energies. However, in general, the basis of the Hilbert 
space at different time r can be chosen to be different. 
Such a time dependent choice of the basis corresponds 
to adding a total time derivative term to the Lagrangian 
£ — >■ £ -I- It is well known that adding a total time 
derivative term to the Lagrangian does not change any 
physical properties. For such more general cases, the log 
of the eigenvalues of U{ti,T2) do not correspond to the 
negative energies, since the eigenvalues of C/(ti,T2) may 
be complex numbers. In those cases, the log of the sin- 
gular values of U{ti, T2) correspond to the negative ener- 
gies. This is why we use the singular values of U{ti,T2) 
instead of the eigenvalues of U{ti,T2)- 

At the low energy fixed point of a gapped system, the 
fixed-point energies are cither or infinite. Thus the 
singular values of the imaginary-time evolution operator 
are either 1 or 0. For a SPT phase without any intrinsic 
topological order and without any symmetry breaking, 
the ground state degeneracy on a closed spatial complex 
is always one. Thus the singular values of the imaginary- 
time evolution operator contain only one 1. 

For the action-amplitude given by a quantized topolog- 
ical ^-term, its corresponding imaginary-time evolution 
operator does have a property that its singular values 
contain only one 1 and the rest are O's. This is due to 
the fact that the action-amplitude for each closed path 
is always equal to 1. So a quantized topological ^-term 
indeed describes a SPT state. 
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D. NLL terms lead to gapless excitations or 
symmetry-breaking boundary states 



state. Therefore the boundary state must be gapless, or 
break the symmetry. Thus, 



On the other hand, if the fixed-point action-amphtude 
in (d + 1) space-time dimension is given by a pure NLL 
term, its corresponding imaginary-time evolution opera- 
tor, we believe, does not have the property that its sin- 
gular values contain only one 1 and the rest are O's, since 
the action-amplitude for different closed paths can be dif- 
ferent. 

In addition, if the pure NLL term corresponds to a 
non-trivial cocycle Vdj^2 in ■^''"'"^[G', C/t(1)], adding dif- 
ferent coboundary to 1^^-1-2 will lead to different action- 
amplitude on closed paths. There is no coboundary that 
we can add to the cocycle Vd+2 to make the action- 
amplitude for closed paths all equal to 1. Further more, 
a renormalization group flow only adds local Lagrangian 
term 5C that is well defined on the space-time complex. 
The renormalization group flow cannot change the NLL 
term and cannot change the corresponding cocycle i'd+2, 
which is defined in one higher dimensions. This leads us 
to conclude that an action-amplitude with a NLL term 
cannot describe a SPT state. Therefore 



An action- amplitude with a NLL term must have 
gapless excitations, or degenerate boundary ground 
states due to symmetry breaking. 



The above is a highly non-trivial conjecture. Let us 
examine its validity for some simple cases. Consider a G 
symmetric non-linear a-model in (1-1-0) dimension which 
is described by an action-amplitude with a NLL term. In 
(1+0) dimension, the NLL term is classified by 2-cocycles 
V2 in H^[G, J7t(1)], which correspond to the projective 
representations of the symmetry group G. So the ground 
states of the non-linear tr-model form a projective repre- 
sentation of G characterized by the same 2-cocycle 1^2- 
Since projective representations are always more than 
one dimension, (1+0)D systems with NLL terms cannot 
have a non-dcgencrate ground state. In (H-l)-dimension, 
continuous non-linear a-modcls with the WZW term are 
shown to be described by the current algebra of the con- 
tinuous symmetry group and are gapless.^'' In Ref. 59, we 
further show that discrete non-linear cr-models with the 
NLL term in (l-l-l)D must be gapless if the symmetry is 
not broken, for both continuous and discrete symmetry. 
The above conjecture generalize such a result to higher 
dimensions. 

We note that the boundary excitations of the SPT 
phases characterized by (l-l-d)-cocycle vi+d a-re described 
by an effective boundary non-linear cr-model that con- 
tains a NLL term characterized by the same (1 -f d)- 
cocycle vi+d- 

As discussed before, a non-linear cr-model with a non- 
trivial NLL term cannot describe a SPT state. Thus the 
boundary state must be gapless, or break the symmetry, 
or have degeneracy due to non-trivial topological order. 
However, the SPT state is a direct product state. The 
degrees of freedom on the boundary also form a product 



A non-trivial SPT state described by a non-trivial 
{l-\-d) -cocycle must have gapless boundary excita- 
tions or break the symmetry at the boundary. 



We would like to stress that the symmetry plays a very 
important role in the above discussion. It is the reason 
why the non-linear a-model field gi takes many different 
values. If there was no symmetry, at low energies, the 
non-linear cr-model field gi would only take a single value 
that minimizes the local potential energy. In this case, 
there were no non-trivial topological terms. 



VIII. CONSTRUCTING SYMMETRIC 
FIXED-POINT PATH INTEGRAL THROUGH 
THE COCYCLES OF THE SYMMETRY GROUP 

In the last section, we argue that SPT phases in d- 
dimension with on-site symmetry G are described by 
quantized topological 0-terms. In this section, we are go- 
ing to explicitly construct quantized topological 0-terms 
that realize those SPT orders in each space-time dimen- 
sion. We will also show that the quantized topological 
6'-tcrms are classified by n^+'^[G,UT{l)]. 



A. (1+1)D symmetric fixed-point action-amplitude 

Let us first discuss (1+1)D fixed-point action- 
amplitude with a symmetry group G. For a (1 -I- 1)D 
system on a complex with a branching structure, a fixed- 
point action-amplitude {ie a quantized topological 9- 
term) has a form (see Fig. 13) 



{ijk} 

= '^2^(91,92, 33)^^2(50, 94, 93)1^2^ {95, 90, 51 J X 

1^2 (51, go, 53) 1^2' ^(55, 90,54) (42) 

where each triangle contributes to a phase factor 
i^2''°i9i,9j,9k), I\{ijk} multiply over ah the triangles in 
the complex Fig. 13. Note that the first variable gi in 
V2{gi,gj,gk) corresponds to the vertex with two out go- 
ing edges, the second variable gj to the vertex with one 
out going edge, and the third variable gk to the vertex 
with no out going edge. Sijk = ±1 depending on the ori- 
entation of i — J — /c to be anti-clock-wise or clock-wise. 

In order for the action-amplitude to represent a quan- 
tized topological 0-term, we must choose 1^2(91, gj, gk) 
such that 



-S({g^}) _ 



n 

{tjk} 



^2"'{9t,93,9k) = 1 



(43) 
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FIG. 14: (Color online) Graphic representation of 

V2{go-,gi,g2)v2{go, 92,93) = i^2{gi,g2, 93)1^2(90, gi, 93) The ar- 
rows on the edges represent the branching structure. 





FIG. 15: (Color online) Graphic representation of 

V2{gi,g2,g3) = (50,51, 52) 1^2 (50, 52, 53) 2^2" Hso, 91,53). The 
arrows on the edges represent the branching structure. 



on closed space-time complex without boundary, in par- 
ticular, on a tetrahedron with four triangles (see Fig. 10): 



-S{{g^}) _ 



n 



^2 '"{91,9 3, 9k) 



{ijk} 

'^2(gl,g2,g3)'^2(ff0,gl,g3) 

'^2{9o, 91, 92)1^2(90, 92, 93) 



(44) 



Also, in order for our system to have the symmetry gen- 
erated by the group G, its action-amplitude must satisfy 



-S{{g^}) _ p-S({gg,}) 



if g contains no T 



(e-^«S'») 



e -S({9.})) ^ Q-s{{gg,})^ g contains one T (45) 



where T is the time-reversal transformation. This re- 
quires 



^^2 {9t , 9 J ,9k) ^ V2 {99t , 99] ,99k)- 



(46) 



Eqn. (45) and eqn. (46) happen to be the conditions of 
2-cocycles 1^2(50,31,32) of G. Thus the action-amplitude 
eqn. (42) constructed from a 2-cocycle 1^2(90, 9i, 92) is a 
quantized topological 0-term. 

If ^2(90, 92, 93) satisfies eqn. (45) and eqn. (46), then 

V2(.9o,92,93) = V2{9o, 92,93) ^ (47 

Ml (50, 52) 

also satisfies eqn. (45) and eqn. (46), for any /iti (50,51) 
satisfying fii{ggo,99i) = f^i{9o,9i), 9 e G. So 
^2 (50, 52, 53) also gives rise to a quantized topologi- 
cal 6'-term. As we continuously deform /ii (50,51), the 
two quantized topological ^-terms can be smoothly con- 
nected. So we say that the two quantized topological 
6'-terms obtained from 1^2(90,92,93) and 1^2(50,52,53) are 




(b) 



FIG. 16: (Color online) (a) The graphic representation of 
*({5i}cdgo) = i^2'^(5i, 52, 53)i^2~^(5i, 53, 55)2^:^^53, 54, 55) 
(b) The graphic representation of ^'Uls^lcdge) = 
1^2(51, 52, 53)2^2(51, 53, 55)i'2(53, 54, 55) The arrows on the 
edges represent the branching structure. 



equivalent. We note that 1^2(90,92,93) and 1^2(90,92,93) 
differ by a 2-coboundary ''^^^'f,fgl'\[]°'^''' ■ So the set of 
equivalence classes of 1^2(90,92,93) is nothing but the co- 
homology group H^[G, C/t(1)]. Therefore, the quantized 
topological 0-terms are classified by H^[G, C/t(1)]. 

We can also show that eqn. (42) is a fixed- 
point action-amplitude from the cocycle condi- 
tions on V2(9i,9j,9k)- From the geometrical pic- 
ture of the cocycles (see Fig. 10), we have the 
following relations: 1^2(90,91,92)1^2(90,92,93) = 
^2(91,92,93)1^2(90,91,93) (see Fig. 14) and 

^2(91,92,93) = 1^2(90, 91, 92)1^2(90, 92, 93)1^2^ (90, 91,93) ■ 
(see Fig. 15). We can use those two basic moves to 
generate a renormalization flow that induces a coarse- 
grain transformation of the complex. The two relations 
Fig. 14 and Fig. 15 imply that the action-amplitude 
is invariant under the renormalization flow. So it is 
a fixed-point action-amplitude. Certainly, the above 
construction applies to any dimensions. 



B. (1+1)D fixed-point ground state wave function 

For our fixed-point theory described by a quan- 
tized topological 0-term, its ground state wave function 
^({5i}) can be obtained by putting gi on the edge of a 
disk and making a triangularization of the disk (see Fig. 
13). We sum the action-amplitude over the gi on the in- 
ternal vertices while fixing the 5i's on the edge (see Fig. 
16a): 



*({ft}odgc) 



pi G internal 



|G|^" 



n ^2"'(9^,93,9k) (48) 
{ijk} 



where X]g eintcrnai sums over gi on the internal vertices 
and i\pnt'=''nai jg ^j^g number of internal vertices on the 
disk. 

Clearly the ideal wave function 5'({5i}cdgc) satisfies 



«''^^H{5*}cdgo) - ^(te ledge), |*({5»}edgo)| = 1, 



(49) 
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FIG. 17: Two solid tetrahedrons poffiff2f/4, <7o5253<74 
and three soHd tetrahedrons §0919293, 90919394-, 91929394 
occupy the same volume, which leads to the graphic 
representation of 2^3(50, 5i, 2, 54)2^3(50, S2, S3, 54) = 
^3 {90 , 51 , 52 , 53 )z^3 (30 ,91,93,94)1^3 (91 ,92,93,94) (see 
eqn. (20)). The arrows on the edges represent the branching 
structure. 




S2 (a) 82 (b) 



FIG. 18: One solid tetrahedron 91929394, and 

four solid tetrahedrons 90919294, 50525354, 5o5i5354, 
51515253 occupy the same volume, which leads to 
the graphic representation of 1^3(91,92,93,94) = 
1^3(90, 91, 92, 94)1^3(90, 92, 93, 94)1^3^ {90, 91, 93, 94)1^3 ^{90, 91, 92, 93) 
(see eqn. (20)). The arrows on the edges represent the 
branching structure. 



which represents a symmetric state. We also note that 
^^({ffilcdgc) can be represented by Fig. 16b, since the 
product of the wave functions in Fig. 16a and Fig. 16b 
is the value of the cocycle on a sphere which is equal to 
1. 



C. (2+l)D symmetric fixed-point action-amplitude 

In (2 + 1)D, our ideal model with on-site symmetry G 
is defined by the act ion- amplitude on a 3D complex with 
gi € G on each vertex: 

e-S({.4)= -Q ^3^--(g„g^,5,,50 (50) 
{i]ki} 

where 1^3(91, gj, gk, gi) is a three cocycle and Ilfiifei} mul- 
tiply over all the tetrahedrons in the complex Fig. 17. 
The 3D complex has a branching structure. The first 
variable gi is on the vertex with no incoming edge, the 
second variable gj is on the vertex with one incoming 
edge, etc . Also Sijki — ±1 depending on the orientation 
of the ij/cZ-tetrahedron. On a close space-time complex, 
the above action-amplitude is always equal to 1 due to 
the cocycle condition on v^{gi, gj ,gk,gi)- Thus the above 
action-amplitude is a quantized topological 0-term. 

The conditions of 3-cocycle lead to the two relations 
in Fig. 17 and Fig. 18. These lead to a renormalization 
flow of the complex in which the above action-amplitude 
is a fixed-point action-amplitude. The fixed-point action- 
amplitude leads to an ideal short-range-entangled state 
(see section IX) that has a symmetry G and is character- 
ized by 1/3 G n^[G,UTil)]. 



D. (d -I- 1)D symmetric fixed-point action-amplitude 

Through the above two examples in (l-l-l)D and 
(2-|-l)D, we see that the {d + 1)D symmetric fixed-point 



action-amplitude is given by 

n (51) 

' ' {^J■■■k} 

where gi is associated with each vertex on the space- 
time complex and is the number of vertices. 
J2{g } sums over all possible configurations of {gi} and 

'^i+d{gi,gj,---,gk) is a {l + d)-cocydem'H^+'^[G,UT{l)]. 

When the space-time complex is closed 
(ie has no boundary), the action-amplitude 
Il{^J...k}K+d''i9^,9j,■■■,9k) is always equal to 1. 
Thus the act ion- amplitude represents a topological 
O-term. 

When the space-time complex has a boundary, the 
action-amplitude will not always be equal to 1 and is 
not trivial. We note that, due to the cocycle condition 
on I'l'^j^^ {gi, gj, ...,gk), such a action-amplitude will only 
depend on g^'s on the boundary of the space-time com- 
plex. Thus such an action-amplitude can be viewed as 
an action-amplitude of the boundary theory. 

As an action-amplitude of the boundary theory, 
Y\{i3...k}K+d''{9t,9j.--,9k) is actually a NLL term, 
which is a generalization of the WZW topological term 
for continuous non-linear tr-models to discrete non-linear 
CT-models. So the boundary excitations of our model 
defined by eqn. (51) are described by a non-linear a- 
modcl with a NLL term composed by the same i^i+d S 
'H^^'^[G,Ut{X)\- We see a close relation between the 
topological 0-term in (d + 1) space-time dimensions and 
the NLL term in d space-time dimensions. An example 
of such a relation has been discussed by Ng for a (l-l-l)D 
model with 50(3) symmetry.^^ 

When vij^d is non-trivial, we believe that the bound- 
ary excitations are gapless if the symmetry is not broken 
on the boundary. (To make sure that the symmetry is 
not broken on the boundary, we may need to add non- 
topological term 5C to the boundary Lagrangian.) 

In the following, we will show that the ground state 
wave function of our model (51) describes a SPT state. 



19 




FIG. 19: (Color online) The graphic representation of 
eqn. (52). The boundary is the complex M, and the whole 
complex Mext is an extension of M. 

IX. TRIVIAL INTRINSIC TOPOLOGICAL 
ORDER IN OUR FIXED-POINT MODELS 

The ground state of our d-dimensional model (51) is a 
wave function m on M, a c?-dimensional complex. It is 
given by (see Fig. 19) 

^M({g»}Af) = ^^'^!n!rn"i' n ^T+d''i9t.93^-:9k) 

(52) 

which generalizes eqn. (48) from (1+1)-D to (d+l)-D. We 
use Afoxt to denote a (d + 1) dimensional complex whose 
boundary is M. {gi\M are on the vertices on M and 
Gintornai sums ovcr the giS on the vertices inside the 
complex Mcxt (not on its boundary M). Also Ilfij k} 
product over all simplices on A/^xt- 

Due to the cocycle condition satisfied by 
ui+d{9i,9],—,gk), we see that, for fixed {gj^, the 
product ri{ij...fc} f'l+d "(ffiiffj: ■■■,9k) does not depend on 
gi^s on the vertices inside the complex Mcxt- Thus 

*M({g.}A/)= n '^Ud'i9^:93,■■■,9*)■ (53) 
{ij...*} 

if we choose Afoxt to be M plus one more vertex with 
label g* (see Fig. 19). The state on M (the boundary of 
Fig. 19) does not depend on the choice of g* . 

Using the above expression, we can show that the 
ground state wave function of our fixed-point model is 
SRE state with no intrinsic topological orders. Let us 
first write the ground state of our fixed-point model in a 
form 

l*A^)= E n ''ltd'(9^,9J,■■■,9*)\{9^}M), (54) 

{ai}M {ij---*} 

where form a basis of our model on d-dimensional 

complex M . The on-site symmetry acts in a simple way: 

9 ■ \{9i}m) |{5.9Jm), g & G (55) 

We note that if we choose the particular form of Aicxt 
in Fig. 19 to obtain state $m on M, the phase fac- 
toi' '^i+d' i9i, 9], ■■■, 9*) can be viewed as a LU 



transformation. We can write \^m) in a new basis 

|{5«}m)' = '^l+d ''^9i,93,■■■,9*)\{9^}M)■■ 

I*m) = J2 \i9^}My■ (56) 

{gi}M 

Thus, on any complex Af that can be viewed as a bound- 
ary of another complex Afoxt, the state on Af can be 
transformed by an LU transformation into a state that 
is the equal weight superposition of all possible states 
KSiIm) on Af. The wave function in the new bases is 
very simple, which is actually a product state. In ap- 
pendix H, we will show that under a dual transformation, 
this product state is equivalent to the canonical form of 
wave function discussed in Sec. IV and V. 

We have used the (l + d)-cocycles in Ti^+'^IG, C/t(1)] to 
construct our fixed-point models which have ground state 
wave functions that also depend on the (1 -I- (i)-cocycles. 
In the above, we have shown that all those states can be 
mapped to the same simple product state via LU trans- 
formations. Does this mean that those states from dif- 
ferent (1 -I- d)-cocycles all belong to the same phase? The 
answer depends on if symmetry is included or not. 

If we do not include any symmetry, those states from 
different (1 -I- (i)-cocycles indeed all belong to the same 
trivial phase. Thus our fixed-point states constructed 
from different (1 + (i)-cocycles all have trivial intrinsic 
topological order. This is consistent with the fact that 
the fixed-point partition function on any space-time com- 
plex has the form 

Z = e-^o^, (57) 

where V is the volume of the space-time complex (say 
V is the number of simplices in the space-time complex) . 
We would like to stress that the above expression is exact. 
So after we remove the term that is proportional to the 
space-time volume, we have Z — 1. This means that 
the ground state is not degenerate on any closed space 
complex, which in turn implies that the ground state 
contains no intrinsic topological order. 

On the other hand, if we include the on-site symme- 
try G, states from different (1 + c?)-cocycles belong to 
the different phases which correspond to different SPT 
phases. This is because the LU transformation rep- 
resented by I\{tj...*}'^itd'i9i,9j,^^;9*) is not a sym- 
metric LU transformation under the on-site symmetry 
G. To see this, we first note that the LU transforma- 
tion I\{,j...*}'^itd''i9i,9j,^^;9*) contains several layers 
of non-overlapping terms. For example, for the (l-l-l)D 
system in Fig. 19, the LU transformation has two layers 

Yl '^2{gi,gj,gk) = ^2(93, 92, g*)'^2{g5, g4, g*)]^- 

{ijk} 

^2(92, 91, 9*)i^2{g4, g3, g*)>^2{gi, g5, 9*)] (58) 

In order for the LU transformation to be a symmetric, 
each local term, such as 1/2(92, 9i, 9*), must transform as 

V2^\92,9i,9*) = V2{9g2,ggi,g*) (59) 
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under the on-site symmetry transformation generated by 

g e G: Although vl^^\g2,gi,g*) = ^2(992, 991,99*), 

in general 1^2^^^92,91,9*) 7^ M992,99i,9*)- In 
fact, only trivial cocycle in 'H^+'^[G, C/t(1)] can satisfy 

Vi+d{9i,92,--,gi+d,9*) = vi+d{g9i,992,--,99i+d,9*)- 
Thus the fixed-point states from different (1 -|-d)-cocycles 
belong to the different SPT phases. 

We have seen that we can use two different basis 
K^iIm) and \{gi\M)' to expand the fixed-point wave 
function \^m)- The old basis \{gi}M) transforms sim- 
ply under the symmetry transformation: |{3i}A/) — >■ 
\{ggi\M)- But the wave function '^M{{9i}) in the old 
basis is complicated. In the new basis, the wave function 
is very simple ^''^f ({.^i}) = 1- But the symmetry trans- 
formation is more complicated in the new basis which 
will be discussed in appendix H. 

In section VI, we discuss the SPT phase by starting 
with a simple many-body wave function, and try to clas- 
sify all the allowed on-site symmetry transformations. 
Such a formalism is closely related to the new basis. 



X. EQUIVALENT COCYCLES GIVE RISE TO 
THE SAME SPT PHASE 

The ground state wave function M{{9i} m) of a SPT 
phase is constructed from a cocycle ui+d as in eqn. (53). 
Let 

1^'1+d ^ cocycle that is equivalent to vi+d- That is 
Vl-^-d and only differ by a coboundary 

i^i+di9o, ■■■,9i+d) 

l+d 

^ iyi+d(go, ■■; gi+d)]^ ^J.d^^ ...) (60) 

1=0 

where fJ-d{9o, gd) is a d-cochain. Then will give 
rise to a new ground state wave function ^j\,/({.9i}Af) 
of a SPT phase. One can show that 4'M({3i}Af) and 
*'M({ffi}M) are related: 

*m({5Jm) = *A/({ffjM) n ^'T' i9^,gJ,■■■) (61) 

{ij-} 

where Y[{ij } multiply over all the d-simplices in M. 
Note that, when we calculate ^'j\/({5i}A/), the terms 
^digi, 9j, g*) containing g* all cancel out. Due to 
the cochain condition eqn. (16) satisfied by /i^, the fac- 
tor Y[{ij }l^d^ {gi,gj,-- ) actually represents a sym- 
metric LU transformation. So the two wave func- 
tions ^ M{{gi}M) and '^'M{{9i}M) belong to the same 
SPT phase. Hence equivalent cocycles give rise to the 
same SPT phase, and different SPT phases are classi- 
fied by the equivalence classes of cocycles which form 

n^+''[GMi)]. 



XI. RELATION BETWEEN COCYCLES AND 
BERRY PHASE 

In this section, from path integral formalism, we will 
discuss some relations between the Berry phase and the 
cocycles that we used to construct topological non-linear 
cr-models. The Berry phase is defined in continuum pa- 
rameter space, so we need to embed the discrete sym- 
metry group G into a continuous group G, with G C G. 
For example, a discrete rotation group can be embedded 
into the 5*0(3) group. The coherent state path integral is 
performed in forms of G. After obtaining the topological 
(?-tcrm, we will reduce the symmetry group back to G. 

Suppose a rotation operator g is a symmetry opera- 
tion g £ G, and g\g)Q oc |(7)o is its eigenstate. The spin 
coherent state is defined as the following 

\g{m)) = R{m)\g)Q, 

where R{m) € G. We can write \g{m)) as \m) for sim- 
plicity. They satisfy the complete relation 

j Am\m) (m| cx I, 

where the integration is performed over the group space 
of G. 

The Berry phase in the spin coherent path integral is 
very important in our discussion. For a non-symmetry 
breaking system, the low energy effective theory can be 
written as the following path integral 

Z = j Dm{x,T)exp{— J d'^xdTCo{m) + iStop}, 

Stop = J d^'xdrCtopiAAo), (62) 

where £0 is the dynamic part of the Lagrangian which 
respects the symmetry group G (it is not important at 
the fixed point), and 5top is the topological 0-term of the 
action, which respects the enlarged symmetry group G. 
The 'gauge' field is defined a.s A = (m(a;, r)|V|m(a;, r)) 
and Aq = {m{x,T)\dr\m{x,T)) . The following is a gen- 
eralization of result of 0(3) nonlinear sigma model dis- 
cussed in Ref. 89. 

At zero temperature, the partition function only con- 
tains the contribution from the ground state. Under pe- 
riodic boundary condition, the ground state is a singlet, 
as a consequence, the Berry phase is trivial (integer times 
2tt). Under open boundary conditions, the Berry phase 
is contributed from the edge states. The topological 9- 
term is dependent on dimension. We will study it case 
by case. 

In (1 + I)D, the topological 0-tcrm is given as 
Stop ^0<f dxdrF, 

JSixSi 

where F = O^Aq — OqA^, and Si x 5*1 is the space-time 
manifold. 9 is an important constant which determines 
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Berry phase Cocycle 

FIG. 20: (Color online) Relations between Berry phase (of 
the end spin) in the loop {\go), \gogi), \gog1g2}) and 2-cocycle 
^^2(50, SoSi, 5o<?i52)- In the group space of 50(3), the two ends 
of a diameter stand for the same group element and can be 
seen as the same point. When 6 = 27r, the Berry phase is equal 
to zero when the the loop intersect the shell even times and 
is equal to vr when the loop intersection the shell odd times. 
If 9 is equal to even times of 27r, the corresponding 2-cocycle 
is trivial. If 9 is equal to odd times of 2n, the corresponding 
2-cocycle is nontrivial. 

the topological properties of the system. Under periodic 
boundary condition, the above integral is quantized and 
is equal to an integer (Chern number) times 27r which 
results in a trivial phase e^'^'^p = 1. However, at open 
boundary condition (where the space-time manifold be- 
comes a cylinder), the integral is not quantized. From 
Stokes theorem, it is determined by the boundaries, 

with (similar expression for S^) 

JSi 

= e i d\{T)A^[\{T)] 

J Si 

= e [ d^XF, (63) 

where S[ is a path in the parameter space (i.e., the group 
space of G, which is parameterized by A), Di is the 
area enclose by S[, and A\, F are the Berry connection 
and Berry curvature in the parameter space, respectively. 
The cyclic path S[ can be chosen as a sequence of symme- 
try operators in the symmetry group G. A closed path 
contains at least three points |5o)il5o5i)i and \gag1g2) 
(see figure. 20). The above integral gives a 2-cocycle or 
a product of 2-cocycles if we choose a proper gauge ( ie , 
multiply a proper coboundary) 

e''^^ = 1^2(90,9091,909192), 

where 90 is an arbitrary symmetry operator in G. 

In (2-|- 1)D, the possible topological 9-term is the Hopf 
term, 

Stop^of d'xdre'^'^AFjk, (64) 

JSixSixSi 



where Fij = diAj — djAi, i,j,k = x,y,T. The space is 
compacted to 5*1 x 5*1 , and the time is compacted to the 
last ^i. The Hopf term can be written as a total dif- 
ferential locally, e^^^AiFjk — e'^^^dk[fij[Ai, Aj)], where 
fij{Ai,Aj) is a (nonlocal) function of Ai and Aj. Thus, 
at open boundary condition, the integral is determined 
by the boundary values. Stop = •S'l — 5'r, here we have 
cut the space along y-direction. 5'l is given as 

Sl = f dydT{fyr{Ay,Ar)- frv{Ar,Ay)) 

JSixSi 

= 0f dAidA2(/AiA2(^Ai, AaJ - /a2Ai(^A2,4ai)) 

Js[xS{' 

= 61 /" d^Xe"^AjFjK (65) 

where I,J,K = Ai,A2,A3 are parameters of the group 
space of G. S[ is the circle formed by parameter Ai(y), 
and S'l is the circle formed by parameter A2(t). Di is 
the area enclosed by the S'^. Above we have mapped the 
two-dimensional integral on the boundary of space-time 
manifold into a three-dimensional integral on the group 
space of G. Notice that the spatial dimension of the 
boundary is ID, the above topological 6'-term is actually 
an effective WZW term of the boundary. 

Since Eq. (65) is a 3-dimensional integral over the 
group space of G, when reducing to the symmetry group 
G, we need at least four points to span the 3-d space 
S'l X Di: \go),\9o9i),\9o9i92), and \9o9i9293)- Thus we 
can identify Eq. (65) as a 3-cocycle or product of 3- 
cocycles under proper gauge choice 

e*"^^ = vz{9q,9o9i,9q9i92,9o9i929z), 

Here 9o, 9i, 92, 93 are group elements in the symmetry 
group G. 

Above arguments can be generalized to arbitrary d- 
dimension. For example, in (H-3)D, we may have 

Stop = el d^xdre'^'-'^F^.F^A. (66) 
JS1XS3 

The topological 0-term (or term in literature) plays 
important roles in various many-body systems. In ref- 
erence 90, the authors came up with a new method to 
calculate the topological 0-term. 

We have shown that the topological term (or the 9- 
term originating from Berry phase) reduce to cocycles if 
we discretize the space and time. The discrete topologi- 
cal term even exist for discrete groups (they are related to 
the 0-term by the embedding argument mentioned previ- 
ously). Although the discrete topological nonlinear sigma 
models constructed from group cocycle are formally close 
to 0-terms in continuous nonlinear model, they actually 
describe quite different physics. The 0-terms in continu- 
ous nonlinear model ignore the physics at cut-off length- 
scale, which should be very important in general, espe- 
cially for those gapped quantum systems, whose fixed 
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FIG. 21: (Color online) A triangularization of (2+l)D space- 
time where each cube represent five tetrahedrons. The shaded 
area represents a [xt] plane, on which each square represents 
two triangles. 

point actions have zero correlation length and quantum 
fluctuation can be non-smooth at arbitrary energy scale. 
Thus, the discrete topological nonlinear sigma models 
can be regarded as the quantum analogous of ^-terms 
in continuous nonlinear sigma model and help us under- 
stand the nature of symmetry protected topological order 
in interacting systems. 

XII. SPT ORDERS WITH TRANSLATION 
SYMMETRY 

In the above we have discussed bosonic SPT phases 
with on-site symmetry G but no other symmetries. Here 
we would like to stress that when we say a SPT phase 
have no other symmetries, we do mean that the ground 
state wave function of the SPT phase has no other sym- 
metries. In fact the ground state wave function of the 
SPT phase can have some other symmetries. What we 
really mean is that when we deform the Hamiltonian to 
construct phase diagram, the deformed Hamiltonians can 
have no other symmetries. 

In this section, we will discuss the SPT phases with 
both on-site symmetry G and translation symmetry. We 
will use the non-linear cr-model approach to obtain our 
results. We have argued that the d-dimensional SPT 
phases with on-site symmetry G are classified by fixed- 
point non-linear cr-models that contain only a topo- 
logical 0-term constructed from a (1-1- (i)-cocycle in 
H^^'^[G,Ut{1)]- The act ion- amplitude (in imaginary 
time) for such a fixed-point non-linear cr-model is given 
by 

e-/^-^(^.+.)= Yi •^t^,'{9.,9,,-,9,). (67) 

{zj...k} 

When the system has translation symmetry, we can in- 
clude additional topological S-terms which lead to richer 
SPT phases. Let us use (2-|-l)-dimensional systems as 
examples to discuss those addition topological 0-terms. 

When we say a (2-|-l)-dimensional system has a trans- 
lation symmetry, we mean that the system has a discrete 
translation symmetry in the two spatial directions. We 



must choose the triangularization of the space-time in a 
way to be consistent with the discrete spatial translation 
symmetry. In this case, we can include a new topological 

0- term: 

e-lcL = e-lcH^s)^-lc^i.r) (68) 

where 

e-lcH^r)=Yl n i'^Tr"'i9^,9J,9k)■ (69) 
[xt] {ijk}e[xt] 

The translation invariant space-time complex can be 
viewed as formed by many 2-dimensional sheets, say, in 
x-t directions (see Fig. 21). We pick a sheet [xt] in x-t 
directions, then U{ijk}e[xt]('^fy'^"'i9i,9j,9k) is simply a 
topological 0-term on the [xt] sheet constructed from a 
2-cocycle i^f* in ■H^[G, J7t(1)]. In the above expression, 
n{y7c}e[xt] multiply over all triangles in the [xt] sheet and 
Y[[xt] multiply over all the [xt] sheets in the space-time 
complex. 

We can include a similar topological ^-term -1^ ^ ^'^^ ^ 
by considering the sheets in y-t directions and using an- 
other 2-cocycle t'f*. A third topological 6'-term can be 
added by viewing space-time complex as formed by many 

1- dimensional lines in time direction: 

e-/^^(^J)=J| Yl {<^lrH9.,9A (70) 

[t] mm 

Here n{ij}g[t] multiply over all segments in the [t] line 
and Y[[t] multiply over all the [t] lines in the space-time 
complex. In fact Y\^ijy^[t]{i^iY'' {gi, gj) is a topological 
0-term on a single [t] line constructed from a 1-cocycle 

We can also try to include the fourth new topological 
(?-term by considering the sheets in x-y directions: 

,-IC-\.r)^ll II {viyr^-{g,,g,,g,). (71) 

[xy\ {ijk}e[xy] 

But such a topological 6'-term corresponds to a 
LU transformation with a few layers. In fact 
Y\{^]k}(^[xv\i^2'Y'''°{9u9j.9k) is a LU transformation 
when viewed as a time-evolution operator. So there is 
no fourth new topological 0-term. 

We see that SPT phases in (2+l)-dimcnsions with an 
on-site symmetry G and translation symmetry are char- 
acterized by one 1-cocycles vl e ^^[6*, [/t(1)], two 2- 
cocycles :^f*,t'|* g 'H^[G', C/t(1)], and one 3-cocycle g 
H^[G,Ut{1)]- If we believe that those are all the pos- 
sible topological ^-terms, we argue that the SPT phases 
in (2-|-l)-dimensions with an on-site symmetry G and 
translation symmetry are classified by [G, C/t(1)] x 
n^[G,UT{l)] X n^[G,UT{l)] X n^[G,UT{l)]. a spedal 
case of this result with z^a = is discussed in Ref. 51 
where the physical meaning of i^J, i^f *, ^"2* is explained 
in terms of ID representations and projective represen- 
tations of G. Certainly, the above construction can be 
generalized to any dimensions. 
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If we do not have translation symmetry, we can still 
add the new topological 0-terms, such as e~ ^ ^'^^ ^ . But 
in this case, we can combine n [xt] planes in to one. If 
7i^[G, Ut{1)] is finite, the new topological 0-term on the 
combined plane can be trivial if we choose n properly. 
So, we cannot have new topological 0-terms if we do not 
have translation symmetry and if H'^[G, C/t(1)] is finite. 



XIII. SUMMARY 

Since the introduction of topological order in 1989, we 
have been trying to gain a global and systematic un- 
derstanding of topological order. We have made a lot 
of progress in understanding topological orders without 
symmetry in low dimensions. We have used the K- 
matrix to classify all Abelian fractional quantum Hall 
states, and have used string-net condensation^^'^° 
to classify non-chiral topological orders in two spatial di- 
mensions, and have constructed a large class of topolog- 
ical orders in higher dimensions. 

The LU transformations deepen our understanding of 
topological order and link topological orders to patterns 
of long range entanglements.^^ Such a deeper under- 
standing allows us to obtain a systematic description 
of topological orders in 2D fermion systems. The LU 
transformations also allow us to start to understand topo- 
logical order with symmetries. In particular, it allows 
us to classify all gapped quantum phases in one spatial 
dimension. We find that all gapped ID phases are SPT 
phases (SPT phases are gapped quantum phases with cer- 
tain symmetry which can be smoothly connected to the 
same trivial product state if we remove the symmetry). 
In ID, the SPT phases can be classified by 2-cohomology 
classes of the symmetry group. 

In this paper, we try to understand topological order 
with symmetry in higher dimensions. In particular, we 
try to classify SPT phases in higher dimensions. We find 
that distinct SPT phases with on-site symmetry G in c? 
spatial dimensions can be constructed from distinct el- 
ements in (1 -I- (i)-Borel-cohomology classes of the sym- 
metry group G. We summarize our results in table I for 
some simple symmetry groups. 

We have used two approaches to obtain the above re- 
sult: the LU transformations and topological non- linear 
tj-models. We generalized the usual topological 0-term 
and the WZW term in continuous non-linear cr-model 
to the topological 6'-term and the NLL terms in discrete 
non-linear cr-models (with both discrete space-time and 
discrete target space). 

Our results demonstrate how many-body entangle- 
ments interact with symmetry in a simple situation where 
there is no long range entanglements (ie no intrinsic topo- 
logical orders). This may prepare us to study the more 
important and harder problem: how to classify quantum 
states with long range entanglements (ie with intrinsic 
topological orders) and symmetry. Those phases with 
long range entanglements and symmetry are called sym- 



metry enriched topological orders. Also, we approach 
can be modified and generalized to describe/classify 
fermionic SPT phases, through generalizing the group 
cohomology theory to group super-cohomology theory. 
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Appendix A: Making the condition eqn. (12) a local 
condition 

We can make the condition eqn. (12) on C7* a local 
condition. Instead of requiring eqn. (12), we may require 
to satisfy 

{U' (g) [/^+"' ® U'+y (g) U'+''+y){P' ® p') (Al) 
= (P^ (g) ® C/^+" ® ® jji+^+y^^pi ^ pi) 

for certain projection operators and with Trp* = 1. 
Here [/*(g) (g)[/*+^(8)C/*+^+2' and P'''®p^ are matrices 
given by 

7172 7374,7172 73 74 Ai A2A3A4,A; Ai, V J 

and 

(^'®?'')^PilA73A4,ai;3^7^AiX (A3) 
pi 

^ Q2a3a4/3l/33/347 l7274Ai A2A3,a2"3"4/9i/53/3l7i72 74Ai A2A3 ' 

The condition eqn. (12) implies the condition eqn. (Al) 
because, in the canonical form, the states on the 
sites ai, /32, 73, A4 and the states on the sites 
02, as, "4, /?!, /33, Pi, 71, 72, 74, Al, A2, A3 are un- 
entangled (see Fig. 8). 

Appendix B: Representations and projective 
Representations 

Let us consider a group G that may contain anti- 
unitary time reversal transformation. We can divide the 
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group elements into two classes: 

s{g) = 1 or - 1, g&G. (Bl) 

The group elements that contain an odd number of time- 
reversal operations have s{g) — —1 and the group ele- 
ments that contain an even number of time-reversal op- 
erations have s{g) — 1. 

Unitary matrices u{g) form a representation of sym- 
metry group G if 

M(5i)"s(gi)(52) = u{gig2). (B2) 

where Ms(ffi)(52) = u{g2) if s(gi) = 1 and Us(g^){g2) = 
Kff2)]* if s(3i) = -l. 

The above relation is obtained from the following map- 
ping 

g u{g), if s{g) = 1; 5 -> u{g)K, if s{g) = -1. 

(B3) 

Here K is the anti-unitary operator 

Ka = a*K, (B4) 

where o is a complex number. For example, if s{gi) — 
5(52) = ^1 and s{gig2) = 1, we require that 

uig,)Kuig2)K = u{g,g2) (B5) 

which leads to eqn. (B2). 

Matrices u{g) form a projective representation of sym- 
metry group G if 

u{gi)us(g,){g2) ^(^{gi,92)u{gi92), 91,92^0. (B6) 

Here uj{gi,g2) G C/(l) and uj{gi,g2) ^ 0, which is called 
the factor system of the projective representation. The 
associativity requires that 

[u(5l)Ws(gi)(52)]Us(gi32)(53) = u{gl)[u{g2)Us(g^)ig3)]s{g^) 

(B7) 

or 

^{91, 92)i^{gi92, 93)u{gig2g3) 
= ^^''^^'H92, 93)^^(91, 9293)u{gig2g3) (B8) 

Thus, the factor system satisfies 

^''^^'H92,93)oj{9i,9293) ^ ^^{9i,92)i^{gig2,g3), (B9) 

for all (?2,53 G G. If uj{gi, g2) — 1, u{g) reduces to the 
usual linear representation of G. 

A different choice of pre-factor for the representation 
matrices u'{g) = f3{g)u{g) will lead to a different factor 
system aj'{gi,g2): 

^'(91,92) = . ^!'^^^^}, . ^{91, 92)- (BIO) 



We regard u'(f/) and u{g) that differ only by a pre-factor 
as equivalent projective representations and the corre- 
sponding factor systems ui'{gi,g2) and ^(31,(72) as be- 
longing to the same class w. 

Suppose that we have one projective representation 
ui{g) with factor system wi (171, 172) of class loi and an- 
other U2{g) with factor system uj2{gi,g2) of class a;2 7 ob- 
viously ui{g)®U2{g) is a projective presentation with fac- 
tor group wi (51 , 52)1^2 (51 , 52) ■ The corresponding class w 
can be written as a sum wi -|- clI2- Under such an addi- 
tion rule, the equivalence classes of factor systems form 
an Abelian group, which is called the second cohomology 
group of G and denoted as H'^[G, Ut{1)\. The "zero" el- 
ement e 'H^[G, is the class that corresponds to 
the linear representation of the group. The best known 
example of projective representation is the spin-1/2 rep- 
resentation of SO{?>). The integer spins correspond to 
the linear representations of 50(3). 



Appendix C: ID representations and projective 
representations of (7(1) x and Uil) xi Z2 

In this section, we are going to discuss the ID repre- 
sentations and projective representations of four groups 
U{1) X Z2, [/(I) X ^2, U{1) X Zj, and U{1) x Z'^ . As 
group, Z2 and are actually the same group. However, 
the generator t of Z2 corresponds to a usual symmetry 
transformation which has a unitary representation. The 
generator T of corresponds to the time reversal trans- 
formation which has a anti-unitary representation. 

Let Ue, d E [0, 27r), be an element in U{1). The four 
groups are defined by the following relations 

U{1) x Z2 : tUe = Ugt; 

C/(l) X : TUe = UgT; 

U{1) X Z2 : tUg = U-gt; 

U{1) X : TUg = U-gT. (CI) 

Their representations are given by matrix functions 
M{Ug) and M{t) (or MiT)K). 

A ID representation of C/(l) x Z2 has a form 

Af([/e) = e"'^ M{t)^ri = ±l, (C2) 

where n e Z. One can check that M{T)M{Ug) = 
M{Ug)M{T), for any n. So the ID representation of 
U{1) X Z2 is labeled by n and 77 (or by Z x Z2). 

A ID representation of U{1) x Z2 also has a 
form eqn. (C2) One can check that M{T)M{Ug) = 
M{U^g)M{T) only when n = 0. So there are two ID 
representation of U{1) x Z2 labeled by Z2 (correspond- 
ing to 77 = il)- 

A ID representation of U{1) x has a form 

M{Ug) = e"'^ M(T)ii: = e'-^iC, (C3) 
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where n e Z and e R. Note that M{T)KM{T)K = 1 
for any </>. We find 

M{T)KM{Ue) = M{T)Ke''^ 

= e-"''^M{T)K. (C4) 

Thus M{T)KM{Ue) = M{Ue)M{T)K only when n = 
0. Also under an unitary transformation e"^, M{T)K 
transforms as 

e'*(e''^ii')e-''^ = e'^^+^-^^i^. (C5) 

So different (f> correspond to the same ID representation. 
Therefore, there is only one ID representation for C/(l) x 

zl. 

From the above calculation, we note that M{T)K\n) cx 
I — n) where \n) is an eigenstate of M{Ug): M([/g)|n) = 
e"'^|n). So the group x U{1) describes the symmetry 
group of a spin system with time reversal and 5*2 spin 
rotation symmetry (the U{1) symmetry). 

A ID representation of U{1) x Zj also has a form 
eqn. (C3). We find 

M{T)KM{Ug) = M{T)Ke'^''> 
= e-'^^'>M{T)K = M{U-e)M{T)K. (C6) 

Thus M{T)KM{Ue) = M{U-e)M{T)K for any n e Z, 
and the ID representations for U{\) x Z2 are labeled by 
Z. 

The above relation (C6) also allows us to show 
M{T)K\n) cx where |n) is an eigenstate of M{Ue): 
M{Ug)\n) = e^'^jn). This is the expected transforma- 
tion of time reversal for boson systems, where n is the 
boson number. Therefore U{1) x is the symmetry 
group of boson systems with time reversal symmetry and 
boson number conservation. 

Next, let us discuss the projective representations of 
the four groups. First, let us consider U{1) x Z2, whose 
projective representations may have a form 

MiUg) ^ (^""J' , Af(T) - (? J) • (C7) 

One can check 

/ n „mie\ 

MiT)MiUe) ^ i^^^.g \ j, 

M{Ug)M{T) = (^^„% ^"J'^ . (C8) 

M{T)M{Ug) and M{Ug)M{T) differ by a total phase 
only when m = n G Z, in that case M{T)M{Ug) = 
M{Ug)M{T). Note that M{T)M{T) = 1. So we have a 
trivial projective representation. If we choose M{T) = 

, we will have M{T)M{T) = -1. But we still have 

a trivial projective representations, since if we add an 
phase factor M{T) = iM(r), we have M{T)M{T) = -1 
Thus U{\) X Z2 has only one trivial class of projective 
representations. 



Second, let us consider the projective representations 
for C/(l) X Z2, which may have a form 

MiUo) - (^'7 Ji,^ , M{T) - (J J) • (C9) 
One can check 

/ n mig\ 

M{T)M{Ug)=i^J,e j, 
M{U^g)MiT} = (^J^,, ''^ . (CIO) 
We have 

M{T)M{Ug) = e-("+")"'Af(C/_e)M(r) (Cll) 

Also note that M{T)M{T) = 1. So we have a non- 
trivial projective representation. If we add a phase factor 
M{Ug) = e'''^M{Ug), we will have 

M{T)M{Ug) = e'^''+"'-'^''y''^ M[U-g)M{T) (C12) 

So the projective representations with different n and m 
belong to two classes m+n = even and m+n = odd. Thus 
U{1) X Z2 has two classes of projective representations 
labeled by Z2. 

The projective representation of U{1) x can have a 
form 



Ug ^ M{Ue) = { 
T M{T)K = 




(C13) 



One can check 



M{T)KM{Ug) = 
M{Ug)M(T)K = 




(CM) 



We have 

M{T)KM{Ug) = c-^"+"''>''^M{Ug)M{T)K (C15) 

Note that M{T)KM{T)K = 1. So we have a projective 
representation when m^n ^ J. . If we add a phase factor 
M{Ug) = e'''H.{[Ug), then 

M{T)KM{Ug) = e-("+™+2''')'*M(J7e)M(T)ivr. (C16) 

So the above projective representations belong to two 
classes: m + n = even and m + n = odd. 

The projective representation of x U{1) may also 
have a form 

M{Ug) = (^'q ' Ji,^ , M{T)K = K. 

(C17) 



One can check 



M{T)KM{Ue) = 
M{Ue)M{T)K = 




(C18) 



We have 

M{T)KM{Ug) ^ e-'^'^+"'^'Hl{Ug)M{T)K (C19) 

Note that M{T)KM{T)K = -1. So we also have a pro- 
jective representation when m^n ^ J. . Those projective 
representations also belong to two classes: m-\-n — even 
and 771 + 71 = odd. So U{1) x has four classes of 
projective representations labeled by Z2 x Z2. 

For the U{\) x group, its projective representation 
may have a form 

M{Ue) = (^^^' ^„L^ , M{T)K = {^^ '^'^ K. 

(C20) 

One can check 



M{T)KM{Ue) = 
M{U^g)M{T)K = 




(C21) 



M{T)KM{Ub) and M{Ue)M{T)K differ by a total phase 
only when m = n € Z. Note that M{T)KM{T)K = -1. 
If we add a phase factor M{T) = e">'M{T), we still have 
M{T)KM{T)K = -1 So we have a non-trivial projec- 
tive representation. Those projective representations for 
different n = 777 all belong to one class. Thus C/(l) x 
has two classes of projective representations (the one dis- 
cussed above plus the trivial one) labeled by Z2. 

Appendix D: Group cohomology 

The above discussion on the factor system of a projec- 
tive representation can be generalized which give rise to 
a cohomology theory of group. In this section, we will 
briefly describe the group cohomology theory.^^ 

1. G- module 

For a group G, let M be a G-module, which is an 
Abelian group (with multiplication operation) on which 
G acts compatibly with the multiplication operation 
(ie the Abelian group structure): 

g-{ab) = {g-a){g-b), g e G, a,beM. (Dl) 

For the most cases studied in this paper, M is simply 
the U{1) group and a an t7(l) phase. The multiplica- 
tion operation ah is the usual multiplication of the C/(l) 
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phases. The group action is trivial: g ■ a = a, g G G, 
a e U{1). We will denote such a trivial G-module as 
M = U{1). 

For a group G that contain time-reversal operation, 
we can define a non-trivial G-module which is denoted 
as Ut{1)- Ut{1) is also a U{1) group whose elements 
are the U{1) phases. The multiplication operation ah, 
a,b G Ut{1), is still the usual multiplication of the U{1) 
phases. However, the group action is non-trivial now: 
g ■ a = a^(f), g e G, a & C/t(1), where s{g) = 1 if g 
contains no anti-unitary time reversal transformation T 
and s(g) — —1 if g contains one anti-unitary time reversal 
transformation T. 

The module defined above is actually a model over a 
ring Z, since we have the following operation Z x Af — >■ M: 

Vn e Z,Va e M, a" G M. (D2) 

A model M can be over a more general ring R if we have 
the operation R x M ^ M: 

WneyaeM, a''eM, (D3) 

such that 

a"6" = (a6)", a^a"* = a'"" = (a"*)", a^" = a, 

(D4) 

if R has multiplicative identity l/j. 

Such a general concept of a module over a ring is a 
generalization of the notion of vector space, wherein the 
corresponding scalars are allowed to lie in an arbitrary 
ring. As we have seen, modules also generalize the notion 
of Abelian groups, which are modules over the ring of 
integers 

2. Algebraic definition of group cohomology 

Let LUn{giT--j g-n) be a function of n group elements 
whose value is in the G-module M. In other words, 
w„ : G" M. Let C"(G,M) = {w„} be the space of 
all such functions. Note that C'^{G,M) is an Abelian 
group under the function multiplication uj[[{gi, g^) = 
uJnigi,.:,gnWni9i,-,9n)- We define a map (i„ from 
C"[G,C/t(1)] to C"+i[G,i7T(l)]: 

{dnUJn){gi, ..■,5n+l) = 

[51 ■ ^nig2, ■■; gn+i)](^'k^^^"^\gi, gn)'x 

71 

) (D5) 

Let 

S"(G,M) = {c^„|^„ = d„_ic^„_i|L^„_i e C"-i(G,M)} 

(D6) 

and 

Z"(G, A/) = {w„M„t.„ = l,c^„ e C"(G,A/)} (D7) 
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B"'{G,M) and ^"(G',M) are also Abelian groups which 
satisfy B''{G,M) c Z"(G,M) where B^{G,M) = {!}. 
The n-cocycle of G is defined as 

n'^{G,M) = Z'\G,M)/B'\G,M) (D8) 

Let us discuss some examples. We choose M — Ut{1) 
and G acts as: g ■ a = a^^^\ g £ G, a £ Ut{1)- In this 
case ujnigi, ■■■,gn) is just a phase factor. From 



(D9) 



we see that 

Z°[G, Ut{1)] = WoH^''^ = oJo} ^ U^fil) (DIO) 

If G contain time reversal, U^{1) ~ {1,-1}. If G 
does not contain time reversal, Uq^{l) — U{1). Since 
6°[G,L/t(1)] ^ {1} is trivial, we obtain 'H°[G,UTil)] = 
C/G(l). 
From 

(diwi)(gi,g2) = a;f^'^(52)wi(5i)/a;i(5ig2) (DH) 
we see that 

Z'[G,Ut{1)] = {c.i|c.i(5iK(''^^(.92) = wi(gi52)}. 

(D12) 



Also 



B'[G, Uril)] - W^M,) = c.^^«^V^o} (D13) 

n^[G,UT{l)] = Z^[G,Ut{1)]/B^[G,Ut{1)] is the set of 
all the inequivalent ID representations of G. 
From 

(d2W2)(gi,52,53) (D14) 

= ^^2 ' (32 , ff3 ) W2 (ffl , 9253 ) /'^2 (51 52 , 53 ) W2 (gi , 32 ) 

we see that 



Z^[G,Ut{1)]^{lo2\ 



(D15) 



^^2(5i,5253)'^2^^'^(52,53) = ^^2(5152, 53)^2(51, 52)}- 



and 



From 

(d3W3)(5i,52,53,54) 

^ ^3 ^ (52 , 53 , 54 )a^3 (51 . 5253 , 54 )^3 (51 . 52 , 53 ) 
^^3 (51 52 , 53 , 54 ) ^^3 (51 , 52 , .93 54 ) 

we see that 

Z3[G,C/t(1)] = {C.3| 



(D17) 



(D18) 



^3 (92,93, 94)1^3 (51 1 5253 , 54)^^3 (51 , 52 , 53 ) ^ -j^i 
W3 (5152 , 53 , 54)^3 (51 , 52 , 5354) 



and 



S'[G,C/t(1)] 



(D19) 



6^[G,t/T(l)] = {0.21^2(51,52) =^i''''(52Vi(5i)M(5i52)} 

(D16) 

The 2-cohomology group 'H'^[G, C/t(1)] — 
Z^[G,Ut{1)]/B^[G,Ut{1)] classify the projective 
representations discussed in section B. 



r I / X ^^2 (52,. 93)^2 51,5253 K 

= {W3 1^3(51, 52, 53) = ^ ^ — 7 ^1, 

^2(5152, 53)^^2(51, 52) 

which give us the 3-cohomology group H'^[G, [/^(l)] = 
Z^[G,UTil)]/B^[G,UTil)]. 

In this paper, we will show that H^^'*[G, C/t(1)] can 
classify SPT phases in d-spatial dimensions with an on- 
site unitary symmetry group G. Here the on-site sym- 
metry group G may contain time-reversal operations. 



3. Geometric interpretation of group cohomology 

In the following, we will describe a geometric interpre- 
tation of group cohomology. First, let us introduce the 
map i/„ : G"+^ — >■ M that satisfy 

9 ■ '^n{90,9l,--,9n) = l^n{990,99l,-;99n), (D20) 

for any g G G. We will call such a map i>n a n-cochain: 

C"{G,M) = {vr,\g ■ Vn{go, ■■■,gn) = J^«(55o, ■••,55n)}- 

(D21) 

LOn discussed above is one-to-one related to through 

w„(5i, ■■■,9n) = Vn{^,9i,9i92, •■•,5i ■■■9n) 

= j/„(l, 51,52, ...,5„) (D22) 

where 5^ = 5i52 ■ • •5i- 

We can rewrite the dn map, d„ : a;„ — uJn+i, as d„ : 

— >• i^n+l- 
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(d„j/„)(l,5i,5i52, ••■,51 ■ ■ -Sii+i) 
= 91 ■ i^n(l,52,52g3, •••,52 • ■ ■ gn+i^n^^^ ^ (1, 31,5132, •••,3i ' ■■ 9n)v~^ O-, 9x92, 9x9293, . . . , 5i • • • g„) X 

Vn{l, 91, 919293, •••,51 • • •ffn)^^- 

= Vn{9l, 9192,919293, ■■■,9192 ■ ■ ■ 9n+l)'^i~^''"^\i, 91, 9l92, ■-, 91 ■ ■■ 9n)Vn^ {I, 9192,919293, ■■■,9l ' ■■9n)x 

Vn{^,9i, 919293, ■-,91 ■ ■■9n)--- 

= Vn{h,h,h,--,9n+l)v^^^''^^\'^,h,h,---,9n)Vn^{^, 92,93, • • • , ( 1 , ffl , 53 , ■■■,gn)--- (D23) 

I 



The above can be rewritten as (after the renaming gi — >■ 

9i) 

{dnVn){90,9l, ■■■,9n+l) 
n+l 

= '[\vl:^^'^\go, ..,gi_i,g^+i, ...,gn+i) (D24) 

which is a more compact and a nicer expression of the dn 
operation. 

When n = 1, we have 

[divi){ga,gi,g2) = J^i (51,52) i'i(ffo,5i)/i^i (50,52) (D25) 
For n = 2: 

^^2(51, 52, 53)1^2(50, 51, 53) 



{d2V2){9G,9l,92,93) = 



and for n = 3: 



^^2(50, 52, 53)1^2(50, 51, 52) 



(D26) 



{d2,V2,){.9o,9i,92,93,9i) (D27) 
^ ^^3(51, 52, 53, 54)1^3(50, 51, 53, 54)^^3(50, 51, 52, 53) 
1^3(50,52,53,54)^^3(50,51,52,54) 

We may represent the 1-cochain, 2-cochain, and 3- 
cochain graphicahy by a hne, a triangle, and a tetra- 
hedron with a branching structure respectively (see Fig. 
10). We note that, for example, when we use a tetrahe- 
dron with a branching structure to represent a 3-cochain 
^3(50,51,52,53), the last variable 53 is at the vertex with 
all the edges point to the vertex (see Fig. 10b). After 
removing the 53 vertex and the connected edges, (72 is 
at the vertex with all the remaining edges point to the 
vertex (see Fig. 10b). This can be repeated. We see 
that a tetrahedron with a branching structure gives rise 
to a natural order 50, 51,52, 53 • In general, a c?-cochain 
can be represented by a c?-dimensional simplex with a 
branching structure. We also note that a d-dimensional 
simplex with a branching structure can have two different 
chiralities (see Fig. 24). The simplex with one chirality 
correspond to and the simplex with the other chirality 
correspond to v'^'^ (see eqn. (El)). 

In this way, we obtain a graphical representation of 
eqn. (D25) and eqn. (D26) as in Fig. 10. In the graphical 
representation, eqn. (18) implies that the value of a 1- 
cocycle vi on the closed loop (such as a triangle) is 1 and 



eqn. (19) implies that the value of a 2-cocycle V2 on the 
closed surface (such as a tetrahedron) is 1. 

Let us choose M = U{1) and consider a 1-form Vli on 
the plan in Fig. 10a. Then the differential form expres- 
sion 

r r9i rg2 rg2 

/ dfii = / fii- ni+ ni (D28) 

•'(90:91,32) J go J go J gi 

give us eqn. (D25) if we set 



(c?ii^i)(5o,5i,52) exp( i / d^i 

'(90,91.92) 



and 



vi{9i,9]) = exp(^i j 17i^ 



(D29) 



(D30) 



Here L ^ is the integration on the triangle 

•'(90,91,92) *= *= 

(50,51,52) in Fig, 10a. Similarly the differential form 
expression 



dO, 



(90,91,92,93) 



(91,92,93) 



(90,92,93) 



'(90,91,93) 

give us eqn. (D26) if we set 
{d2V2){9o,9i,92,93) = exp(^i 



^2 (D31) 



(90,91,92) 



(90,91,92,93) 



and 



M9t,93,9k) expM / ^2 
J(gi,gj,gk) 



(D32) 



(D33) 



This leads to a geometric picture of group cohomology. 
For example, if 0,2 is a closed form, dQ,2 = 0, the corre- 
sponding 1^2(91, 9 j, 9k) will be a cocycle. If 5^2 is an exact 
form, D,2 = dfii, the corresponding 1^2(91, 9j, 9k) will be 
a coboundary. 



4. Cohomology on symmetric space 

We would like to mention that cohomology can also be 
defined on symmetric space G/H where H is a subgroup 
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of G. However, cocycles on the symmetric space G/H 
can also be viewed as cocycles on the group space G (the 
maximal symmetric space) and we have Z'^{G/H,M) C 
Z'^{G,M) . As a result, the SPT phases described by 
the quantized topological 0-terms on the symmetric space 
G/H can all be described by the quantized topological 9- 
terms on the maximal symmetric space G. So classifying 
quantized topological 0-terms on the maximal symmetric 
space G lead to a classification of all SPT phases. 



Appendix E: Branching structure of a complex 
1. Branched simplex and its geometric meaning 

In geometry, a simplex is a generalization of the no- 
tion of a triangle or tetrahedron to arbitrary dimensions. 
Specifically, an n-simplex is an n-dimensional polytope 
which is the convex hull of its n + 1 vertices. It can also 
be viewed as a complete graph of its n+1 vertices. For ex- 
ample, a 2-simplex is a triangle, a 3-simplex is a tetrahe- 
dron, and a 4-simplex is a pentachoron. An n-simplex is 
the fundamental unit cell of n-manifolds, any n-manifold 
can be divided into a set of n-simplexes through the stan- 
dard triangulation procedure. It is obvious that any in- 
variant under the re-triangulation of n-manifolds would 
automatically be a topological invariant. 

One of such examples is the famous state sum in- 
variants of 3-manifolds first proposed by Turaev and 
Viro^^. The basic idea in their construction is associat- 
ing a special data set(e.g., 6j-symbol) with each tetrahe- 
dron, and then showing the states sum invariants under 
re-triangulations. However, their construction requires a 
very high tetrahedral symmetry for the data set, based 
on the assumption that all the vertices/edges/faces in a 
tetrahedron are indistinguishable. Indeed, in a more gen- 
eral set up, labehng the vertices/edges/faces is important 
because they are actually distinguishable objects. 

A nice local scheme to label an n-simplex is given by 
a branching structure. A branching is a choice of an 
orientation of each edge of an n-simplex such that there 
is no oriented loop on any triangle. For example. Fig. 
22 (a) is a branched 2-simplex and (c) is a branched 3- 
simplex. However, (b) is not allowed because all its three 
edges contain the same orientations and thus form an 
oriented loop, (d) is also not allowed because one of its 
triangle contains an oriented loop. Actually, a consistent 
branched triangulation can always be induced by a global 
labeling of the vertices (We notice any labeling of the 
vertex u*, i = 0, 1, 2, • • • , w" will imply a nature ordering 
-y' < a i < j) . This is because any global ordering will 
induce a consistent local ordering for all the triangles of 
an n-simplex. If we associate an orientation from i to j if 

< , it is obvious that there will be no oriented loop 
on any triangle. 

A branched n-simplex will have the following proper- 
ties: 

(a). Any given branching structure for an n-simplex 




(a) ib) 




FIG. 22: Examples of allowed ((a),(c)) and unallowed 
((b),(d)) branching for a 2-simplex and a 3-simplex. 




(a) {b) 

FIG. 23: (a): If a 3-simplex contains a vertex with no incom- 
ing edge, we can label this vertex as u° and canonically label 
the vertices of the remaining 2-simplex as v^,v'^,v"*. Such a 
scheme can be applied for arbitrary n-simplex if n — 1-simplex 
has a canonical label, (b): If a 3-simplex contains no ver- 
tex without incoming edge, then there must be a vertex with 
one incoming edge. (Because canonical ordering is true for 2- 
simplex.)If we label this vertex as v'^, the vertex connect to 
through a incoming edge must contain no incoming edge, 
otherwise the branching rule will be violated. The above argu- 
ment is true for n-simplex if n — 1 simplex can be canonically 
ordered. 



will uniquely determine a canonical ordering of the ver- 
tices. For example. Fig. 22 (a) is a branched 2-simplex 
with three vertices, one of them contains no incoming 
edges, one of them contains one incoming edge and the 
rest of them contains two incoming edges. Thus, we can 
canonically identify the vertex corresponding to each d', 
i — 0,1, 2. Such a canonical labeling scheme can be ap- 
plied to any n-simplex, due to the fact that the n -I- 1 
vertices of any n-simplex will be uniquely associated with 
0, 1, 2, • • • , n incoming edges. 

Proof. Assuming the above statement is true for an n- 
simplex (The statement is true when n = 2, see Fig. 
22.), let us prove it is also true for (n -I- l)-simplex. See 
Fig. 23, if the (n -I- l)-simplex contains a vertex with no 
incoming edge, we can drop this vertex and apply the 
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is locally identical to their orientations in (a). 



FIG. 24: (Color online) (a) and (b): 2-simplex has two dif- 
ferent chiralities, depending on the clockwise or anticlockwise 
ordering of the vertices, (c) and (d): The chirality of the 3- 
simplex can be determined by the chirality of the 2-simplex 
which is opposite to the vertex u". Similarly, the chirality of 
n-simplex can be determined by the chirality of n — 1 simplex 
which is opposite to u". 

statement for the remaining n-simplcx. If we label the 
n + 1 vertices of the n-simplex as 1, 2, • • • ,n+ 1, it is 
clear the vertex with no incoming edge can be labeled 
as 0. In the following we will prove a branched {n + 1)- 
simplex must contain a vertex with no incoming edge. 
If an (n -I- l)-simplex does not contain any vertex with 
no incoming edge, it must contain a vertex with one 
incoming edge. This is because if we remove an arbitrary 
vertex (denoted as vq) of the rt-|-l-simplex, the statement 
is true for the rest n-simplex. Hence we can always find a 
vertex with one incoming edge. Let us denote this vertex 
as f ^ and it is clear that the orientation of the edge that 
connects v'^,v^ must be outgoing towards u°(Otherwise 
is a vertex with no incoming edge). However, in this 
case, the edges that connect and other vertices must 
be outgoing from w*^, if the branching rule is not violated. 
Thus, is the vertex with no incoming edge. □ 

(b) Although the branching rule of n-simplex uniquely 
determines the ordering of the vertices, it could not 
uniquely determine an n-simplex. This is simply be- 
cause the mirror image of a branched n-simplex is also 
a branched n-simplex with the same vertices ordering. 
Thus, any branched n-simplex has a unique chirality ±1. 

Proof. It is clear that 2-simplex has two different chiral- 
ities (see Fig. 24). Assuming that an n — 1-simplex has 
a unique chirality, let us deform the boundary of an n- 
simplex (which can be divided into n— 1-simplex) into an 



n — 1 sphere. Due to the fact that there is one and only 
one vertex v*^ of an n-simplex without incoming edges, 
we can make a canonical convention and determine the 
chirality of the n simplex by the chirality of the n — 1 
simplex opposite to v'^ . Such a definition is sufRcient be- 
cause mirror reflection will always change the chirality of 
the boundary of any n-simplex. Indeed, we can define 
the chirality of the n simplex by the chirality of the n — 1 
simplex opposite to any up to a global sign ambiguity 
(e.g., reversing the chiralities for all n-simplices.). □ 

The above two properties allow us to use the branched 
n simplex to represent an n-cocycle: 

I'n' ''{gt,gj, ■ ■ ■ ,gk), (El) 

where Sij,,,k = il are determined by the chirality of the 
simplex and gi, gj, - ■ ■ , gk are defined on the canonically 
ordered vertices v^,v^, - ■ ■ , w". 

Finally, let us briefly mention the geometric meaning 
of the branched tetrahedron in 3 dimensions. See in Fig. 
25 (a), in a dual picture, the orientations of the edges of 
tetrahedron correspond to the orientations of the region 
of the simple polyhedron. A branching on a simple poly- 
hedron allows us to smoothen its singularities and equip 
it with a smooth structure as shown in Fig. 25 (b). At 
a more rough level, it can be shown that branched tetra- 
hedron can be used to represent the Spin'^-structures on 
the ambient manifolds^^. 



2. Basic moves 

To show the topological invariance of the amplitude: 

^=^^I[''^'"^''(9^,gJr■■ ,gk), (E2) 
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we need to generalize the Turaev-Viro moves to their 
branched versions in arbitrary dimensions. Because each 
move will have many different branched versions, it is not 
easy to check all the branched versions case by case. In 
the following, we will introduce a simple way to look at 
the basic moves. 



a. Graphic representation of {dni'n){go, Qi, 
basic moves 



and 



In last section we have shown a branched n- 
simplex can represent an rt-cocycle {gi,gj, ■ ■ ■ , gk ) 
or its inverse v~^{gi,gj, ■ ■ ■ ,gk), depending on the 
chirality of the branched n-simplex. Here we 
want to show the boundary of a branched n + 
1-simplex can represent (d„t'„)(5o, 5i: ■ ■ • j 5n+i) = 



n 



i=0 "'^ 



(50, • ■ ■ ,gi-i,gi+i, ■ ■ ■ ,gn+i)- Since any n - 



1-simplex branched simplex has a canonical ordering for 
its n + 2 vertices and its boundary contains n + 1 n- 
simplices (We can label these rt + 1 n-simplices as Sn{vi), 
where Vi is the vertex opposite to the n-simplex.), it is 
not surprising if we use the n-simplex Sn{vi) to represent 
i^nigo, ■ ■ ■ ,gi-i,gi+i: ■ ■ ■ ,gn+i) or its inverse. However, 
the key difficulty is that we need to show that the chi- 
rality of the n-simplex Sn{vi) is determined by ±(—1)*, 
where the global sign ± depends on the chirality of the 
n -|- 1-simplex. 

Proof. It is easy to check that the above statement is true 
for n ~ 2. Thus, we can represent (d2;^2)(5o, 5i, 32, gs) as 
a branched tetrahedron. Its boundary 2-simplex 6*2 (w*) 
has opposite chirality for even and odd i. If we assume 
the above statement is true for n — 1, let us proof it is 
also true for n. First let us remove the vertex v'^ from the 
n+ 1-simplex that represents (d„t'„)(go) 5ij • • ' ^gn+i)- 
By applying the statement to the rest n-simplex, whose 
boundary contains n n — 1-simplices {i — 

1, 2, • • • , n) with chirality ±(— 1)\ However, according to 
the definition, the chirality of any 5'„(w') (i = 1, 2, • • • , n) 
can be defined by the Sn-iiv^) simplex which is opposite 
to v^, thus we prove S'„(f *) (i = 1, 2, • • • , n) will also have 
opposite chiralities for even and odd i. To prove that the 
above statement is also true for Sn{v'^), we can remove 
any vertex j ^ and apply the same scheme. Although 
there can be a global sign ambiguity for the chirality of 
any n-simplex S'„(u*) with i ^ j, it is sufficient to show 
Sn{v^) (* = 0, 1, • • • , n) will have opposite chirahties for 
even and odd i. Thus, S'„(u*) (i = 0, 1, • • • ,n) will have 
chirality ±(—1)' with the global sign ± determined by 
the chirality of the n -|- 1-simplex. □ 

Based on graphic representation of 
(d„i/„)(c/o,5i: • • • ,gn+i), it is easy to check that all 
the basic moves are actually induced by the identity: 



,(-1)" 



{go, 



(E3) 
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gl g 




FIG. 26: (Color online) (d2!^2)(go, <7i, 52, 33) can be repre- 
sented as the boundary of a 3-simplex. (a) and (b) correspond 
to two different basic moves of 2-simplexes. 



Due to this identity, the chirality of the n -I- 1-simplex 
is not important because if we inverse both sides we will 
end up with the same identity. Thus, we can pick up 
any n -I- 1-simplex to represent (d„^'„)(go,5i, • ' ' ,5«+i) 
and project it into the n-plane from opposite directions. 
The shadows of these two projections are n-manifold with 
exact the same vertices. However, they may correspond 
to different ways of triangulations. Thus, each side of the 
equations of the basic moves will correspond to the two 
different ways of projection. For example. Fig. 26 (a) 
represents 2 -H> 2 moves: 



j^2(5o,5i,53)j^2(gi,g2,g3) 



'^2(go,5i,g2)i^2(go,52,53) 

(E4) 



and Fig. 26 represents 1 -H> 3: 

V2{gi,g2,gz) = Mao^ 9i, g2)'^2{go, g2, g3)'^2^ (go, gi, ga) 

(E5) 

However, all these two equations will be equivalent to the 
identity: 

1^2(51, 52, 53)^^2^^(30,52, 53)^^2(50, 51,53)^^1^2(50, 51,52) = 1 

(E6) 

We also notice the projection from opposite directions 
will induce opposite chiralities for the boundary of the 3- 
simplex, that's why we need to change the chiralities of 
the 2-simplex in one side of basic moves. Such a change 
corresponds to inverse 1^2 when we move it from left side 
to right side of Eq. (E3) , which is consistent with multi- 
plication rules of the complex number i>2- It is also clear 
that the above two different moves correspond to projec- 
tions in different ways, hence all of them are equivalent 
to the above identity. Similar argument is true for arbi- 
trary dimensions. In conclusion, the identity Eq. (E3) 
will induce the correct 2 -s-^ n and 1 -H- n + 1 moves in n 
dimensions. 
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FIG. 27: (Color online) Examples of unallowed branched 
moves in 2D and 3D. 



b. Some final details 

It looks like we have successfully generalized the basic 
moves to their branched versions, however, there are still 
some subtle issues here, especially in 2D and 3D. This 
is because the n + 1-simplex representation of the ba- 
sic moves relies on the assumption that any basic moves 
can be associated with a consistent branching structure, 
which is not generically true in 2D and 3D. Fig. 27 (a) 
is an example of unallowed 3 — > 1 move in 2D and (b) 
is an example of unallowed 3 — > 2 move in 3D. Although 
a global labeling scheme will not allow triangulations to 
contain any local pieces like the left part of Fig. 27, how- 
ever, those local configurations can be generated dur- 
ing local moves because each simplex still satisfies the 
branching rule. In this case, we can not directly apply 
these unallowed local moves to the local pieces. For- 
tunately, in 2D and 3D it has been proved^^ that any 
branched triangulations can still be connected through 
all the allowed moves. In high dimensions, we can show 
that there are no unallowed moves like these. In the fol- 
lowing let us prove this statement. 

We notice that the 2 — >■ n move can be realized by 
adding one more edge while 1 — )• n + 1 move can be 
realized by adding one more vertex and n + 1 edges, let 
us show that these two moves are always allowed, by 
adding proper orientation(s) on the edge(s). 

Proof. It is trivial to show that the 1 — > n + 1 move 
is always allowed, by adding vertex without incoming 
edge. The existence of 2 — n move can be slightly more 
complicated. One can easily check that it is true when 
n = 2. Let us assume that it is true for 1-simplcx now. 
We label the two unconnected vertices as w", v'^ and label 
other vertices as w'. (Notice we don't require a,b,i have 
an ordering here.) To show there always exist a proper 
orientation for the edge ab, we only need to show that 



any triangles made by v"', and v'^ will not violate the 
branching rule. If there exists a vertex v'' containing two 
incoming edges from v°',v^ or containing two outgoing 
edges towards v"" , we can remove this vertex and apply 
the statement to the rest n — 1 simplex. It does not 
matter what's the orientation on ab, the triangle abi will 
not violate branching rule. If such a vertex does not exist, 
we can show that there can be only two cases: either 
contains no incoming edge, contains no outgoing edge 
or the opposite case. In both cases, we can find a proper 
orientation for the edge ab. □ 

The inverse of the above moves, namely, the n — > 2 
move can be realized by removing one edge, and the n + 
1 — >■ 1 move can be realized by removing one vertex. 
Now let us show that these moves are always possible 
when n > 3. 

Proof. To prove that these moves are always possible in 
dimensions n > 3, let us understand why sometimes they 
are impossible in 2D and 3D. Actually, this is simply be- 
cause three edges of an oriented triangle which violates 
the branching rule can belong to three different simplices 
before we apply 3 — ?■ 1 or 3 — )■ 2 move. However, after 
we apply the move, they belong to the same triangle and 
hence violate the branching rule. In high dimensions, 
when we apply these inverse move, we always start from 
a complete graph and the number of simplices is always 
larger than 3. Thus, any triangle must belong to one of 
the n-simplex and will not violating the branching rule. 
If there is no triangle violate the branching rule in a com- 
plete graph, of course there will be no triangle violating 
the branching rule by removing edge or vertex. □ 

Appendix F: (1-|-1)D solutions of eqn. (12) 

1. U^[g) is a linear representation 

To show that U^{g) defined in eqn. (27) is a linear 
representation of G, let us compare the combined actions 
of W^ig) and U^{g'g~^) with the action of U^{g') which 
are given by (see Fig. 28) 

U'{g'g-')U\g)\ai,a2) (Fl) 
= f2{ai,a2,g,g*)f2igai,ga2,g'g~'^,g*)\g'ai,g'a2) 

and 

[/*(g')|ai,a2) ^ f2{ai,a2,g',g*)\g'ai,g'a2) (F2) 
We see that 

/2 («! , "2 , 5, 3* )/2 7 5^2 , g'g^ ^,g*)f2^{ai,a2,g' ,g*) 
^ '^2{ai,g'^g*,g*) V2{gai,gg'~^g* ,g*) V2{a2,g'~^g* ,g*) 

'^2{a2,g'^g*,g*) V2{ga2,gg'~^g* ,g*) V2{ai,g'~^g* ,g*) 

^ V2{ai,g^^g* ,g*) V2{ai,g'^^g* ,g^^g*) 1^2(^2, g'^^g* , g*) 

V2{a2,g-^g*,g*) 1^2(0^2, g'^^g*,g^^g*) V2{ai,g'-^g* ,g*) 

(F3) 




FIG. 28: (Color online) The evaluation of the 2-cocycle 
on the above two complexes with branching structure gives 
rise to two phase factors in eqn. (Fl) and eqn. (F2), which 
shows that the ratio of the two factors, eqn. (F3), is equal to 
1, since the complexes in (a) and (b) overlap. 

a p P y y a 



FIG. 29: (Color online) The ID state eqn. (25) on a ring. The 
degrees of freedom form maximally entangled dimer states. 

The above expression can be represented as Fig. 28 which 
indicates that the expression is equal to 1. Thus C/*(g) 
defined in eqn. (26) form a unitary representation of G. 



2. U\g) satisfies eqn. (12) 

The action of ^U'^{g) on the ID state on a ring in Fig. 
29 is given by 

(g>iU'ig)\a,p;P,r,l,a) (F4) 
= f2{a,f3,g,g*)f2{l3,-f,g,g*)f2{j,a,g,g*)x 
\9a,gl3;g/3,g^;g'y,ga). 

From (27), we see that 

f2{a,f3,g,g*)f2{f3,'y,g,g*)f2{'l,a,g,g*) 

^ V2{a,g~^g* ,g*) v2{P,g~'^9* ,g*) '^2(7,9'^ 9*, 9*) 
'^2{l3,9~'^9*,g*) '^2{7,9~'^9*,9*) ^2{a,g-'^g* ,g*) 
= 1 (F5) 

We find that 

(g)iL/^(5)|Q;,/3;/3,7;7,a) = \ga,gP;gf3,g'y;g-f,ga). 

The state jvI/pSRE) on a ring is invariant under the sym- 
metry transformation. So, defined in eqn. (26) is in- 
deed a solution of eqn. (12). We can obtain one solution 
for every cocycle in H^{G,U{1)) and each solution cor- 
respond to a SPT phase in 1 dimensions. 
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FIG. 30: (Color online) A segment of ID chain with open 
ends. The degrees of freedom not on the end form maximally 
entangled dimer states. 



3. States at the chain end form a projective 
representation 

Now let us consider the action of on-site symmetry 
transformation ®iU'^{g) on a segment with boundary (see 
Fig. 30): 

®i?7^(g)|ai,/3;/3,7;7,a2) (F6) 
= f2{ai,P,g,g*)f2{f3,i,g,g*)f2{'i,a2,g,g*)x 

19^1, 9P; 9/3, 91; 91, 90-2) 
^ '^2{ai,g^^g*,g*) 

V2{a2, g^^g*,g*)' 

or 

®iC/'(<?)l«i,«2)o = '''i"''^X'^*i lffai,ga2)o, (F7) 
V2{a2,g ^g*,g*) 

where 

Eqn. (F7), is the same as eqn. (26) and eqn. (27). Thus, 
®iU'^{g) form a linear representation of G. 

Note that |Q:i,a2)o is the ground state of our fixed- 
point model on a segment of chain, where all the internal 
degrees of freedom form the maximally entangled dimers 
(just like the ground state on a ring), while the boundary 
degrees of freedom are labeled by ai and 02 on the chain 
ends. «! and 012 label the effective low energy degrees of 
freedom |ai,a2)o- Those low energy degrees of freedom 
form a linear representation of the symmetry transforma- 
tion as expected. Eqn. (F7) describe how the boundary 
low energy degrees of freedom |ai,a2)o transform under 
the symmetry transformation. 

On the other hand, the symmetry transformation 
(^iU^{g) factorize (see eqn. (F7)), also as expected. This 
is because the degrees of freedom labeled by ai and a2 
are located far apart and decouple. We have (on the end 
whose states are labeled by ai) 

<»iU\g)\ai)o = '^2{ai,g'^g*,g*)\gai)o. (F9) 
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FIG. 31: (Color online) (a) The graphic representation 



of 



— ; , , . — . (b) The graphic repre- 

sentation of i^2{g'~^ g* , g* , g") which allows us to show 



'^2(ai.9' ^9*)'^2(ai.g '"9*. 9*) 

i'2(ai,9'^^g*,g*) 



= 1^2 ( 



g ,g g ,< 





(b) 



IfT* 



FIG. 32: (Color online) (a) The graphic representation of 
the phase factor eqn. (Gl). (b) The graphic representation 
of the phase factor eqn. (G2). The graphic representations 
indicate that the two phases are the same. 



Such transformation satisfies (see Fig. 31) 



V2 {gai , gg'~'^g* , g* )i^2 (ai , g~^g* , g* ) ^ 

'^2{ai,g'~^g*,g*) 
'^2{ai, g'~^g* , g^^g*)'^2{ai,g^^g* , g*) 

V2{ai,g'-^g*,g*) 

^2{g'-'g\g-'g\g*)®i u\g')\a,)^. 



$i f/'(g')|ai>o 
®i C/'(g')l«i)o 

(FIO) 



We see that the degrees of freedom on one end form a 
projective representation labeled by the 2-cocycle V2, the 
same 2-cocycle V2 that characterize the symmetry trans- 
formation of the SRE state. 



Appendix G: (2+l)D solutions of eqn. (12) 

1. U^{g) is a linear representation 

To show that defined in eqn. (30) is a lin- 
ear representation of G, let us compare the action of 
two symmetry transformations: U^{g)U^{g~^g') with 
the action of Lf^^g'), which changes |ai, 012, 0^3, 04) to 



a 




a 



FIG. 33: (Color online) A 2D |^pSre) state on a torus. In 
I^pSre), the linked dots carry the same index a,/3,7, ... 




g-'g' 



FIG. 34: (Color online) The graphic representation of the 
phase F3 in eqn. (G3). F3 is the value of a 3-cocycle z^a on the 
above complex with a branching structure. Note that the top 
pyramid and the bottom pyramid each form a solid torus (due 
to the periodic boundary condition) and the whole complex 
is a sphere. So F3 = 1. Note that the two pyramids on top 
and blow each small square represent the phase factor /s in 
eqn. (30). 



Ig'ai, g'a2, g'as, g'a4| >. One has a phase factor 

'^3{ai,(^2,g^^g*,g*)'^3{(^2,a3,9^^9*,9*) 

1^3(04, a3,g-'^g*,g*)i'3{ai,a4, g~^g*,g*) 
V3{gai,ga2, gg'~'^ g* , g*)v3{ga2, gas, gg'~^g*,g*) 

1^3(9^4, 90^3, gg'"\9*, g*)i'3(gai, ga4, gg'~^g*,g*) 

_ V3{ai,a2,g~^g* ,9*)i'3{oL2,a3,g-^g* ,g*) 



V3{a4,a3,g '^g*,g*)i^3{ai,a4,g ^9*, 9*) 



(Gl) 



'^3(Q^i,«2,g' g*,g g*)t^3(«2,a3,g' 9*, 9 9*) 

j/3 (^4 , aa , g'- ig* , g-ig* )i/3 (ai , a4 , g'" \g* , g-ig* ) 
and the other has a phase factor 

'^3{ai,a2,g'~^g*,g*)i^3{a2,a3, 9'~^g*,g*) 



z^3(a4,a3,g' ig*, g*)i^3(ai, 0:4, g' ^g*,g*) 



(G2) 



From their graphic representations Fig. 32, we see that 
the two phases are the same. Thus J7*(g) form an unitary 
representation of the symmetry group G. 



2. U'{g) satisfies eqn. (12) 

Following a similar approach as for the (l-l-l)D case, 
we can also show that the state I^'psre) on a 2D complex 
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FIG. 35: (Color online) A 2D I^I'pSre) state on an open 
square. In I^I/psre), the linked dots carry the same index 
Qfi, a2, /3, 7, ••• The indices on the boundary are given by 
ai,a2, ... The indices inside the square are given by ... 



where n(ij) ^ product over the nearest neighbor bonds 
{uli I* ~ j\ — Ij around the boundary. The direction 
« — > j is the direction of the bond and Sij — I f i > j , 
Sij = —1 I i < j. Since is independent of the indices 
/3,7, ... that are not on the boundary, we find 

(8)i t/'(5)|{aJ)o ^Y{i^3" {ai,aj,g~^g* ,g*)\{ga,})o 

(G6) 

where |{ai})o is the SPT state with a boundary which 
depends on the indices {ai} on the boundary: 

|{a*})o= Yl |ai,a2,/3,7,-)- (G7) 

/3,7,...eG 



FIG. 36: (Color online) The graphic representation of the 
phase F3{g,g*;ai,a2,l3,-y,...) in eqn. (G4). Compare to the 
complex in Fig. 34, the above complex do not have the peri- 
odic boundary condition. 

(see Fig. 33) that is a boundary of another graph is 
invariant under the symmetry transformation (8)iC/* (see 
Fig. 34): 

(^iU'l^pSRE) = FsI^pSRe) = I^pSRe) (G3) 

So, defined in eqn. (30) is indeed a sohition of 
eqn. (12). We can obtain one solution for every cocy- 
cle in ?^^(G', [/t(1)) and each solution correspond to a 
SPT phase in 2 dimensions. 

3. The action of ®U^{g) on I'I'psre) with boundary 

Now let us consider the action of ^iW^ on a state in 
Fig. 35 with a boundary (see Fig. 36): 

®i C^*(.9)|ai,a2,^,7, ••■) 
= ^3(.g,5*;ai,a2,^,7, ••■)l5"i,5a2,ff/3,57, •■•) (G4) 

From the Fig. 36 and the geometric meaning of the co- 
cycles, we find that 

■f3(5,3*;ai,a2,/3,7, •••) = H ^a'" ' 5~^5* , 3*) 

(G5) 



We see that the action of ®iU^{g) on |{ai})o is very 
similar to the action of a single U'^{g) on a single site 
(compare Figs. 12 and 36). Using a similar approach, 
we can show that (S)iU^{g) indeed form a linear represen- 
tation (see Fig. 32), when viewed as an operator Ub{g) 
acting on the boundary state |{ai})o. 

To summarize, we discussed the form of on-site symme- 
try transformations ®iU'^{g) in a basis where the many- 
body ground state is a simple product state. We find 
that different on-site symmetry transformations can be 
constructed from each 3-cocycle in 'H^[G, /7t(1)]- 

We would like to stress that, in such a simple basis, the 
symmetry transformation ®iU'^{g) on the boundary (G6) 
has a very unusual locality property: Due to the non- 
trivial phase factor W^ij) 1^3^' i<^i,C(ji 9~^9* , 9*), we can- 
not view Ub{g) (acting on the boundary state |{Q;i})o) as 
a direct product of local operators acting on each bound- 
ary sites \ai). (Note that we can view the boundary state 
|{Q;i})o as |{aj})o = (Sijsboundarylaj)-) Therefore, Ub{g) is 
not a on-site symmetry transformation on the boundary. 

In the above, we have viewed i as effective sites on 
the boundary with physical states ja^) on each site. We 
see that the symmetry transformation is not an on-site 
symmetry transformation. If we view, instead, each near- 
est neighbor bond (ij) as an effective site with physical 
states \aiaj) on each site, then the symmetry transfor- 
mation will be an "on-site" symmetry transformation, 
but the states on different bounds are not independent 
and |{Q!i})o 7^ «) 

(i^) (z boundary I ^i^j) • 

Thus in a basis where the many-body ground state 
is a simple product state, although ®iU'^{g) is an on- 
site symmetry transformation when acting on the bulk 
state, it cannot be an on-site symmetry transformation 
when viewed as a symmetry transformation acting on the 
effective low energy degrees of freedom on the boundary 
when the 3-cocycle is non-trivial. This is the non- 
trivial physical properties that characterize a non-trivial 
SPT phase in (2+l)D (see appendix I and Ref. 59 for 
more details). 
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FIG. 37: (Color online) The graphic representation of the 
product of the phase factor Yl{ij) ^3^^ i'^^^^j ^ 9* ' 9*) 
eqn. (G6). 



Appendix H: Two symmetry representations in our 
fixed-point model 

In the old basis in the path integral formalism [see 
eqn. (54)], the wave function is complicated, but the 
many-body on-site symmetry transformation has the fol- 
lowing locality structure 



(HI) 



where U'^{g) is the symmetry transformation on the i*'' 
site 



U'igm = \99i 



(H2) 



This is the definition of the so called on-site symmetry 
transformation. 

In the new basis[see eqn. (56)], the wave function is 
simple but the many-body symmetry transformation is 
no longer an on-site symmetry transformation. It has 
the following form 



\{99i}M)' 



®^ U\9)\{9i\M)' 
Y\{^J...*}^l+d'i9u93,■■■,9*) 

Yl{tj...*}^i+d'{99i.99j.-.9*) 
X{{ij...*}Vi+d'i9^^9j^-^9*) . ., 

l^{^j...*}^l+d {9^,9j,■■■,9 ^9*) 

\{99^}My (H3) 



where W" (g) is the symmetry transformation that acts on 
the i*'' and {i ± 1)*'' sites [in (1+1)D for example] 

U''i9)\9i-i,9i^9i+iy = h{9i-i, 9i^ 9i+ii 9)\9t-ii 99ii 9i+i)' ■ 

(H4) 

Here the phase factor /2 is given by the 2-cocycles 



f2{9i-i,9i,9i+i,9) 

1^2(91-1, 9i, 9*)'^2{9i, 9i+i, 9*) 
^2 (ffj- i,99i.9*)v2 {99t , 9t+i , 9* ) 



(H5) 



where 9* is an fixed element in G. For example we may 
choose 9* = 1. 



FIG. 38: (Color online) The graphic representation of the 
product of the phase factor — ^{ij.,...}^ — ( ^. j. — ^£__) — 

n{ij....} "i+d " (Q.,aj,...,9 19*) 

eqn. (H6) for a 2D complex (ai,...,a7) with a boundary 
(ai, ae). 



We have seen that, in the new basis, we still have 
'^iU''{9)\{9i}My = |{5sJm)' if the state |{gjA/)' is de- 
fined on a complex M which is the boundary of another 
complex Mcxt- It is hard to see the non on-site struc- 
ture of (E)iU^{9)- To expose the non on-site structure 
of <E)iU^(9) in the new basis, let us consider the action 
of <E)iU'^{g) on a state defined on a complex that has a 
boundary. In this case, we still have 



^^ C/Xff)|{ft}Af)' 

U{^J...*}l'l+d'i9^,9j,■■■,9* 



n 



{ij- 



'{9i^9j^ ■ 



,9-'9* 



-\{99^} 



(H6) 



But now, the phase factor is not equal to 1. 

In Fig. 38, we give a graphic representation of the 
above phase factor - — (">."j.---^g ) ^ 

\A{ij...,)''i+d' {ai,aj,.--,g Iff*) 

complex with a boundary. We see that the complex in 
Fig. 38 and Fig. 37 have the same surface. So the phase 
factor represented by Fig. 38 equal to that represented 
by Fig. 37. So eqn. (G6) is the same as eqn. (H6). 

We have discussed two ways to classify SPT phases. 
The first way to classify SPT phases is to classify sym- 
metry transformations that act on simple wave func- 
tion I^'psre), which lead to eqn. (G6). The second 
way to classify SPT phases is to classify fixed-point 
action-amplitude (the topological terms) which lead to 
eqn. (H6). The above analysis indicates that the two 
ways to classify SPT phases are equivalent. 

The equivalence between the two formalisms eqn. (G6) 
and eqn. (H6) will become more clear after a duality 
transformation. 

In the following we will show that the ground state 
wave function (56) in the Lagrangian formalism is dual 
to the ground state wave function (25) in the Hamilto- 
nian formalism discussed in Sec. IV and V. Furthermore, 
after the duality transformation, the the symmetry rep- 
resentations (55) are the same as that defined in eqn. 
(26) or eqn. (G6). 

Let us illustrate the above result in ID. Firstly, we 
introduce the dual transformation which maps a state to 
its dual wave function living on the dual lattice. In the 
dual transformation, the bases [(/;) at site i correspond to 
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g\ g\ g\+lg'i+l g'i+2g''i+2 g'.+3 g\+3 

i i+1 i+2 i+3 

gi gi+1 gi+2 

i i i 



FIG. 39: (Color online) The dual transformation in the new 
bases in ID. 



the bond \gl, g\^i) in the dual lattice (see Fig. 39), where 
gl — g\+i — di and the amplitude of the configuration 
remains unchanged. In this way, we obtain the dual wave 
function -^diiglgl}) of vl/({gj). 

Now we introduce the new bases K^', through the 
LU transformation introduced in eqn. (56), 

M.anr = \{^2{g^,g^+ug*)M,gn) 

i 

= Wv2{gl,gl+i,g*)\{g\,gl}) 

i 

= n h(5Ui>5r+i,5*)l5^+i,5[+i)] (H7) 

i 

In the new bases, the fixed point state in the dual lat- 
tice becomes a direct product of bonds. Notice that the 
previous local unitary transformation in eqn. (56) be- 
comes on-site unitary transformatioir. Furthermore, in 
the new bases the symmetry representation also becomes 
on-site and is fractionalized into two 'projective' opera- 
tions: 

®rU\g)\{glg:})' 
= W_^2 (ff'+i ' .91+1 ' 9* ) |{.9.9!+i , gal+i}) 



f: — 7f — ; 

6i+y tei+x+y 



n 



. i^2\gi+i, gi+i, g g ) 

n'^2{gi+i,g^^g*,g*) .r i ^ , . 

—r^ _i , , J m+i,. 9.9,;+ il H8 
. '^2[gl+i,g ^g*,g*) 

Above formula is the same as eqn. (27). 

Similarly, we can illustrate the result in 2D. Now the 
basis \g^) correspond to \gl , gf+x, gf+^+y, gf+y) in the dual 
lattice (see Fig. 40). After the dual transformation, the 
wave function ^'({f/i}) becomes "if d{{gj , gf , gf , gf}) (here 
gl = gf+a: = 9f+x+y = gf+y = gi)- Again, we introduce 
the LU transformation 

Mglglglgh)' 

n'^3{gi,gi+x,gi+y,g*) ,r i 2 3 4^\ 
—^\{9i : 9z : gt ) 9i } 
^ >^3{9^+x,g^+y,gi,g*) 

n '^3{gj,gf,9j,g*) , ^ ^ , „3 ^ttqn 
, >^3{gf,g{,9i,9*) 



i+y i+x+y 

I ^ &']fx+y I ' 

3 i „4 i+x 



1^ f: — ^ 
i i i 

^ — >l i — 4 — i 



FIG. 40; (Color online) The duality transformation in 2- 
dimension. The green dots represent the dual lattice of the 
red dots. In the new bases, the wave function in the green 
lattice is the same as the one introduced in Sec. Ill and IV. 



where i — i + x + y. The LU transformation between 
the old bases and the new ones is an on-site. In the new 
bases, the fixed point wave function is a direct product 
of plaquettes. The symmetry operation now becomes 

®^U\g)\{glglglgt}y 

nMgl,gf,gi,9*) i 2 3 4i\ 
, 4 2 1 .^ {ggi^ggi^ggh 99^}) 
, M9j,gj,gi,9 ) 

^ T-r ^3(gf,g^gf,g*) -pr M9f,9lgj,g-'g*) 
^}''^3{9f,gf,gl9*)^}'M9l9f,glg-'g*) 

>^\{ggj,9gi,ggl,ggt}y 
_ j^ M9h9t,g~^g*,g*>3i9t^9h9^^9*,g*) 
M9l,9i,g~^g*,g*)M9i,9j,g~^g*,g*) 

^\{99l99!,g9l99t}y (HIO) 

Above equation agree with eqn. (31). 

From the above examples, we can see that after the 
'dual transformation' the ground state wave function and 
its symmetry representation in the Lagrangian formal- 
ism are the same as the Hamiltonian formalism as we 
discussed in Sec. VI. 



Appendix I: (2 + 1)D SPT states constructed from 
3-cocycles and matrix product unitary operator 

Based on the 3-cocycles vsigo, gi, g2^ 93) of group G, 
we can construct short range entangled models with 
SPT order as discussed in section VI C (also discussed 
in appendix G). In order to assess the non trivialness 
of the SPT order of a certain model, in Ref. 59 we de- 
veloped the tool of matrix product unitary operators 
(MPUO) and used it to show that the particular model 
we gave in that paper-the CZX model-has very special 
boundary properties and hence nontrivial SPT order. In 
this section, we are going to apply the MPUO method 
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to the general models constructed in section VIC and 
show that for the model constructed from a 3-cocycle 
^'3(S'0)ffi)52>53), the effective MPUO on the boundary 
transform with the same 3-cocycle. Therefore, according 
to the result in Ref. 59, models constructed from non- 
trivial vsigo, gi, g2, gs) must either break the symmetry 
or have gapless excitations if the system has a bound- 
ary. Moreover, we can show the contrary for models con- 
structed from trivial lysigo, ffi, <72, 53)- That is, for models 
constructed from trivial ^'3 (go j <7i j 52 , 53 ) we are going to 
explicitly construct a short range entangled symmetric 
state for the effective symmetry on the boundary. For 
basic definition and properties of MPUO, see Ref. 59. 

We consider in this paper models with on-site symme- 
try of group G. SPT order exist in models whose ground 
states on a closed manifold are short range entangled and 
symmetric under the on-site symmetry. The ground state 
is unique and gapped. If the system has boundary, on the 
other hand, there are low energy effective degrees of free- 
dom are on the boundary. The effective symmetry on the 
boundary however, may no longer take an on-site form. 
In general, the effective symmetry on the ID boundary of 
a 2D model can be written as a matrix product unitary 
operator 

[/= J2 Tr(ri'*'ir='*^=...r"'^'")|i'ii'2...i^)(iii2-.iw| 

(11) 

where i and i' are input and output physical indices and 
for fixed i and i', T*'' is a matrix. 

For the models defined in section VIC, the effective 
symmetry U{g) on the boundary takes the form(see ap- 
pendix G) 

U{g)\{a^})=l[4'^{a,,a„g-'g*,g*)\{ga,}) (12) 

where Yiij is a product over the nearest neighbor bonds 
I* ~ j| — 1; around the boundary. The direction 
i — J is the direction of the bond and = 1 if z > j, 
Sij = —1 if j < j. This symmetry operator on a ID chain 
can be expressed as a MPUO. If the bond goes from 
to ai+i 

T""^"'{g) = X^a.+i '^3'^("»'"*+i'5"^3*'3*)l"»)("i+il 

(13) 



This can be reduced to 
1^,(31.92)""^^^^"' 



^ i/g ^(ai,ai+i,52 ^5*,5*)|a*)(ai+i| (K) 



other terms are zero 
If the bond goes from a^+i to a^, 

^^"^"'(3) = Ea.+i '^3(^1+1, at, g-'^ g* , g*)\a,) {a^+i\ ,Va, 
other terms are zero 

(14) 

Now we compose multiple MPUOs and find their reduc- 
tion rule. We will see that the reduction rule is related to 
the same v^. First, the combination of Ti[g2) and Ti(gi) 
gives (if the bond goes from to a^+i) 

Ti(5i,ff2)"-^^''="' = z^3-i(a„a,+i,52-V,5*) 
v*'''^^^\g2a, , a -+1 , g^^g* ,g*)\a„ g2a^) {a,+i , a 1 

(15) 



by applying the following projection operator to the right 
side of the matrices 

P'' (17) 

31,92 y 
= X! '^3^i^^+^^92^9l^9*,g2^9*^9*)\0i^+l,92a^+l){a^+l\ 

and the hermitian conjugate of 



91:92 

^^'^3^(oi„g:^'^g^^g*,g2^g*,g*)\a,,g2a,){a,\ (18) 

to the left side of the matrices. This is because, 

1^3 ig2ai , 52^1+1 , .9^^9* , 5* ) 
= 4^''\a,,a,+i,g2'g^'g*,g2'g*) (19) 
and the 3-cocycle condition of i^^ 

iy3{a,, ai+i, g^^ g* , g*)iy3{ai, a,+i, g2^ gi^ g* , g2^ g*) 

g2^9T^g* , 92^ g* , g*)^^ ("i+i , g2^gi^g* ,g2^9*,9* 
= i'3{ai,a^+i,g2^gi'^g*,g*) 

(110) 

It is easy to check that the same reduction procedure 
applies when the bond goes from oii+i to a^. The above 
definition of and P^ has picked a particular gauge 
choice of phase for P' and P^ . 

Next we consider the combination of three MPUOs and 
find the corresponding 3-cocycle associated with different 
ways of combining the three MPUOs into one. If we 
combine T{g2),T{gi) first and then combine T{gig2) with 
T{g3), the combined operation of Pg^.ga and Pgig2,s3 
(we omit the site label i) 

{Pgua.®I)Pgi92,9. (I^) 

^V3,{a,g^^g2^g^^g* ,g^'^g2^g\g^^g*)-x. 

a 

i^3(a, g3^g2^9i^g* , g^^g* , 5*) !«, 53a, 5253a) ("I- 

On the other hand, if we combine T{g3),T{g2) first and 
then combine 7(17253) with T{gi), the combined operator 
^92, 93 ^91,9293 is 

{I ® Pg2,g3)Pgi,g293 ~ 

J2a ^3(a,5r^5^^5*,53^^5*,5*)x 

i^3(a, g3^g2^9i^g* , 93^g2^g* , g*)\a, gsa, 5253a) («! 

(112) 

These two differ by a phase factor 

^^3(53 52 5i 5 ,53 52 5 ,53 9,9) (US) 
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Hence we see that, the reduction procedure of T's is as- 
sociative up to phase. The phase factor is the same 3- 
cocycle that we used to construct the model. From the 
result in Ref. 59 we know that if is nontrivial, the 
model we constructed has a nontrivial boundary which 
cannot have a gapped symmetric ground state. It must 
either break the symmetry or be gapless. Therefore, the 
model constructed with nontrivial 3-cocycles belong to 
nontrivial SPT phases. 

On the other hand, if the model is constructed from a 
trivial 3-cocycle, the boundary effective symmetry does 
allow SRE symmetric state. Actually, the SRE symmet- 
ric state on the boundary can be constructed explicitly 
for the models discussed here. If is trivial, it takes the 
form of a 3-coboundary 



J^3(ff0,3l,ff2,53) 



M2(gl,g2,ff3)M2(gO,gi,ff3) 

M2(5o, 52, 53)^2(50, 51, 52) 



(114) 



where ^2 is an arbitrary 2-cochain. Note that it is not 
necessarily a cocycle. The effective symmetry on the 
boundary can hence be written as 



1. Canonical choice of cocycles 

Let us consider an n-cocycles v'^ which satisfies the 
condition Av'^ — 1. By a proper transformation by a 
coboundary 6„: = v'^Jd"^ , we can choose a particular 
cocycle Vn in a given cohomology class that satisfies 

J^n(gO,gO,5l, •■•,5n-2,5n--l) = 1 (Jla) 
'^n(5l,g07gO,.--,.gn-2,5n-l) = 1 (Jib) 

J^n(5l,52,53---,5n-l,gO,go) = 1- (Jlc) 

To this, let us focus on eqn. (Jla). To prove that the 
choice eqn. (Jla) is valid in general, it is equivalent to 
prove that we can always choose a cocycle in a cohomol- 
ogy class that satisfies 



^'n( gO, gO^ , 5l, 9n+l-m) — 1, 
m terms 



(J2) 



U{9)\{0^:}) = 



i-J \M2(Q!,g ^g' ,g*)fj-2{ai 



(115) 



The fi2i<^i, g~^g* , g*) terms cancel out in the prod- 
uct of phase factors, and the remaining terms can 
be grouped into two sets Yiij 1^2^ i'^it'^j^ g*) 



M2 ' 



js(g) 



{gai,gaj,g* 



Define e(,g) = H.^ Ea.,a, M2'' (a^, , 5*)|aiaj) (a^aj I- 
&{g) is a product of local unitaries. It is easy to see 
that 



where m is the number of repeating index go with 2 < 
m < n + 1 and gi ^ go- 

Firstly, we will show that a cocycle can satisfy eqn. (J2) 
for m = rt + 1, which means 



^n{go,go, ■••,5o) = 1- 



(J3) 



n+l 



If n is odd, eqn. (J3) can be easily shown from the cocycle 
condition 



U{g)^e\g) 0(5) 



(116) 



(a complex conjugation operation needs to be added if 
U{g) is anti-unitary). The term in the middle is an on- 
site operation which permutes the basis. It has a simple 
symmetric state which is a product state ^iC^^. \c(i))- 
Therefore Q^ig) (8)i {J2a- ^ symmetric state of 

U{g). Because (^i{J2a- ^ product state and &^{g) 

is a product of local unitaries, this is a short range en- 
tangled state. Therefore, we have explicitly constructed 
a short range entangled symmetric state on the boundary 
if the model is constructed from a trivial 3-cocycle z/3. 



Appendix J: Calculations of group cohomology 

In the section, we will calculate group cohomology for 
some simple groups. We will first present a direct calcu- 
lation from the definition of the group cohomology. Then 
we will present some more advance results. 



(dt^n)(5o,.9o,5o , ■■■,5o ) = t^n ( 5o , 50y ■ , 5o ) = 1- 

n+2 n+l 

If n is even, then we can introduce a cobound- 
ary bn{ go,go , ^...,go) = {dnn-i){ go,goy,go ) = 

n+l n+l 

/^n-i(507 50)7 where ^„_i is a cochain. If we re- 

n 

quire that M«-i( 5o, 5o ) = '^n ( 5o , 5oy ■ , 5q ) , then 

n n+l 

after the gauge transformation Vn = v'^i>~^ we have 
i'„((7oj 5o, •••7 5o) — 1- Thus, we have proved the validity 

n+l 

of eqn. (J2) in the case m = n + \. 

Now we show that a cocycle can satisfy eqn. (J2) for the 
case m = 2fc + l (1 < fc < ["./2], where [n/2] is the integer 
part of n/2, i.e., [n/2] = n/2 if n is even and [n/2] = (n — 
l)/2 if n is odd), namely, Vn { go , ■ ■ ■ , 50 , 51 , 52 , . . . , 5«-2fc ) = 

2fe+l 

1. Here we assume gi 7^ go. Again, we introduce the 
gauge transformation j/„ = v'^bn^ with 6„ = (d/x„_i). 
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We requires that the cochain satisfies 



2 k 

Ki9o,9o, •■•,50,51,52, ■••,5n-2fc) 



2fe+l 



A^n-l(50,50, •■■,50,52,53, ■•■,5n-2fe) 

' V ' 

2fc+l 

/"ri-l(50,50, •••,50,51,53, •••,5n-2fc)^ 

2fc+l 



/^„_i l5o,5o, •••,5o,5i,52, •••,5n-2fe-iJ 

V ' 

2fe+l 



(J4) 



for fc = [n/2], fc = [7i/2] - 1, fc = 1 in 

sequence. Notice that there is one-to-one corre- 
spondence between ^'n( 5o,5o , •■•,5o ,5i,52, ■••,5n-2fc) and 

2fc+l 

^„_i(go, ■■■,50,51,52, •••,5n-2fc)- In the square bracket, 

2fc 

the terms which have odd number of successive index 
go are free variables. Equation (J4) can always be sat- 
isfied by letting Mn-i(5o, •••,5o,5i,52, •••,5n-2fc) equal to 

2 k 

the right-hand side of eqn. (J4) (we will illustrate it by 
several examples) . In other words, eqn. ( J4) is a constrain 
for the components of /k„_i which have even number of 
successive index go, and at the same time the components 
of fin-i which have odd number of successive index go can 
be chosen arbitrarily. Thus, we can always find a 
that satisfies eqn. (J4). 

From eqn. (.14), we have 



CO cycle condition 

{dUn){go,go, ...,50,51,52, ...,5n+l-2fc) 

V ' 

2k + l 

= Vnj gOi •-, 50 ,51,52, •■•,5ri+l-2fc) 
2k 

'^n ^( 50,50y,50 , 52, 53^-^, 5«+l-2fc)-^- 

2fe+l 

^ (50,50, •••,50,51,52, •••,5"-2fe) 

V ■' 

= 1, 

we obtain i^„(5o, •••, 5o, 5i, 52, •••, 5n+i-2fc) = 1, which 

2k 

proves eqn. (J2) in the case m — 2k. 

Above proof includes the cases of 1 < k < [n/2]. 
Together with eqn. (.13), we have finished the proof of 
eqn. (J2), or equivalently, eqn. (Jla). Notice that the 
only requirement in the proof is eqn. (.14). 

To prove eqn. (Jib) in general, it is equivalent to prove 
the following equations 

'^«(5l,gO, •••,§□, 52, 53, •••,5n-2,5«+l-m) = 1, (J5) 



here m is the number of repeating index go with 2 < m < 
n and gi ^ go (the case g\ = go reduces to eqn. (J2) and 
has been proved already). 

Let begin with the case m — n, namely. 



'^n(5i,5o,5o, •••,5o) = !• 



is odd, we introduce — 



(J6) 



If 

^n(5i,5o,5o, •••,5o) 



<6„\ with 



(d/i„_i)(gi,.go,5o, •••,5o) = 



1 ( 50 , 50y • , 5q ) Mn-i^ 1 (51 , 50 , • - , 5o ) • If we require 

n n—1 

A*"-i(5i,5o,5o, •••,50) 



fcn( 50,50y,50 , 51,52, •••,5n-2fc) 
2fc+l 

(d^„_i)(go,.go, •••,.go,5i,52, •••,5«-2fe) 

v ' 

2/c+l 

i^n(5o,5o, •••,50., 51,52, .■.,g„-2k)- 

2k+l 



Consequently, we obtain 

'^n{go,go--,go,9i,92,--,9n-2k) = 1 after the gauge 

V ' 

2fc+l 

transformation Vn = ^rfin^- From this result and the 



= t'n (51 , 50 , • - , go) M«- 1 ( 50 , 50y • , 50 ) , 
n — l n 

then we obtain fe„(gi, go, 5o, •••, 5o) = t^,'i(5i,5o,5o, •••,5o) 
■> „ ■' V ■' 

n n 

and consequently m„((7i, go, 5o, •••, 5o) = !• If ft is even, 

n 

then from the cocycle condition 

(dj^n)(5i,5o,5o,5o, •••,5o) 



n+l 

i'n(5o,.9o,5o, •••, 50)^^,7^(51,50,50, 

1, 



,5oJ 



n+l 
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we obtain Vn{gi, cjq, go, go) = 1, here we have used 

n 

eqn. (J3). 

Now we prove eqn. (J7) for the case m = 2fc + 1 (1 < 
k < [{n — l)/2], where [{n — l)/2] is the integer part of 
{n- l)/2), namely, z^n(gi, go, ■^,go^ , g2,g3, ■■■,gn-2k) = 1 

2k+l 

with gi,g2 7^ go- Again, we introduce the gauge trans- 
formation i/n = with bn = (d/i„_i). We requires 
that the cochain satisfies 

M«-i(gi. go. - .go i 32,93, ...,g„-2fc) 

2fc 

= ^'n(gi, gO,gOy,gO , 92,93, ■■■,gn-2fc) 
2fc+l 

Mn-i(go,go, ■■•,go,g2,g3, ■■■,gn-2k) 
^ ' 

2fc+l 

Mn-i(gi,go,go, ■•■,go,g3,g4, ■■■,gn-2k)--- 

V ' 

2fc+l 

mI-V (gi,go,go,---,go.,g2,g3,---,gn-2fe-i) (J,7) 

2/c+l 



for fc = [(n - l)/2], fc = [(n - l)/2] - 1, 
fc = 1 in sequence. Again, there is one-to-one cor- 
respondence between iy'^{gi, go, go, go, g2, gs, gn-2k) 

V ' 

2/c+l 

and A*«-i(gi, go, ■^■,go , 92, g3, ■■■,g«-2fc)- Equation (J7) 

2fc 

can always be satisfied by constraining the value of 
f^n-iigi^ goj_-j_go, g2, gs, gn~2k) to equal to the right- 

2 k 

hand side (we will illustrate it by several examples). 
From eqn. (J7), we have 



^w(gi, go,goy,go , g2,g3, ■■■,gn~2k) 

2k+l 

= (dAi,i-i)(9i,9o,9o,---,go,g2,g3, ■•■,gn-2fc) 

" V ' 

2k+l 

= ^'n(gi,gO,gO, ■••,g0,92,93, ...,9n-2fc)- 

V ' 

2fe + l 



Consequently, after the gauge transformation Vn — 
^n^n^ we obtain t'„(9i, 9o, 9o, 9o, g2, gs, gn-2fc) = 
1. From this result and eqn. (J2) and the cocycle condi- 



tion 



(dt^n)(9i, go,9o , ■■■,go , 92, gs, ■•■,gn+i-2fe) 

2k+l 

'^n( go,gOy,gq ^, g2,93, ••■,9n-|-l-2A:) 
2fc+l 

i^rT^ (91, 90, - ,.90 , 92,g3---,9n+l-2fc) 

2fc 

'^n(gi, go, goy ■, gg ^, 93, 94, •■•,.9n+i-2fc)--- 

2fc + l 

^ (gi,go,go, ■•■,go,g2,g3, ■••,gn-2fc) 



2fc+l 



1, 



we obtain Vn{gi, go-, ■^■,9q , 92,93, ■.■,9«-n-2fc) = 1, which 

2fe 

proves the case m = 2k. 

Above proof includes the cases of 1 < fc < [(n — l)/2]. 
Together with eqn. (J6), we have finished the proof of 
eqn. (J5), or equivalently, eqn. (Jib). Notice that in 
the proof we have used two conditions eqn. (J4) and 
eqn. (J7). Obviously, they can be satisfied simultane- 
ously. 

The remaining part of eqn. (Jl) can be proved by the 
same procedure and will not be repeated here. We stress 
that all of the equations in eqn. (Jl) can be satisfied 
simultaneously, because in proving different equations we 
are fixing the values of different classes of components of 
the {n — l)-cochain 

As examples, let us illustrate that eqn. (J4) and 
eqn. (J7) can be satisfied simultaneously for n < 4. When 
n = 2 (k = 1), eqn. (J4) becomes 

Mi(go,go) = j^2(go,go,go), 

which can be satisfied obviously. We do not need to con- 
sider eqn. (J7) for n — 2. 

When n = 3 (fc = 1), eqn. (J4) becomes 

Ai2(go,go,gi) = t'3(go,go,go,gi)M2'^(go,go,go), 

which can be satisfied by constraining the value of 
/j,2(go, go, gi) to be equal to the right-hand side. On the 
other hand, eqn. (J7) becomes 

A'2"^(gi,go,go) = v'sigi, ga, go, gQ)^J■2^ (go, go, go), 

which can also be satisfied by properly choosing the value 

of M2(gi,go,go)- 

Finally, when n ^ 4, there arc two cases k — 2 and 
k — 1. For k = 2, eqn. (J4) becomes 

Ai3(go,go,go,go) = v'^{go,go,go,go,go), 

which can be satisfied obviously. For k ~ \, eqn. (J4) 
becomes 

m(go,go,gi,g2) = ^'4(.9o,go,go,gi,92) 

x[M3'^(9o,go,go,g2)M3(go,go,go,gi)]" 
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which can be satisfied by constraining the value of 
^3(507 5o, 5i: ff2 to be equal to the right-hand side. 

One the other hand, when n — A and k — 1, eqn. (J7) 
becomes (We do not need to consider eqn. (J7) for k = 2) 

Ai3'^(5i'5o,3o,ff2) = ^^4(91750,50,50,52) 

[M3(5o,5o,5o,52)M3'^(5o,5o,5o,5i)]~ 

which can be satisfied by restraining the value of 

^^3(51, 50,50,52)- 

Now let us see what happens for the term 

^'n(go,5i,52, •■•,5n-i,go) with gi,gn-i 7^ go- Consid- 
ering the coboundary 

^1(50, 51, 52, ---,5n-i,5o) 
= (d^„_i)(go,5i,52, ---,5n-i,5o) 
= Ai«-i(5i,52, ---,5n-i,5o)Ai,T-i(5o,52, ---,5«-i,5o) 

■■■fJ-n-l I50,5l, ---,5n-lj, 

Notice that the two cochains /i„_ 1 (gi , 52 , - - - , 5n- 1 , 5o ) and 
mI-V (50,51, ---,5n-i) may cancel each other in some 
condition. In that case, 6ri((7o,5i,52, ---,5n-i,5o) has less 
degrees of freedom than v'^igo, 91, 92, ■-, 9n-i, 9o), so we 
CAN NOT always set t'„(go, 5i, 52, 5n-i, 5o) = 1- 

2. Group cohomology of Z2 

Using the properties obtained above, we will show that 
for the group Z2 — {E, a} (where E is the identity and 

'y'-E), 

H2™-1[Z2,[/(1)] =Z2, 

H^™[Z2, [/(!)]= Zi, m>l (J8) 

and for the time reversal Zj = {E, T} group, 

n^'^lZl, Ut{1)] = ^2, m> 1. (J9) 

We note that Z2 and Z2 are the same group. However, 
the generator in has a non-trivial action on the mod- 
ule. Also C/(l) and Ut{1) are the same as Abelian group. 
The subscript T in the module Ut{1) is used to indicate 
that the group Zj has a non-trivial action on the module. 

Let us begin with eqn. (J8). Firstly, 
i^2m-i(5o,5i, ---,52m+i) havc even number of 
group indices. From eqn. (-Jl), we can set 
i'2m-i(5o,5i, ---,52m+i) = 1 if any two neighboring 
indices are the same. So the only possible nontrivial one 
is when the group indices vary alternatively, namely, the 
component iy2m^i{E,a, E,a) = iy2m-i{(^, E, a, E). 
Considering the cocycle condition 

{dv27n-i){E,a,...,a, E) 
= V2m-i{(^,E, ...,a,E)i^^,l_^{E,E,a, ...,a,E) 

V2m-i{E,a, a,E, ...,cr, E)...V2m-i{E,(T, ...,E,a) 
= [iy27n-ii(7,E,...,a,E)]^ 

= 1. (JIO) 



here we have used eqn. (Jl). So we have 

l'2m-l{E,a, ...,E,(t) = V2ra-l{(y,E, ...,(J,E) = ±1. 

Now we need to show that these two solutions are not 
gauge equivalent. Consider the coboundary 

b2m~l{(7,E, ...,(T, E) 

= (d/i2m-2)(o', -E, ...,a-, £:) 

= fJ.2m-2{E, (J, CT, E)fi2^_2{a, a, E, CT, E)... 

fJ.2m-2{cr, E, CT, E, E)n2m-2i'^, ■■; E, (t) 
= M2ri-2('^,^,-^, ...,a,E)fi2m-2{(^,E,E, ...,a,E)... 

...IJi2m^2{o,E,...,a,E,E). (Jll) 

Notice that ^2m-2iE,<J, ...,a,E) is canceled by 
^2m-2^'^,^, ■■■,^,'^)- I'^ all of the remaining com- 
ponents, a pair of neighboring group indices are the 
same. The values of these components have been fixed 
in the gauge choice eqn. (Jl). Consequently, the value 
of the coboundary b2m~i{o', E, ...^ a, E) is also fixed. 
But there are two cocycles satisfying eqn. (JIO), so they 
must belong to two different classes. 

Secondly, V2mi9o, 9i, 92m) contains odd number of 
group indices. The only possible nontrivial one is 
i'2miE,a, ...,a,E) = i^2m.{o-,E, ...,E,a). Considering the 
coboundary 

b2m{E,a, ...,a,E) 
= {dn2„i~i){E,a,...,a,E) 
= iJ,2m-i{cr, E, cr, E)^i2^_^{E, E, a, ct, E) 

Ai2m-l(-B, CT, cr, E, cr, E) ...fl2,n-l{E , CT, ...,E, cr) 
= [/i2m-l(CT,S,-,CT,S)]'- (J12) 

Notice that the component fi2m-i{'^, E, ...,a, E) is free 
since it is not fixed by the gauge choice eqn. (Jl). 
So the b2m{E,a, ...,(T,E) has the same degrees of free- 
dom as i>2miE,(7, ...,a,E), and we can always set 
t^2m{E, a, a, E) — 1 with the gauge transformation 
i.„ = I/; 6-1 . Consequently, we have n^"'[Z2,U{l)] = 
Zi- 

Conditions are on the contrary for the time reversal 
group Z2 , because of the relation VniTgo, Tgi, ...,Tgn) = 
Vn^{ga, gi, ■■■,gn)- Corresponding to eqn. (JIO), we have 

{Av2m){E,T,...,E,T) 
= P2m{T, E, E, T)v2^{E, E, T, E, T) 
U2,n{E,T,T,E,...,E,T)...i,2,l{E,T,...,T,E) 

= [v2m(T,E,...,E,T)f 

= 1, (J13) 

which result in 

V2m{T,E,...,E,T) = V2,l{E,T,...,T,E) = ±1. 

Similar to eqn. (JIO), these two solutions are not gauge 
equivalent. Consequently, [/^(l)] = Z2. Simi- 
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larly, corresponding to eqn. (J12), we have 

b2m-l{E,T,...,E,T) 
= {d^l2,n-2){E,T,...,E,T) 

(T, E, E, T)^^^^_^{E, E, T, ...,E, T) 
^,2m-2iE, T, T, E, ...,E, T)...Ai2-„\_2(£;, T, T, E) 

= [^l2m^2{T,E,...,E,T)f... (J14) 

The free component ii2m-2{T, E, E,T) guaran- 
tees that the b2m-i{E,T, E,T) has the same de- 
grees of freedom as V2„^_i{E,T, ...,E,T), so we 
can set i^2m-i{E,T, E,T) = 1 and consequently 
H2'»-1[Z2^,[/t(1)] = Zi. 



3. Group cohomology of Zn over a generic 
Zn-module 

The cohomology group 'H'^[Zn, M] has a very simple 
form. To describe the simple form in a more general 
setting, let us define Tate cohomology groups Tk'^[G, M]. 

For d to be or —1, we have 

'H°[G,M] = M^/lmg{NG,M), 
n-^[G, M] = Ker(iVG, M)/IgM. (J15) 

Here , Img^Na, M), Ker(7VG, M), and IgM are sub- 
module of M. is the maximal submodule that is 
invariant under the group action. Let us define a map 
Ng ■■ M ^ M as 



a — > g • a, a G M. 

geG 



(J16) 



Img(A^G': M) is the image of the map and Ker(iVG, M) is 
the kernel of the map. The submodule IqM is given by 

IgM = { n (-9 • a)"' I 5] rig = 0, a e M} (J17) 

gea geG 

In other words, IgM is generated by {g ■ a)a^^, V g G 
G, a e M. 

For d other then and —1, Tate cohomology groups 
'H'^[G,M] is given by 



n'^ [G, M] =n'^[G, M] , for d>0 
■H'^[G,M]=H-d-i[G,M], ford<-l. 



(J18) 



For cyclic group Z„, its (Tate) group cohomology over 
a generic Z„-module M is given by^'^'^°° 



n''[Zn,M] 

where 



?i''[Z„,Af] if<i = 0niod2, 
?i-i[Z„,A'f] if<i=lmod2. 



(J19) 



n°[Z^,M] = Af^"/Img(iVz„,M), 
M] = Ker(iVz„ , M)/Iz^M. 



For example, when the group action is trivial, we have 
M^" = M and Iz„M = Zi. The map Nz„ becomes 
Nz,^ : a ^ a'\ For M = Z, we have Img(iV2„,Z) = 
nZ'and Ker(7Vz„,Z) = Zi. For M = U{1), we have 
Img(Afz„, Z) = C/(l) and Ker(iVz„, Z) = Z„. So we have 

if = 0, 

-H''[Z„, Z] = <( Z„ if = mod 2, d > (J21) 
[i if d = 1 mod 2. 



and 



H'^[^„,C/(1)] 



C/(l) ifd^O, 

!i if d = mod 2, d > (J22) 
-„ if d = 1 mod 2. 



which reproduces the result mentioned in Ref. 92, and 
the result obtained in the last subsection for d = 2. 

What does a non-trivial cocycle in H'^[Zn, U{1)] looks 
like? Since 'H^[Zn,U{l)] describes the ID unitary rep- 
resentation of Zn — {0, 1, fc, n — 1}, we find that 
the m'^ 1-cocycles in H^[Z„, [/(I)] = Z„ are represented 
by complex function u}["'\k) = iy["'\o,k) = e'^'^'^Tr/n^ 

k G Zn- 

If a group operation T acts on Zn by inversion: 
TkT~^ = —k mod n, k G Zn, then T act on an 1-cocycle 
iVmik) in 'Hi[Z„, Z„ as T • ^^(/c) = a;„i(-fc) = 

w-m mod «(fc). Since n^[Zn, U{1)] = H^Zn, Z] , wc find 
that 

T-a^~a, aen^[Zn,Z]. (J23) 

A similar result can also be obtained for U{1) group: 

T-a^~a, aen^[Z^,Z]. (J24) 

Such a result will be useful later. 

We can also use the above approach to calculate some 
other cohomology groups. To calculate C/t(1)], 
we note that the invariant submodule [Ut{1)]^^ = Z2, 
and the map iV^r becomes a — > 1. So Img[7V2T, [/t(1)] = 
Zi and KeT[NzT,UTil)] = Uril)- Also iItUt{1) = 
Ut{1). Thus 

U^[ZlUril)]^\l^ ^r^"d^' 
1/1 II d = 1 mod 2, 

which reproduces the result obtained in the last subsec- 
tion. 

Now let us calculate H'^iZ^ , Zt], where = {E,T} 
has a non-trivial action on the integer module Z^: 



T ■ n = — n, E ■ n = n, n G 



(J26) 



(J20) 



We note that the invariant submodule [Zy]^^ = Zi, and 
the map Nzt becomes n — )■ 0. So Img[iV2r,ZT] — Zi 
and Kei[NzT, Zt] = ^t- Also /^tZt = 2Zt. Thus 

[1 if d = 0, 

-^"[^2^ ,Zt] = <( Zi ifd:=:0mod2, d>0 (J27) 
^2 if d = 1 mod 2. 
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Next, let us consider H'^[Z2 x Zp,U{l)] where p is 
an odd number. Notice that Z2 x Zp = Z2p = 
{1, zt, z"^ , z^t, zP~^ ,t, z, ...} where t generates Z2, z 
generates Zp, and zt = tz generates Z2p. So we have 



C/(l) 



H''[Z2 X Zp,U{l)] = { Z 



1 

2 X 



if d = 0, 

if d = mod 2, 

if d = 1 mod 2. 



d > 



(J28) 



Last, we consider H'^[Z|J„ C/t(1)]. Z^p = {1, z, z'^ , ...} 
where contains a time-reversal operation and z^" 

contains no time-reversal operation. Thus s(z") = 

acts non-trivially on [/t(1): z" : a a^^^^"!, 

a e ?7t(1). So we have [Uril)]^^" = Z2, and the map 
N^T becomes a — >■ 1. So lmg[NzT ,UT{i)] = Zi and 

Kev[NzT,UTil)] = Uril). Also /^t C/t(1) = C/t(1). 
Thus 



I2 if d = mod 2, 
[1 if = 1 mod 2. 



(J29) 



When p is odd, = x Zp, and we have 



4. Some useful tools in group cohomology 

To calculate more complicated group cohomology, such 
as 'H'^[Zm X Z„, t/(l)], we would hke to introduce some 
mathematical tools here. 



a. Cohomology on continuous groups 

In the above discussion of group cohomology, we have 
assumed that the symmetry group G is finite. For con- 
tinuous group, one can also define group cohomology. 
One may naively expected that, for continuous group, 
the cochain Vd{{gi\) should be a continuous function of 
QiS in G. Such a choice of cochain indeed give us a defi- 
nition of group cohomology for continuous groups, which 
is denoted as 'Hi[G,U{l)].^^'^'^ 

However, continuous cochain is not the right choice 
and 'Hi[G, is not the right type of group cohomol- 

ogy. Although 'H\[G,U{1)\ does classify all the ID rep- 
resentations of G, 'H1[G, U{1)] only classifies a subset of 
projective representations.'^'*'^''' In fact, 'Hl[G,U{l)\ only 



classified topologically split group extensions 
1 ^ C/(l) G-> 1 



(J31) 



such that, as a space, E = C/(l) x G.^"* However, a generic 
projective representation of G can be viewed as an C/(l) 



extension of G where the extension, as a space, can be a 
principal U{1) bundle over G. 

So we need to come up with a generalized definition of 
group cohomology, such that the resulting H^[G, ?7(1)] 
classifies the projective representations of G. In fact, 
there are many different generalized definitions of group 
cohomology for continuous groups. ^^"^^ What is the right 
definition? We note that the cochain Vd{{gi}) is re- 
lated to the action-amplitude e"'^^^^'^) (see eqn. (40)) 
that describe our physics system. In general, the action- 
amplitude e~'^*^^^*^^ is a continuous function of gi. How- 
ever, the cochain J^d{{gi}) is actually a fixed-point action- 
amplitude which is a limit of the usual continuous action- 
amplitude away from the fixed point. So the fixed-point 
action-amplitude, and hence the cochain i^(i({5i})i may 
not be continuous function of gi. For example, as a limit 
of continuous functions, it can be piecewise continuous 
fimction. 

We also note that the cochain appears in the path in- 
tegral. So only the integration values of the cochain over 
sub-regions of G are physical. Two cochains are regarded 
as the same if their integrations over any sub-regions of 
G are the same. 

The above considerations suggest that the proper 
choice of the cochain i^d{{gi}) is that the cochains should 
be measurable functions. Measurable functions are 
more general than continuous functions, which can be 
roughly viewed as piecewise continuous functions. Such 
a choice of cochain defines a group cohomology called 
Borel cohomology.^^ We will use ■H^[G, C/(l)] to denote 
such a group cohomology. The SPT phases with a con- 
tinuous symmetry are classified by the Borel cohomol- 
ogy group 'H^[G, U{1)\. It has be shown that the second 
Borel cohomology group H^[G, U{\)\ classifies the projec- 
tive representations of G,^^ which classifies the ID SPT 
phase with an on-site symmetry G. 

On page 16 of Ref. 95, it is shown that 'H%{G,R) = 
ni{G,R) (there 'H%{G,R) is denoted as H^^^JG, R) 
and Hi{G,R) as n'^ioh,ciG,^))- It is mentioned that 
'H^(G, R) = Zi. It is also shown on page 16 that 
^foc,c(G,;7(l)) = nUG,U{l)) and ?^f,,,,(G, [/(I)) = 
ffd+i^j^Q^ Z) (where G can have a non-trivial action on 
Z and H'^^^{BG, Z) is the usual topological cohomology 
on the classifying space BG of G). Here Hj^j. ^.(G, U{1)) 
is the group cohomology where the cochain Vdigo, ■■■Tgd) 
is a continuous function near gi — 1. Therefore, we have 

ni{G,U{l)) = H'^+\BG,Z), Hi(G,R) = Zi, d>0. 

(J32) 



These results will be useful later. These results are valid 
for both continuous groups and discrete groups. Also G 
may have a non-trivial action on the modules U{1) and 
Z. In this paper, we use 'H'^{G,M) to denote the Borel 
group cohomology class H%(G, M). 
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ABC 



FIG. 41: (Color online) The graphic representation of A ^ B 
and A ^ B ^ C. 













AbT 







ABC 



FIG. 42: (Color online) The graphic representation of 1 — >■ 
A ^ B ^ C ^ 1, where Ab = 1, A ~ B^ = Be, and 
= C. 

b. Relation between group cohomology classes with different 
modules 

Let A, B, C be G-modules related by an exact se- 
quence: 

(J33) 

Then there is a long exact sequence in cohomology: 

-^n"{G, A) n°iG, B) n°{G, c) 

n\G, A) n\G, B) n^G, G) ... (J34) 
We also have 

... '^n-\G,A) ^n-\G,B) ^n-\G,c) 
n"{G, A) n"{G, B) n°{G, c) ^ 

n\G, A) n\G, B) n\G, G) ... (J35) 

Here represents the trivial module Zi with only one 
elements and 

H°(G, A)^A^ = {a\g ■ a ^ a, g e G, a e A}. (J36) 

An arrow A ^ B means that A maps into a submodule 
in B, B^ C B, where a submodule in A, Ag, maps into 
the identity 1 £ B. In other words, B"^ is the image 
of the map and is the kernel of the map. Those 
maps preserve the operations on the modules. An exact 
sequence A ^ B ^ G means that B^ = Be (see Fig. 
41). SoO-^A^-B^-C— ^0, basically, is another way 
to say G = B/A. 
Since 

^ Z ^ R ^ U{1) 0, (J37) 



we have 

... ^•H°(G, Z) ^ •H''(G, R) ^ ?i°(G, U{1)) 

n^G, z) ^ n^G, R) n\G, u{i)) 

n^{G, z) ^ n^{G, R) n\G, u{i)) ... (J38) 

Since 'H'^{G,R) = Zi [see eqn. (.132)] when G is a finite 
group or a compact Lie group, we have 

0^-H'*(G,C/(l))^-H'^+i(G,Z)^0 (J39) 

So [also using eqn. (.132)] 

W^iG, U{1)) = W^+^G, Z) = H'^+\BG, Z) (J40) 

for any finite group G, and for any compact Lie group 
G (if d > 1). This allows us to use topological-space 
cohomology H'^~^^{BG, Z) to calculate group cohomology 
H''(G,t/(l)). 

We may also assume that G has a non-trivial action on 
the modules IR and Z (which are renamed Rt and Zt): 

g ■ X — s{g)x, X £ Rt, (J41) 
where 5(17) = ±1. We have 

Zt ^ Rt ^ Rt/Zt ^ 0, (J42) 

and 

... ^^"(G, Zt) ^ ^"(G, Rt) n"{G, Ut{1)) ^ 

n\G, Zt) ^ n\G, Rt) ^ n^G, Ut{1)) ^ 

ii^iG, Zt) ^ ii^G, Rt) ^ n^{G, Ut{1)) -> ... 

(J43) 

where we have used Rt/Zt = Ut{1). Since n'^{G, Rt) = 
Zi when G is a finite group or a compact Lie group, we 
have 

0^?i'*(G,[/T(l))^-H'*+^(G,ZT)^0 (J44) 

So -H''(G,[/t(1)) = n'^+\G,ZT) for any finite group or 
compact Lie group G. (On page 35 of Ref. 96, it is 
stated that ?^'^(G, ;7(1)) = n''+^{G,Z) = H'^+\BG,Z) 
for compact Lie group G.) 

c. Module and G-module 

In the next subsection, we are going to describe 
Kiinneth formula for group cohomology. As a prepara- 
tion for the description, we will discuss the concepts of 
module and G-module here. 

The concept of module is a generalization of the notion 
of vector space. Since two vectors can add, two elements 
in a module M , a,h £ M , also support an additive + 
operation: 

a + h£M. (J45) 
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The + operation commute and has inverse. So + is an 
Abehan group multipUcation. The module M equipped 
with the + operation is an Abehan group. 

Vector spaces also have a scaler product operation, so 
do modules. The coefficients n that we can multiply to 
elements a in a module form a ring R. We will use * to 
describe the scaler product operation: n * a € Af , n Cz R 
and a € M. 

A ring i? is a set equipped with two binary operations: 
addition + : Rx R —?' R and multiplication •: i?x i? — > _R, 
where • may not have a inverse. A ring becomes a field 
if • does have an inverse, except for the additive identity 
"0". 

The scaler product * satisfies, for rt, m G i? and a, 6 G 
M, 

n* (a + b) — {n* a) + (n * 6), 
(n + m) * a = {n* a) + (m * a) , 
(n • m) * a — n* (m * a), (J46) 
^B * o. — a, if i? has multiplicative identity 1r. 

We will call the structure defined by (M, R, + , *) a. mod- 
ule M over R. 

A module over R is i?-free if the module has a ba- 
sis (a linearly independent generating set): there ex- 
ist elements xi,X2,--- G M, such that for every ele- 
ment a G M, there is a unique set rii € R such that 
a = {ni* xi) + (^2 * X2) + 

Here are some example of modules. A ring i? is a 
module over itself. Another simple example is Z over Z. 
The module Z over Z is a free module. The basis set 
contains only one element 1. The third example is the 
module Z2 over Z. Such a module is not free, since if we 
choose 1 as the basis, the element G Z2 can have several 
expressions: = 0x1 mod 2 = 2x1 mod 2. However, 
the module Z2 over Z2 is a free module with basis 1. 

The module that we are going to use in this paper is 
formed by pure complex phases M = U{1). It is a module 
over Z. The + and * operations are defined as 

a + b = ab, n*a = a", a, 5 G C/(l), n G Z. (J47) 

Just like vector spaces, we can define direct sum of two 
modules. Mi over R and M2 over R, which produces a 
third module M3 over R. As a space, is given by 
M3 = Ml X M2. The + and * operations on M3 = 
Ml X M2 are given by 

(01,02) + (61, 62) = (fli + ^1,02 + 62), 
n * (fli, 02) = (n * fli, n * 02), 
ne R, ai,bi G Mi, a2,62 e M2. (J48) 

The resulting module is denoted as Mi ® M2. However, 
in this paper (see Table I) , we will use Mi x M2 to denote 
Ml® M2. Say, for module Z„ over Z, Z„ x Z„ = Z„®Z„ 
and Z^ = Z„ X Z„ = Z„ ® Z„, etc . 

We can also define tensor product which maps two 
modules. Mi over R and M2 over R, to a third module 



M3 over R: M3 = Mi M2. If the two modules Mi 
and M2 are free with basis {xi} and {yi} respectively, 
then their tensor product M3 — Mi (g)/{ M2 is simply a 
module over R with {xi (g) yj} as basis. If the modules 
Ml and/or M2 are not free, then their tensor product is 
more complicated: M3 — Mi M2 is module over R 
whose elements have the form 

[ni * (ai bi)] + [n2 * (02 <»r 62)] + 

Ui^R, a,£Mi, 6, GM2, (J49) 

subject to the following reduction relation 

((ai + 02) ®R b) - (ai (8)_R b) - (02 <»r b) = 0, 

(a (^i? (61+62)) -{a®Rbi) - (a (8)fl 62)) = 0, 

n* {a®Rb) ^ {{n * a) ®r b) = (a ®r {n * 6)). (J50) 

Such a definition allows us to obtain the following result: 

Z (K)z M = M (8)z Z = M, 

Z„ (8)z M = M (8)z Z„ = M/nM, 

1-m ®T_ Z„ = Z(„i_„), 

{A X B) <g)R M = {A<g)R M) X {B ®r M), 

M ®r{Ax B) ^ (M ®R A) X (M ®R B); (J51) 

where (m, n) is the greatest common divisor of m and n. 
In the above x really represents ®. 

We can also define torsion product Torf (,), which 
maps two modules, Mi over R and M2 over R, to a third 
module M3 over R: M3 = Torf(Mi,M2). We wiU not 
discuss the definition of the torsion product. We just list 
some simple results here: 

Tori (Z, M) = Torf (M, Z) = Zi, 

Torf (Z„, M) = Tor^ (M, Z„) = M/nM, 

Torf (Z„j,Z„) = Z(„^„), 

Torf (A X B,M) = Torf (A,M) x Torf (S,M), 

Torf (M, Ax B)= Torf (M, A) x Torf (M, B). (J52) 

Again x really represents ®. 

A G-module over i? is a module over R that also admits 
a group G action: g • a G M for a G M and g G G. The 
group action is compatible with the + and * operations: 

g ■ (a + b) ^ {g ■ a) + {g ■ b), g ■ {n* a) ^ n* {g ■ a). 

(J53) 

In the group cohomology H'^{G,M), M is a G-module 
over R. In fact, W^iG, M) is also a G-module over R. 

d. Kiinneth formula for group cohomology 

Now, we are ready to describe the Kiinneth formula for 
group cohomology. Let M (resp. M') be an arbitrary G- 
module (resp. G'-module) over a principal ideal domain 
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R. We also assume that either M or M' is i?-free. Then 
we have a Kiinneth formula for group cohomology^'^^'^'^^ 

d 

o^Yi '^^^^^ ^) ®R n'^-p{G', M') 

d+1 

WTor^{W{G,M),H'^~P+\G',M')) (J54) 

p=0 

If both M and M' are i?-free, then the sequence splits 
and we have 

n'^iG X G',M(S)B.M') 

d 

= \Y{W{G,M)®rH'^-p{G\M') X 

d+l 

[ W Torf (HP (G, M),H'^-P+^ (C, Af'))] • (J55) 

If i? is a field if, we have 

n'^{GxG',M(E)K M') 

d 

= [ Yl nP{G, M) ®K H'^-P(G', M')] . (J56) 

p=0 

For the cases studied in this paper, we have R — AI — 
Z {ie G acts trivially on Z) and M' = J.t (ie G' may 
act non-trivially on Z^). So M A/' = Z^, on which 
G acts trivially and G' may act non-trivially. Also the 
sequence splits. 

e. Cup product for group cohomology 

Consider two cochains f„j S C"i(G, Mi) and I'nj £ 
C"'^(G, M2). From Vn-^ and Vn2, we can construct a third 
cochain i/„i+„2 S C"i+"2(G,Mi (g)z M2): 

) (J57) 

The above mapping C"i(G, Afi) x C"=^(G,Af2) ^ 
C"i+"2(G, Az/i (g)z Afz) is called the cup product and is 
denoted as 

^'ni+n2 = ^ni U Vn^. (J58) 

The cup product has the following nice property 
where 



{di'n){go, ...,.g„+i) 

+ z"=0^[(-)'*^^n(50,-,5t,-,5n+l)]- 



where the sequence go, gi, gn+i is the sequence 
50,-. -,571+1 with gi removed. So if Vm and are co- 
cycles i^„j U i'„2 is also a cocycle. Thus the cup product 
defines the mapping 

H"KG, Afi) x n"HG,M2) H"i+"^(G, A/fi (»z M2). 

(J61) 

As a map between classes of cocycles, the cup product 
satisfies^"^ 

{Vm U I/„J U t'„3 = U {Vn2 U I'nj). (J62) 



Let 



H*(G,M) = H°(G,M)®H^{G,M)> 



(J63) 



ii* (G, A4^) has an additive operation + inherited from the 
module M. The cup product provided an multiplicative 
operation on H*{G,M). So H*{G,M) with the additive 
operation + and the multiplicative operation U is a ring 
which is called the group cohomology ring for the group 
G. 



5. Group cohomology of Zm x Z„ 

Let us first calculate 'H'^ [Z„i x Z„, Z]. Using eqn. (J55), 
we find that 

n'^iz^ X z) = (j[n\z^; z) 0z n^'-^z,,; Z) j 

X ( nTorf(H''(Z„;Z),H^+i-f(Z„;Z)) ) . (J64) 



\p=o 



The above can be calculated using eqn. (J21) and the 
simple properties of the tensor product <^z and torsion 
product Torf in eqn. (J5I) and eqn. (J52). 
Thus 



H"{Zrn X Z„;Z) 
^ (Z (g)z Z) X Torf (Zi, Z) x Torf (Z, Zi) 



(J65) 



•Hi(Z„ X Z„;Z) 

(Zi ®z Z) X (Z (8)z Zi) X Torf (Z, Z„) 

X Torf(Zi,Zi) X Torf(Z™,Z) 
Zi. (J66) 



H2(Z„ X Z„,. 



(J60) 



(Z„, ®z Z) X (Z (»z Z„) X Torf (Z, Zi) x Torf (Zi, Z„) 

X Torf(Z™,Zi) X Torf(Zi,Z) 
- X Z„. (J67) 



-rn ^ "-n- 
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Torf (Z, Z„) X Torf (Z„., Z„) x Torf (Z„, Z) 

2(m,n)- (J68) 



H4(Z™X^„;Z) 

(Z(8)z Z„) X (Z„ (»z Z„) X (] 

2(m,ri) X Z„i X Z„. 

(Z„ X Z„;Z) 



(J69) 



^ Torf (Z™, Z„) X Toif (Z™, Z„) 

= Z(„j_„) X Z(„ „). (J70) 

Using the relation H''(Z,„ x Z„;C/(1)) = ^''+^2™ x 
Z„; Z), d > 0, we find 



-C/(l) 

H^[z„xz„,[/(i)] = <'z;'/f_„) 



d = even 

^™ X Z„ X Z^'^^^y^ d = odd. 

m a (m,n) 



(J71) 



This agrees with eqn. (J28). 



6. Group cohomology of Z2 x Z„ 

Let us first calculate T-L'^[Zj x Z„, Zt] where acts 
on Zt non-trivially (see eqn. (J26)). Using eqn. (J55), 
we find that 

n'^iz^ X z„; Zt) = (^n^'(^2''' ^"""X^"; 

X (nTorf(H''(Zj;ZT),H'^+^-P(Z„;Z)) I . (J72) 



\p=0 



Thus 



H°(Z^ X Z„;Zt) = (Zi®zZ) = Zi. (J73) 
^^(Z^ X Z„;Zt) - (Z2®zZ) (J74) 
n^Z^ X Z„;Zt) =Torf(Z2,Z„) = Z(2,„). (J75) 
H'(Zj X Z„; Zt) = (Z2 ®z Zn) (Z2 ®z Z) 



^(2,») X Z2. 



(J76) 



n^Z'^ X Z„;Zt) = Torf(Z2,Z„) x Torf(Z2,Z„) 

(J77) 



Using the relation n''{Z'^ x Z„;Ut{1)) = n''+\Z^ x 
Z„; Zt), d > 0, we find 

, rr, (ZDxZfl^-, d~ even 

H'^[Z2^xZ„,C/t(1)]= L+n/2^'") , (J78) 



y{d+l)/2 
-(2,n) 



d = odd. 



7. Group cohomology of D2h 

The group I?2ft is the same as Z2 x Z2 x Zj. So we can 
use the Kiinneth formula and the results in subsection 
J 5 to calculate 'H'*[I?2/i, t^T(l)]- Let us first calculate 
-W^iZm X Z„ X Z2^,Zt]: 

n°[Z^ X Z„ X Z^, Zt] = Zi. (J79) 
H^[Z„, X Z„ X Zj, Zt] ^ (Z ®z Z2) = Z2. (J80) 

n^[Z,n X Z„ X Z2^,Zt] 

^Torf(Z,„ X Z„,Z2) 

= Torf (Z„„ Z2) X Torf (Z„,Z2) 



-(2,m) X Z(2,„). 



(J81) 



H''[Z„, X Z„ X zJ,Zt] 

(Z ®z Z2) X [(Z„, X Z„) 0Z Z2] X Torf (Z(„,„),Z2) 
Z2 X Z(2,™) X Z(2^„) X Z(„^„). (J82) 

n^[Z.n, X Z„ X Z2^,Zt] 

= (Z(™,„) (8)z Z2) X Torf (Z™ x Z„, Z2) 
X Torf(Z„ X Z„ X Z(„,,„),Z2) 

= Z(2,rn,„)X42.™)^4,«)- (J^S) 

n^[Zm X Z„ X Z2^,Zt] 
= (Z ®z Z2) X [(Zrn X Z„) «)z Z2] 

X [(Z™ X Z„ X Z(„^„)) 0Z Z2] 

X Torf(Z(,„,„),Z2) X Torf(Z2„ „),Z2) 
= Z2XZ^2_„.„. xZ2 xZf,„.. (J84) 



This gives us 



H'[Z,n X Z„ X Z2^,C/t(1)] 



-(2.m) X Z(2,„). 



n^[Z„, X Z„ X Z2^,C/t(1)] 
^ Z2 X Z(2,™) X Z(2,„) X Z 

X Z„ X Z2^,;7t(1)] 

— ^(2,m,n) X ^(2,m) ^ ^(2,ri)- 

n^[Zm X Z„ X Z2^,t/T(l)] 

Z2 X Z^2 X Z^2,™) X Z^2,„)- 



(J85) 



(J86) 



(J87) 



(J88) 



When n is odd, this agrees with eqn. (J30). 



The resuhs H^[Z2 x Z2 x Z^ , Ut{1)] = H\D2h, Ut{1)] = 
Zl and U^\Z2 X Z2 X Z^ , [/t(1)] = U^D^h, Vt{X)\ - Z| 
agrees with those in Ref. 53,77 obtained through direct 
calculations. 
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8. Group cohomology of U{1) 

To calculate H''[C/(1), C/(l)] (the Borel cohomology) di- 
rectly from the algebraic definition is very tricky since 
U{1) has infinite uncountable many elements. Here, we 
will use a physical argument to calculate it by first cal- 
culating H'^[Zn,U{l)], and then let n oo. This way, 
we find 



n'[uii),u{i)] = 



U{1) 



if = 0, 

if = mod 2, d > 
if d = 1 mod 2. 

(J89) 



In Ref. 105, it is stated that 



Since W^iZ^ x U{1),Rt] = ^ 
C/(1),C/t(1)]=H'^+M^I^xC/(1), 



we have 'H'^iZ^ x 
1, and 



H''[ZJx[/(1),C/t(1)]- 




d = even, 
d = odd. 



(J92) 



This can be obtained from eqn. (J78) by taking n — 
CXI and choosing lim„_j.oo ^(2,n) = ^2 for d = even, and 
lim„^oo 2(2, „) = Zi for d = odd. 



10. Group cohomology of 17(1) xi 



if d = mod 2, 
if d = 1 mod 2. 



(J90) 



This is consistent with eqn. (J89) since n'^[U{l), U{1)] = 
H'^+i[C/(l),Z]. 

We note that the ID representations of U{1), M{Ue) — 
e"'^, where the U{\) group elements are denoted as Ue^ 
are labeled by n € Z. Also, U{\) has no non-trivial pro- 
jective representation. This is consistent with the above 
results: V}[U{1),Ut[1)] = Z and 'H^[U{1),Ut{1)] = Zi. 



9. Group cohomology of Z2 x U{1) 

As pointed out in the section C, a spin system with 
time reversal and U{1) (say generated by Sz) symmetries 
has a symmetry group C/(l) x Zj. To find the SPT phases 
of such a bosonic system, we need to calculate 'H'^\U{1) x 
Z2 ,Ut{^)\- We simply need to repeat the calculation 
of 'H'^[Z'^ X Z„,Zt] by replacing Z„ with Z. We note 
that Z2 ®j_ Z„ = Z(2,„) and Tor^ (Z2, Z„) = Z(2,„), while 
Z2 (g)z Z = Z2 and TorJ(Z2, Z) = Zi. So we find 



As pointed out in the section C, a bosonic system with 
time reversal symmetry and boson number conservation 
has a symmetry group U{\) xi Zj. To find the SPT 
phases of such a bosonic system, we need to calculate 
'H'^\U{1) X Z^, [/^(l)]. In the last subsection, we calcu- 
lated W^lUil) X Zj, C/r(l)]- In those two cases, has 
a non-trivial action on the module M — Ut{1)- 



The other two related cohomology groups are 
-W^lUil) X Z2,U{1)] and H'^[[/(l) x Z2,U{1)], where the 
Z2 is a usual unitary symmetry. In those two cases, Z2 
has a trivial action on the module M = U{1). 

In this section, we will use spectral sequence method 
to calculate the above four cohomology groups from the 
following facts: 



H^P[Uil),Z] = Z, 
H2p[z„,Z] =Z„ 



H2P+1[Z2,Z] =0, 

H2p+i[Zf,ZT] =Z2. (J93) 



W^iZj X [/(1),Zt] 



d = 0, 
d = even 

d — odd. 



(J91) 



^^2 I 



Let G be one of the [/(I) x Z2, U{1) x Z2, J7(l) x 
and J7(l) X Zj. From the exact sequence of the groups 
1 ^ C/(l) ^ G ^ Z2 
H''[Z2,H9[[/(i)^M]] page: 



1, we have the following £'2''' 



n°[Z2,H^[U{l),M]] n^[Z2,H^[U{l),M]] 

■H°[Z2,n^[u{i),M]] n^[Z2,n^[uii),M]] 

n°[Z2,n^[u{i),M]] n^[Z2,H^[u{i),M]] 
n°[Z2,n°[u{i), M]] [Z2, h^[u{i), m]] 



H^[Z2,n^[U{l),M]] H^[Z2,H^[U{l),M]] 

n^[Z2,n'[u{i),M]] n'[Z2,n'[u{i),M]] 
n^[Z2,n'[uii),M]] n^[Z2,n'[u{i),M]] 
n^[Z2,n"[u{i),M]] n^[Z2.n''[uii),M]] 



To calculate 'Hp[Z2,'H''[U{1), M]], we need to know how First we consider U (1) x Z2 group and module M = Z. 
Z2 acts on H'^[U{l),M] through how Z2 acts on C/(l) In this case, Z2 acts on f7(l) group trivially and it acts on 
group and M module. M trivially. As a result, Z2 acts on 'H'^[U{1), Z] trivially: 
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T-a^a,aG W[U(1),Zt]. Note that Z] = Z, 

n'^[l2, Z] = Z2 for d = even and ■^''[[/(l), Zt] = Zi = 0, 
n'^[l2, Z] = Zi = for d = odd. We obtain the following 
-Ef'' ~ 'HP[Z2, 'H'^\U{1),J.]] page in the spectral sequence: 

00000000 
Z Z2 Z2 Z2 
00000000 
Z Z2 Z2 Z2 
00000000 
Z Z2 Z2 Z2 
00000000 
Z Z2 Z2 Z2 

In the £'2''' page, we have spectral sequence 



(J95) 



(J96) 



-a, a e W[U{1),Zt]- We obtain the following E^''^ 
'HP[Z2,'H'^[U{1),Zt]] page in the spectral sequence: 



(JlOl) 



Again, the £"2'*^ page stabilizes: E^^"^ — iJj''' If we assume 
that 7i'^[U{l) X ZJ,Zt] contains only 2-cyclic groups, 
then we have 



n 


n 


n 


n 


n 




n 

yJ 


n 















Z2 





Z2 





























Z2 





Z2 





Z2 





Z2 





























Z2 










Z2 





Z2 





























Z2 





Z2 





Z2 





Z2 



-1, 

d+l 



d 



mod 2, 

1 mod 2, 



(J102) 



generated by ^2''^ 



= 0. So the cohomology of dj'' produces the 



such that which gives 



£f'' page 



kerdf'Vimd^^''''^^ We note that all 
the non trivial terms in the £2 page are connected to an 
incoming zero and an outing zero. Thus ^2'^ — 0, which 



leads to kerdf''^ — '"^ and im^ 



In the page, we have a spectral sequence 



0. So EP'^ = Ef 



(J97) 



and again all the non trivial terms in the E^''^ page are 
connected to an incoming zero and an outing zero and 
= 0. So El'^ = El'\= El'"). This way, we can 
show that the El^'" page stabilizes: E^^ = E^'". Using 
£P^9 with p + q = n,we can calculate W[U(1) x Z2, Z], 
since [[/(!) x Z2, Z] has a filtration 

= i/^-^+i c i?;'... c i/j^ c HS = 7{"[c/(i) X Z2, z], 

(J98) 

such that H'p'/H^^^ = EP^'^-p. If we assume that 
'H'^[U{1) X Z2,Z] contains only 2-cyclic groups, then we 
have 



W'iUil) X Z2,; 



which gives 



yd/2 
-2 ! 



d mod 2, 
d = 1 mod 2, 



(J99) 



n^pii) X Z2,t/(i)] = 




d = 0, 

d = mod 2, 

d = 1 mod 2. 

(JlOO) 



Next we consider U{1) x group and module M — 
Zt- In this case, acts on U{1) group trivially and 
it acts on M non-trivially: T ■ n ^ —n, n G Z^- As 
a result, Z2 acts on H"[U{1), Zt] non-trivially: T • a — 



n''[uii)xzi,UTii)] 




mod 2, 

1 mod 2. 

(J103) 



This agrees with eqn. (J92). 

Third, we consider U{1) x Z2 group and module M = 
Z. In this case, Z2 acts on U{1) group non-trivially 
TUgT = U-0 and it acts on M trivially. As a result, 
Z2 acts on W[U{1),Z] non-trivially: T ■ c 



a e 'H"[U{\), Z]. To obtain the above result, we note T 
act on the g-cocycles in T-L''\U{1), Z] in the following way: 



T : yq{go, ...,gq) VqiTgoT 



(J104) 



Through some explicit calculations, we find that a 
T transformed 2-cocycle, i'2{TgoT~^ ,TgiT-\Tg2T-^), 
is same as — 1^2(50 j5i 152) up to a 2-coboundary (see 
eqn. (.124)). Thus T • a2 -> -a2 for ^2 e n^[U{l),Z]. 
The generator ck2p G H'^p[U{1),Z] can be obtained 
from the generator a2 G 'H^[U{1),Z] by taking the cup 
product^^ a2p = Q2 U q;2 U ... U 02- For example, the cup 
product of two 2-cocycles, a2 and a2, gives rise to a 4- 
cocycle 04: a4{go, gi, gi) = (q;2 U a2)(go, 517 54) = 
a2(go,5i:52)a2(52,53,54)- Therefore T • a2p = (-)Pa2p. 

We obtain the following E^'" = ^^[^2, -H«[?7(l), Z]] 
page in the spectral sequence: 



(J105) 



But now, we can no longer show that the E2'" page sta- 
bilizes. 
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51 



So, we can only obtain a weaker result 



n'^[u{i) X Z2,z] < < 



Z X Z|, 




mod 4, 


d-l 

7 ^ 


d = 


1 mod 4, 


d+2 




2 mod 4, 


d+1 

7 ■* 




3 mod 4. 



(J106) 



(In fact, we can show E^^'^ = E^''' when p + q <2. So, < 
becomes = for d = 0,1, 2.) The above gives 



(U{1) 



W^Pil) X Z2,Uil)] < < 



-2 ' 

d + 3 

^2^ : 

d + 2 



d + 1 

-2 : 



d = 0, 

d = mod 4, 
d = 1 mod 4. 
d = 2 mod 4, 
d = 3 mod 4, 



(J107) 



Last we consider x group and module M — 
Zt- In this case, acts on U{1) group non-trivially 
TUqT — U-g and it acts on M non-trivially: T ■ n ^ 
—n, n e Zt- As a result, acts on H'^[?7(l), Zt] non- 
trivially: T-a ^ a e -H9[C/(1), Zt]. We 
obtain the following E^''^ = W[Z2,W[U{1),Zt]] page 
in the spectral sequence: 





Z 



Z2 



Z 



Z2 





Z2 





Z2 






(J108) 



Again, the above E2 page may not stabilize. 
So we have 



-2 ! 

d+3 



: x . 

d+1 



d-2 
7 4 

-2 ■ 



d — mod 4, 

d = 1 mod 4, 

d = 2 mod 4, 

d — 3 mod 4, 

(J109) 

d = mod 4, 

d = 1 mod 4, 

d = 2 mod 4, 

d = 3 mod 4. 

(JllO) 

We note that [G, Ut{1)] classifies the ID representa- 
tion and 7i^[G, Ut{1)] classifies the projective representa- 
tion of G. The ID representation and the projective rep- 
resentation for groups U{1) x Z2, U{1) x Z2, U{1) x Zj, 



X Z2^,Zt] < < 



which gives 



W^lUil) X Z{,UTil)] < < 



d+4 

Z2' : 

Zxl 

d + 2 

Z2' : 

d+1 



and U{1) X Zj are discussed in section C. They agree 
with n^[G,UT{l)] and ■^^[G, C7t(1)] calculated here. In 
particular, < becomes = in eqn. (J107) and eqn. (JllO) 
for d = 0,1,2. 

Ref. 103 gives a calculation and obtains a more com- 
plete result for H'^[f7(l) x Z2, Z] and H'*[?7(l) x Zj, Zt]: 



W^pil) X Z2,l] = < 



-1 

'-2 

U X 
-2 



d = 
d= 1, 

d = 2,3,5 (Jill) 
d = 4, 
d = 6. 



and 



•H'^[C/(1) X Z2^,Zt] 



U d 
I2 d 
Z or Z2 X Z d 
d: 





1,3,4 

2, 

5. 

(J112) 



From that we can obtain 'H'^\U(\) x Z2,U{\)] and 
•H''[t/(1) X Z2,C/t(1)]: 



H'^[C/(1) X Z2,C/(1)] = 




and 



-H^iUil) X Z2^,t/T(l)] 



d = 0, 

d= 1,2,4 (J113) 
d = 3. 



I2 d==0,2,3 
Z or Z2 X Z d = 1, 
I d^4. 

(J114) 



Since ■H^[C/(1) x ZJ,Ut{1)] classifies the ID represen- 
tation of C/(l) X Zj , from the calculation in subsection 
C, we find n^[U{l) X Z^,Ut{1)] = Z. < becomes = in 
eqn. (J107) and eqn. (JllO) for d = 0, 1, 2, 3, 4. 



11. Group cohomology of Z„ xi Z2 

In this section, we are going to use spectral se- 
quence method to calculate 'H'^[Zn x Z2, U{1)], 'H''[Zn x 
Z^,Ut{1)], n''[Zn X Z2,U{1)] and H'^[Z„ x Z2,Ut{1)]. 
The group Z„ x Z2 contain two subgroups Z„ = {Uk,k = 
0,1, - 1} and Z2 = {1,T}. We have TUkT = 
U-k mod n and = 1. Just like the U{1) x Z2 cases 
studied in last section, for those four groups, the E2 page 
of the spectral sequence do not obviously stabilize. How- 
ever, it turn out that the E2 pages do stabilize, can we 
can calculate 'H'^ directly from the E2 page. 

We need to first calculate 'H''(Z2,Z„) and 
n''{Z2,ZT,n)- To calculate -H''(Z2,Z„) using eqn. (J19), 
we note that Zf^ = Z„ and = '^i- The 

map becomes : a 2a for M — Z„. We 
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have Inig(iV22 , Z„) = 2Z„ = Z„ when n = odd and 
Z„/2 when n = even. This gives 
Kei{Nz^,I-) — Z(2^„). So we have 



\ iid = 0, 

'-(2,n) if > 



(J115) 



To calculate 'H'^(Z2, Z^.n) where Z2 acts non-trivially 
as T • a = —a, a G Z„, we note that Z^^^^ = ^{i.n) 



and 



Iz2^T.n — 2Z„. The map becomes Nz2 : a — > 0. So 
we have Img(iV22,Z„) = Zi. and Ker(iV2„,Z) = Z„. So 
we have 



H'*[Z2,Zt,„] =Z 



(2,n) 



(J116) 



First, we consider Z„ x Z2 group and module M — J.. 
In this case, Z2 acts on the Z„ subgroup trivially and it 
acts on M trivially. As a result, Z2 acts on ■H'[.Z„,Z] 
trivially. This allows us to obtain the following £2''^ = 
'HP[Z2,'H'^[Zm Z]] page in the spectral sequence: 



-Bf = TiJ'[Z2 , H'^[Zn, Zt]] page in the spectral sequence: 



























(2,n) 




Z(2,n) 




2(2,n) 




2(2, n) 




2(2,ri) 




2(2,n) 




2(2,n) 




2(2, n) 




(2,n) 




Z(2,n) 




2(2,n) 




2(2, n) 




2(2,ri) 




2(2,n) 




2(2,n) 




2(2, n) 




(2.n) 




Z(2,n) 




2(2,ji) 




2(2, ri) 




2(2,n) 




2(2,ri) 




2(2,n) 




2(2, n) 







Z2 





Z2 





Z2 





22 

(J120) 



Again, we can show E"^' = i?2'' for n = odd, but not for 
n = even. So we only have 



H'^[Z„ X Z2^,Zt] < 



which gives 



d = mod 2, 



^2XZ(2^_„), d-lmod2. 



(J121) 



























z„ z 


■(2,n) 


2(2,n) 


2(2. n) 


2(2,n) ^ 




-(2,n) 


2(2. n) 


























z„ z 


■(2,n) 


2(2,n) 


2(2. n) 


2(2,n) ^ 




-(2,n) 


2(2. n) 


























z„ z 


■(2,n) 


2(2,n) 


2(2. n) 


2(2,n) ^ 




-(2,n) 


2(2. n) 
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22 





22 





22 



















(J117) 


We can 


show EP^'J 


= EP2'' 


when n 


= odd. 


But for n = 



even, the above E2 page may not stabilize. So we only 
have 



n''[z„ X Z2,z] < 



which gives 



Z2 X Z„ X 

d-l 
7 2 

-(2,n)' 



d-2 
? 2 

-(2,n)' 



d = 0, 

d = mod 2, 

d = 1 mod 2, 
(J118) 



rt/(l), d = 0, 

n''[Zn X Z2, C/(l)] < J z|_„), d = mod 2, 

[Z2 X Z„ X Z^"^), d=lmod2. 

(J119) 

Second, we consider Z„ x group and module 
M = Zt- In this case, Z2 acts on the Z„ subgroup 
trivially and it acts on M non-trivially. As a result, 
Z2 acts on T-L'^[Zn,1j'] non-trivially: T • a — —a, 
a e 'H*[Z„,Zt]. This allows us obtain the following 



H'^[Z„ X Zi^Uril)] < 



-2 '^(2,n) 

d + l 



xZ,^, d = 0mod2, 

d = 1 mod 2, 

(J122) 



-(2,«) 



Third, we consider Z„ x Z2 group and module M — Z. 
In this case, Z2 acts on the Z„ subgroup non-trivially and 
it acts on M trivially. As a result, Z2 acts on H'^[Zn, Z] 
non-trivially: T-a -> a G H«[Z„, Z]. To obtain 

the above result, we note that T act on the g-cocycles in 
H'^[Zn, Z] in the following way: 



T : Vq{gQ, ...,gq) -> Vq{TgoT' 



,TgqT ^), G Z„. 

(J123) 



Through some explicit calculations, we find that a T 
transformed 2-cocycle, V2{T g^T'^ ,T giT'^ ,T g2T-^) , is 
same as —V2{9ot 91,92) up to a 2-coboundary. Thus 
T • a2 — —0-2 for 02 G 'H^[Z„,Z]. The generator 
ct2p G 'H^^[Z„,Z] can be obtained from the generator 
a2 G y?[Zn,J.] by taking the cup product^'^ a2p = 
a2 U a2 U ... U Q!2- Therefore T ■ a2p = (— )^ck2p- 

This allows us obtain the following £'2'^ — 
'HP[Z2,'H^[ZmZ]] page in the spectral sequence: 



























(2,n) 




2(2. n) 




2(2,n) 




2(2. n) 




2(2.n) 




2(2. n) 




2(2, n) 




2(2. n) 




z„ 




2(2,n) 




2(2,n) 




2(2, Tl) 




2(2,ri) 




2(2,n) 




2(2,n) 




2(2, n) 




(2,n) 




2(2,n) 




2(2,n) 




2(2, Tl) 




2(2,n) 




2(2,n) 




2(2,n) 




2(2, n) 




z 





22 





22 





22 




(J124) 
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We have 



even. The above gives 



'-2 X ^ 

7 2 

-(2,n)' 
l2 X A 



-(2,n)' 



■(2,n)' 



d = mod 4, 

d = 1 mod 2, 

d = 2 mod 4, 
(J125) 



which agrees with the result obtained in Ref. 106 for n = 
even. The above gives 



W'iZ,, X Z2,U{1)] < { 



d 

^(2,n)' 

Za X Z 
Z2 X Z 



d-l 
7 2 

"(2,«)' 



d = 0, 

d = mod 2, 

d = 1 mod 4, 

d — i mod 4, 
(J126) 



This agrees with a result for the symmetric group on 
three elements, 6*3 = Z3 x ^2^°'' 



U{1) ifd = 0, 

11 if d = mod 2, d > 

12 if d = 1 mod 4, 
Iq if d = 3 mod 4, 

(J127) 



if we replace < by =. 

Last, we consider Zn x group and module M = 1.t- 
In this case, Z^ acts on the Z„ subgroup non-trivially 
and it acts on M non-trivially: T ■ n ^ — n, n e Z^. 
As a result, acts on 'H'[Z„,Zt] non-trivially: T • 
a ^ a € H^iZn, Zt]. We obtain the following 

iJ^'"? = W\Z2, 'H'^[Zn, Zt]] page in the spectral sequence: 
























Zn Z(2.„) I 




(2,n) 




2(2.„) 2 



(2,n) 2 




-(2.ri) 




2(2,n) 




2(2.n) 




(2,n) 2(2. „) Z 




(2,n) 




2(2.„) 2 



(2,n) 2 




-(2.n) 




2(2,n) 




2(2.n) 




Zn 2(2. „) 2 




(2,n) 




2(2.„) 2 



(2,n) 2 




-(2.n) 




2(2,n) 




2(2,n) 




Z2 





I 


-2 





Z2 





22 

(J128) 


We obtain 
























)' 


d = 


mod 4, 


WiZn X Z^ 


,Zt] 


'1 


22 X 


d-l 

^(2,n) 


d = 


1 mod 2, 








.2„ X 


d-2 
7 ^ 

^{2,n) 


d = 


2 mod 4, 



(J129) 

which agrees with the result obtained in Ref. 106 for n = 



'-2 X 



W^lZn X Z^,Ut{1)] < < Z„ X Z 



X 

d + 1 

(2,«)' 



(2,«)' 

d-l 

(2,«)' 



d = mod 2, 

d = 1 mod 4, 

d = 3 mod 4, 
(J130) 



Although we cannot prove it, in fact, the i?f pages do 
stabilize: E^^ = E^'"^ , for all the four groups discussed 
here. For n = odd, we can show that the i?^''' pages stabi- 
lize. So we can replace < by = in eqn. (J119), eqn. (J122), 
eqn. (J126) and eqn. (J130). For n = even, the results 
obtained before using the Kiinneth theorem imply that 
we can replace < by = in eqn. (.1119) in eqn. (.1122), and 
the results in Ref. 106 imply that we can replace < by 
= in eqn. (.J126) and eqn. (.1130). This suggests that the 
-Ej''^ pages stabilize even when n = even. 



12. 



Group cohomology of U{n), SU{n), and Sp{n) 

SU{n), and Sp(n) 



The group cohomology ring of C/( 
are given by^'^'' 



n) 



■H*[U{n),] 
n*[SU{n),] 
n*[Sp{n),] 



Z[ci, c„], 
Z[C2, c„], 
2bl, •••,Pn], 



(J131) 



where c, G H'^^[U{n),Z] or c, e 'H'^'[SU{n),Z], and 
Pi e ■H**[S'p(n), Z]. Here Z[x,y, ...] represents a ring of 
polynomials of variables x, y, ... with integer coefficients. 

For example 'H*[J7(1),Z] = Z[ci] mean that the el- 
ements in 'H*[U{1),Z] has a form uq + riici + n2cf + 
nacf + ... = no + "ici + "-2C1 U ci + 713C1 U ci U ci + .... 
Note that ci is a two cocycle and rtici is a two cocycle 
labeled by rii € Z. Also niCi is the only 2 cocycle in the 
expression niCi + n2Ci U ci -I- naci U ci U ci -I- .... Thus 
■H^[J7(1), Z] = Z. Similarly n2CiUci is the only 4-cocycle 
in the expression riiCi + n2Ci U ci -I- n^Ci U Ci U ci -I- .... 
Thus -H*[;7(1),Z] = Z. There is no odd cocycles in 
niCi+n2CiUci+n3Ci\JciUci + .... Thus 'H''[t/(1), Z] Zi 
when d = odd. 

From 'H*[5[/(2),Z] = Z[c2], we find that the elements 
in H* [SU (2), Z] has a form hq + niC2 -1- ri2C2 U C2 -f ^302 U 
C2U C2 + .... Thus 



n'^[SUi2),Z] 



and 



n'^[su{2),u{i)] 



d = mod 4 
!i d 7^ mod 4. 



[ d = 3 mod 4 
-1 otherwise. 



(J132) 



(J133) 
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13. Group cohomology of 50(3) and 50(3) X Zi' Using (J55) and •H''[S'0(3) x 2'J,f/T(l)] = 

'H'^+^[SO{^) X Zj, Zt], we obtain 

The group cohomology ring of SO{i) is given by^°^'^''^ 
n* [SO{-i), Z] = Z[v, c]/(2w), (J134) 
where z; € H^[SO{i),I.] and c e ■H4[S'0(3), Z]. Here 

/(2w) means that the expression 2v in the polynomial is , 
regarded as 0. The elements in %* [SO{2,) , Z] have a form [ ^2 d = 0,3 

J^nindv^. Note that c'z;-' is a 4i+3j cocycle. The lowest ^/^^^c'/^^o^ rr J '^=1; fnoa\ 

cocycle m the expression i^nijC v-' is a 3 cocycle noiw. | Z2 d = 2, 

Since 2ti is regarded as zero, there are only two 3 cocycles d — 4^ 

and V labeled by noi = 0, 1. Thus n^[SO{3),Z] = Z2. ^ ^ 

The expression 'Y^Ui^jdv^ contains only one 4-cocycle 
nioc labeled by nm € Z. Thus -H4[S'0(3), Z] = Z. This 
way, we find that 

fz d = 0,4 

(J135) 





I Z d = 


0,4 


•H'*[S'0(3),Z] = 1 


Zi 


1,2,5 




[Z2 d = 
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